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ABSTRACT: 

The method developed by Bauer, Garabedian, and Korn [2] for the design 
of wing sections which exhibit shock-free flow at high subsonic speeds has 
been applied to the problem of cascade design. A cascade is a periodic dis- 
tribution of two dimensional blade cross sections which serves as a basis 
for the design of axial flow compressors and turbines. 

Those aspects of the design procedure which are used for both airfoil 
and cascade design are reviewed briefly in this paper. Emphasis has been 
made on the differences, namely in the form of the singular solution in the 
hodograph plane. The singular solution for flow around an airfoil consists 
of a pole and a nearby branch point in the hodograph plane. The singular 
solution for the cascade is formed by two logarithmic singularities resembling 
a source and sink. 

A listing of the computer program used to design transonic cascades 
with this method has been included. Details omitted from the text of this 
paper can be found by careful examination of the camputer program. In order 
to make it possible to mn this program a description of the input needed to 
execute this program has also been included. The output listing for the blade 
shown in Figures 2 and 3 which has been included contains the list of input 


parameters necessary to regenerate this case. 
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1. Introduction 

In recent years there has been a growing concern about the availability 
and utilization of our energy resources. Considerable effort has gone into 
the development of new resources and into methods of more efficient use of 
existing energy supplies. We will describe here a procedure which can be 
used to design more efficient high speed turbomachinery which could result 
in significant fuel conservation. In particular, this technology could be 
the basis for the design of improved axial flow compressors used for the 
production of enriched uranium at gaseous diffusion plants. Also, it can be 
used to design lighter weight and more efficient airplane engines. 

In this paper we used a technique developed by Bauer, Garabedian and 
Korn [ 2 ] which has been used to design airplane wing sections and apply it 
to the problem of designing transonic turbine and compressor blades used in 
axial flow pin Senerecnly The method is applicable to two-dimensional problems 
and therefore our model for design will be a periodic distribution of two 
dimensional blade cross-sections called a cascade. As the distance between 
blades increases our solution must approach that of the isolated airfoil and 
therefore we can view this paper as an extension of the airfoil design method. 

As with the airfoil problem, our concern is to find smooth, transonic 
steady state solutions to the equations of irrotational motion of an inviscid 
compressible fluid about some object. In this case the object will be a 
periodic distribution of blades. The equations of motion are the same as those 
used for the airfoil design problem, therefore, the same technique of using 
conjugate characteristic coordinates extended into the complex domain can be 
employed. A more detailed description of this method can be found in reference 
[ 2]. 


The main differences between designing a compressor blade and an airfoil 
are the locations and types of singularities of the solution in the 


hodograph plane. In the airfoil case, the basic singularity was a pole in 
the hodograph plane at some prescribed free stream velocity. For the 
cascade problem the velocity far upstream of the cascade, called the inlet 
velocity, is different from the velocity far downstream of the cascade, 
called the exit velocity. In fact the purpose of the cascade is to change 
the inlet velocity in some prescribed way. Thus, in the hodograph plane 
there are two velocities which correspond to infinity in the physical plane. 
The appropriate singularities in this case are logarithms at the inlet and 
exit velocities. As the distance between these logs decreases the spacing 
between blades increases approaching the isolated airfoil case. The 
logarithmic singularities behave like a dipole and, therefore, approach the 
correct singularity for flow around an airfoil in the limit. 

For the cascade design problem, as with the airfoil design problem, we 
solve the inverse problem; that is, we generate solutions to the equations of 
motion with the prescribed singularity and examine the streamline passing 
through the stagnation point. This streamline describes the shape of the body. 
To generate solutions to the flow equations we choose a set of parameters which 
define physical quantities and describe the initial analytic function. We 
also prescribe where in the hodograph plane we would like the blade to lie. 

The program finds a solution to the equations of motion whose stagnation stream 
line matches the prescribed one in a least square sense. If this curve which re-— 
presents the body is not closed or if it is not physically realizable in the phy- 
Sical plane,for example negative thickness, the input parameters must be 


changed and the program rerun. 


In Section 2 the equations of motion are presented and rewritten in 
characteristic coordinates which are extended into the complex domain. A 
formilation of a method to generate solutions with the appropriate singular- 
ity for a cascade appears in Section 3. In Section 4 we discuss the problem 
of choosing the initial function. In Section 5 a description is given of 
the numerical technique which has been programed on a CDC 6600 computer. 

We have enumerated many of the modifications of the design procedure of 
reference [ 2 ] and presented some results in Section 6. Finally, a listing 
of the FORTRAN program used to generate the results of Section 6 has been 
included, which should help explain many of the details. A sample run has 
been included along with glossaries describing input specifications. This 
should make it possible for someone to copy and execute this program. We 
intend to send a copy of this program for inclusion in the library at the 
Argonne National Laboratory. 

The author would like to thank the many people who in one way or 
another assisted me in this work. The idea of extending the design procedure 
to cascades was suggested independently by Milton Plessit and A. J. Acosta 
of the California Institute of Technology and Julian Prince of United Aircraft. 
I am deeply indebted to Paul Garabedian for the aid and encouragement which 
was given to me at all stages of development of this work. Eldon McIntyre 
has assisted in the preparation of this paper by carefully reading an early 
draft of this paper and correcting numerous errors in many of the formulas. 
Frances Bauer read a later draft and made numerous helpful suggestions. 
Farntella Graham patiently typed and retyped this manuscript. Finally, I 
would like to thank my family whose patience and understanding has not gone 


unnoticed. 


2. Equations of Motion 
Let x and y be Cartesian coordinates and let u and v be the corresponding 


velocity components. Then the steady state equations of motion for an inviscid 


polytropic gas with isentropic, irrotational flow is 


(2.1) (c? - u?)9,, — 2uv 6, + (c? - v2)¢,, = 0 


where $ is the velocity potential, oy =u, 6 =v, andc is the local sound 


y 
speed. 
Let p denote pressure and p density. The local speed of sound,c , 


is defined by 
2= dp = 2 
(2) a do Y 5 


where y is a gas constant. Conservation of energy is expressed by Bernoulli's 


equation 
(2.3) e=c2-Tig ; q? = u2 + v? 


where Cy is the sound speed at stagnation velocity. Equation (2.3) shows 
that c, hence o and p are functions only of the speed q. 

Introducing characteristic coordinates — and n, equation (2.1) can be 
written as a system of four linear equations for u, v, x, and y [cf.2 ]. 


Computing the characteristic directions 


(2.4) 4, = wt evg?oc? 


+ 


we can write the equations of motion in characteristic form as follows: 


(2.5a) v, =A,u 3 vi Sstrxhu 


g gant = n =7 
(2.5b) AVp + Xp = Ogee Y, + ad 0 


We define the local Mach number M as q/c. The equations of motion 
are hyperbolic when A + are real and unequal which occurs when the flow 
is supersonic, M>1. For M< 1 the flow is subsonic and the complex 
conjugate roots in i + mean that the equations of motion are elliptic. 

The set of points where M = 1 is called the sonic line and is a circle 
in the hodograph plane q = constant. This constant is called the critical 
speed, which we denote by dy, . 

Since we are interested in solutions to these equations in a region 
of the hodograph plane where both subsonic and supersonic flow is to occur 
we must use a method applicable to both regions. The method of complex 
extension [ 4] is ideally suited since in the complex domain there is no 
distinction of type, in fact, A + need not be real or complex conjugate. The 
only theoretical difficulty which arises in using complex extension is along 
the set of points where Ay =)_, i.e. M=1, since equations (2.5b) are 
coincident there. It has been shown that this problem is easily treated 
(er. £ 2; 9)). 

With the method of complex extension the independent variables € and n 
are each extended into the complex domain so the solution u, v, x and y mst 
be allowed to take on complex values. Since both & and n have real and 
imaginary components this amounts to enlarging the independent domain from 
two dimensions to four dimensions. We are interested in the solution only in 


the real hodograph, We would like to compute the solutio in as small a 


subset of the extended space as possible in order to obtain the solution 
in the real hodograph plane. It turns out to be possible to choose 
characteristic coordinates — and n so that the set of points — =n 
includes the real subsonic hodograph plane. We call this choice of 
characteristic coordinates, conjugate characteristic coordinates as des- 
eribed in [ 2]. 

Equations (2.5a) are independent of x and y and therefore u(&, n) 
and v(€, n) depend only on our choice of characteristic: coordinates, not 
on the solution to our equations. We may therefore view & and n as being 
generalized hodograph variables. We can actually write closed form 
expressions for conjugate characteristic coordinates in terms of the hodo- 
graph variables u and v and thus solve (2.5a) explicitly. By integrating 
(2.4) (cf. [11]) we find that s and s* are conjugate characteristic 


coordinates by choosing 


(2.6) See) eae See 


where 
(2.7) rife) nehedhi(e ae cgay net ae aac 


where wu? = (y-1)/(y+1) , 6 = tan! v/u and K is a constant chosen so that 
h(1) = 1. Note that as the critical speed qy becomes infinite h(q) approaches 
q. Thus for incompressible flow s approaches w = u-iv and it follows that 
s is an extension for compressible flow of the complex velocity. The function 
h is real when q is real and subsonic so that s and 3" are complex conjugates 


in the subsonic region. We choose coordinates & and n and consider mappings 


of the form 


(2.8) s=F(E) 3s" = F(R) 

where F is an arbitrary analytic functim. For any function F , 

—E and n are characteristic coordinates since any function of a 

characteristic coordinate is itself a characteristic coordinate. Also, 

it is clear from (2.8) that whenever s = 3" we have — =n so that & 

and n are also conjugate coordinates for any analytic fumction F. 
In previous work on airfoil design we have found it sufficient to 


choose a two parameter family of transformations 
(2.9) s = F(E) = 1 - Be - &? 


for the parameter B which is suggested by considering incompressible 

flow around an elliptic cylinder. This led to a class of solutions which 
had a single branch point in the hodograph plane at — = -B. Moreover, design 
in the —-plane is simpler than in the w-plane because a simple closed curve 
representing the body streamline is transformed into the two-sheeted hodo- 
graph. Turbine and compressor blades generally have a high degree of camber 


and may require more than one branch point. Therefore we have chosen 
(2.10) s 6 1 +-Cye + Coe? + Caer ee 


which can have as many as three branch points. Locating the branch points 
can be done by examining their location for incompressible flow around 2 


cambered Joukowski profile. 


Work is presently underway to investigate a rational method for 
the selection of the fumction F to achieve the desired design objectives. 
The selection of a suitable transformation F for turbine design will be 


the subject of a future report [10]. 


3. Singular Solution for the Cascade Problem 

We apply the method of complex extension to equations (2.5) in order 
to generate solutions which represent flow around a cascade. We will do this 
by setting up an appropriate characteristic initial value problem in the ex- 
tended hodograph plane. 

For the characteristic initial value problem in the extended variables 
we chose two characteristic planes &€ = Eo and n = No and assign initial data 
for either x or y along each of them. Equation (2.5b) is a set of ordinary 
differential equations along a characteristic which can be integrated to get 
the other independent variable. The solution at points not on the characteristic 
initial planes can be found by the numerical technique described in Section 5. 
The solution will be analytic as long as the initial data is analytic. 

For flow through a cascade we would like a solution which has an inlet 
velocity at —- and an exit velocity at + in the physical plane. Let Wp = 
Up - iv, and Wy = Uy - iv, denote the complex upstream and downstream velocities 
respectively. Let (Ea Np) and (En ; Ny) be the corresponding values of the 
characteristic coordinates which are found by solving (2.6) and (2.8) once the 
transformation function F has been prescribed. It follows that ny = Ey 
and n, = Ea since — and n are conjugate coordinates and uy , Uy » V, and vp 
are real quantities. We require that the solution in the physical plane have 
infinities at the two points (En , Ey) and (Ep - En) . In other words we have 
a source at Ep and a sink at EA in the real hodograph plane. We further require 


that the solution be periodic in the physical plane. 


Any closed loop in the hodograph plane surrounding both Ea and 5B corre 
sponds to a closed loop in the physical plane enclosing a blade. We will 
modify this enclosure requirement for our wake model of the trailing edge. 


One can verify that a solution of the form 


(3.1a) x = Re[X' log(— - &,) + X* log(E - &) + x*] 


(3.1b) y= RelY” loge ="85) + Y les(G Eee 


satisfies the above requirements for suitable choices of the regular functions 
x? , x? , x®° , y! , Y? and Y® . The assumption that X° and Y® are regular 
lead to the requirement that X’,Y? and x?, Y? be solutions to the partial 
differential equations (2.5b) with initial conditions 


(3.2a) A, yit+iggte 


1 
oO 


| 
oO 


(3.2b) ANA ea 


on & = Ea and € = ER respectively. This can be verified by substituting 
equations (3.la-b) into (2.56). From this construction we can see that x? , Y? 
and X* , Y? are essentially Riemann functions for equations (2.5b) . 

Equations (3.2a-b) can be solved with (2.5b) to obtain X’ along — = &, and 
X? along — = En - We must define xX! and X? along the corresponding character- 
istics of the other family n = Nao N= Np to obtain a well-posed characteristic 
initial problem for the determination of X! , X? , Y? and Y? everywhere. We 
can save computation if we define the data on the other family of character- 
istics by reflection, i.e. X!(E, Ny) = XI (Ey rey Gee leng Np) = totes » Oe 


etc. This is not always possible since X! , x? , Y! , Y? need not be real at 
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En 2 Ty and ER > However, since equations (2.5b) admit constants 


as solutions we can find solutions for X! , X? , Y! and Y? of the form 


(3.3a) wet 1X, 3 


Be ae, 


(3.3b) y! = y! + iY, 3 y? = y? + iy, 


where X, , Xj, Y, and Y, are real constants and x4. Oe Ws andl * same 
symmetric solutions to equations (2.5b) satisfying (3.2a-b) along the 
initial characteristics. 


The equations for X? and Y° obtained by inserting (3.la-b) into 


(2.5b) are 
(XA Y4X?) + 40. Y4K,)  (A_Y24K2) + 1O0Y. +X) 
(3.4a) ee en ies i ee ed Wa ik Bia Ek 
= E E En — E eB = 5 
(3.4p) Re oe mo 


where the regularity of the right hand side of (3.4a) is assured by the initial 
conditions (3.2a-b). In order to formate a well-posed problem for the 
evaluation of X? and Y° it is necessary to assign initial data along a pair 

of initial characteristics. We choose initial characteristics, & = Eo > 

n=MNq » SO that their intersection is a point in the real subsonic domain. 
This means that No = Eo - To avoid difficulties in the numerical calculation 
of the right hand side of (3.4a) we choose Eo safely away from both EA and 


EB - We assign 


(3.5) X* (Eas n) =G(n) 5 X°(E, ng) = G(ng) 


lj 


on the initial characteristics & = Eo and n = No respectively where G 

is an analytic function. The selection of G is the subject of the next 
section. Since the solution for X and Y uses only the real parts of xX? 
and Y* we can prescribe many combinations of initial data which will yield 


the same real parts. In particular, the symmetric set 


(3.6) X%(Eqs n) = E[G(n) + Ging] $ X%(Ey ng) = F (ng) + G(EDI 


has been used to reduce computation in the subsonic portion of the flow. 
We have the four real constants Xq a Yq ; Xp and YR to determine. 
First we consider the requirement that the solution be single-valued for 


any real closed path encircling a blade, We find 


XX x24 

(3.72) fac= Ref (= ce + apn )aE] = —2nhy + 3) 
WHi¥, YHY, 

(3.7b) p dy = Re} oa Pare = ~en(Y, + Yp) 


where the path of integration is around both & A and ER in the plane—& =n. 


Therefore , X = —Xy and Ye = -Y, if we are to have a single-valued solution. 


We further require that the stream function w~ defined by 
(3.8) dy = (pu)dy - (pv)dx 


be single-valued for any similar path. Since pu and pv are real in the real 


domain we have 


(3.9) pb = em[p pUpY y + Ppp s - PaVaXy = PaVpXs] =O. 
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Making use of the closure conditions (3.7a-b) and (3.9) we have 
(3.10) (O,Uy - Ppp) te (ony - PpVp)Xp 


From (3.la—b) we see that the solution is periodic with a period of 
X and Ys in x and y respectively. We define the stagger angle 8 and the 


gap g by 


= 7 ae ae 2 
£3,511) B=tan xr 3 & Xy + Yy 


Given the inlet and exit velocities the stagger angle is computed from (3.11). 
Therefore, it is convenient to measure the inlet and exit velocity vectors 
relative to 8. If we require that the blades be sta¢ked vertically, 8 = 1/2, 
we can prescribe only three of the four parameters for the inlet and exit 
velocity. For example we may prescribe the inlet and exit Mach numbers and 
the turning angle, and then determine the inlet and exit flow angle from (3.11) 
using 8 = 71/2. To have physical significance the gap must be measured relative 
to some suitable distance measurement such as the blade chord. Since this 

is not known until the solution has been found we have instead chosen to 
normalize the solution by requiring that the gap g have magnitude En - & 4 
As a result of this normalization the chord length will change slowly as 

EA + EB - Moreover in the limiting case the isolated blade will have the same 
singularity as the one used in our previous work on airfoil desig. 

The wake model of the trailing edge which was formulated in [ 1 ] to 
correct for boundary layer displacement can easily be incorporated into this 
cascade design procedure. This is done by requiring that the residues of 
(3.7a-b) be non-zero. In fact, they should be proportional to x and Yp 
respectively so that the blades will have a vertical gap when stacked 


vertically. The constant of proportionality is the ratio of the thickness of 


is 


the trailing edge outside the boundary layer to the gap and is chosen 
to be small. It should be noted that equation (3.10) will have to be 
modified accordingly. 

The problem of flow through a cascade of blades can therefore be 
solved by generating the solution of two characteristic initial value 
problems for the homogenous equations (2.5b), one along &— = En sn= Ey ' 
the other along € = EB sn Ep » and by solving a characteristic initial 
value problem for the inhomogeneous equations (3.4a-b) for arbitrary initial 
data chosen along some initial characteristics § = Eo sn= Pe The initial 
data for the homogenous equations is determined by solving ordinary 
differential equations along the initial characteristics reminicent of the 
Riemann functions. The four constants X,, Xp ; Yq and YR are completely 
determined by three closure conditions and normalization. The solution for 
x and y is formed from the solution of these three problems using (3.la-b) 


and has the singularities necessary for flow through a cascade. 
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4, The Selection of Initial Data 

Solutions to the equations of motion of the form (3.la-b) will have 
the appropriate singular solution for flow through a cascade but will not 
necessarily generate solutions with a physical interpretation for arbitrary 
initial data assigned along the characteristic initial planes & = Ec and 
N=Nc.- Once the transformation from the hodograph variables to 
characteristic coordinates has been selected the design problem centers on 
choosing the initial data so that the resulting design will achieve as many 
of the design objectives as possible. Many design objectives, such as 
boundary layer control, cannot be stated in precise terms within the context 
of this inviscid model. Other design objectives may be impossible to achieve 
within the class of shockless solutions. Therefore we should expect the 
problem of choosing initial data to require some degree of trial and error. 

To find the shape of the blade we generate a solution and examine the 
level curve jp = WS where ps is the value of the stream fimction at stagnation. 
In general, this shape will not resemble a blade and in fact does not have 
to be a closed curve. The shape of this curve will depend upon the initial 
data we choose. We now discuss a procedure for selecting the initial data 
which attempts to minimize the chances of obtaining extraneous solutions. To 
do this we make use of most of the ideas that have been used for the airfoil 
design problem described in [1, 2]. 

We pick the initial characteristic & = Ec away from Ea and 3 to avoid 


computational difficulties in evaluating the right hand side of (3.4a). 


The characteristic n = Eo is selected for symmetry reasons. Since we 
need only the real part of the solution to equations (3.4a-b) we can 
symmetrize the equations in a manner that saves half the computation in 
the subsonic region. Using equation (3.6) the initial data on n = Eo 
will be the complex conjugate of the initial data along & = Ea which we 
prescribe. In keeping with the procedure followed for airfoil design we 


choose the initial function G on the initial characteristic — = Ea as 


NL 
(4.1) G(n) = 2f (Ly +4 Lyo) log (n-L;. -i Ly) +En+E 


2 3 
j=l ae Feeu" 


tsq 
The constants that appear determine the resulting solution. The last three 
terms are polynomial terms while the remaining NL terms are logarithms 
located at n= L,tibsy . They can be placed anywhere inside the blade to 
simulate sources or sinks. Their influence is primarily local when their 
coefficients are small and can therefore be used to effect small changes 
without much change to the overall solution. 

Selection of the linear parameters has been automated by a least 
squares process. Instead of prescribing the linear parameters we select a 
curve or set of arcs in the hodograph plane where we would like the blade to 
lie. We choose the linear constants in (4.1) so that the solution it generates 


minimizes 
(4.2) fw y? dé 


where the integration is taken along the prescribed arcs and where the weighting 


function W can be chosen to achieve a better fit along any given arc. 
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By prescribing the automation arcs in the characteristic planes we 
are in effect prescribing the body streamline in the hodograph plane. This 
is very convenient for boundary layer control since we can select arcs for 
the least square automation which have sufficiently small pressure gradients 
in order to avoid separation. If we prescribe an arc so that it approaches 
the trailing edge along a line of constant speed we can eliminate boundary 
layer separation or delay it until the last per cent of chord in most cases. 
In addition we narrow the choice of initial data by placing constraints on 
the solution. We prescribe the curvature at the leading edge and the angle 
between the inlet velocity and the central streamline. We choose a velocity 
for the trailing edge and require the stream function to assume its stagnation 
value there. We further require that its first and optionally its second 
derivative vanish. This is our model of the Kutta-—Joukowski condition in the 
hodograph plane. 

The above constraints and the least square conditions greatly reduce the 
effort in selecting meaningful initial data. It should be noted that we have 
overprescribed the problem and therefore cannot expect good agreement between 
the prescribed arcs and the resulting streamline unless they are chosen in a 
manner consistent with the given compression and turning. For a description 
of how to locate the logarithms the reader is referred to [1, 2] since there 
appears to be little difference between compressors and airfoils in this aspect 


of the design problem. 
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5. Numerical Method of Solution 

The same numerical technique that was successfully used to design air- 
foils described in [2] has been used with minor modifications for the cascade 
design problem. In the real domain this technique is often called the method 
of characteristics or Massau's Method [ 3 ] and is essentially a predictor 
corrector scheme. The method is applied in the complex domain merely by 
performing all operations using complex arithmetic. 

In Section 3 we explained how to construct the appropriate singular 
solution. In doing so we must find the solutions for X!, Y!, x?, Y? and x3, 
Y? . The functions X', Y! and X?, Y? satisfy the equations of motion (2.5b) 
with initial data given by (3.2a-b) on the initial characteristics while xX’, 
Y? satisfy the inhomogeneous equations (3.4a-b) whose right hand side is 
computed from X!, Y? and xX?, Y* . Therefore, we will apply the method of 
Massau to each of these three sets of unknowns. The hodograph variables may 
be obtained explicitly from equations (2.8), (2.6) and (2.7). However, it is 
faster to use this numerical method to solve for them as well,using these 
equations only to generate characteristic initial data. 

The characteristic initial value problem for X! and Y! starts on initial 
characteristics En 2M ° The characteristic initial value problem for X? , 
Y? starts on initial characteristics ER > Tp and the characteristic initial 
value problem for X?, Y° uses initial characteristics Ea 2G along which it 
is necessary to know X? , Y? , X? and Y* in order to evaluate the inhomogeneous 
terms. To apply the method of Massau we must select a path along each of the 


initial characteristics, as illustrated in Figure la. Each path starts at the 
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point of intersection of the two initial characteristics. Let r and t 
be real parameters. We assign initial data along a path n = n(t) from 
Ny on the initial characteristic & = En and along another path — = E(r) 

from En on the other initial characteristic. The initial data along these 
paths determine the solution on the two-dimensional surface consisting of 

the points of intersection E(r), n(t) . To find the solution at a given 
point EP » Nip we choose a path on — = aN from n = Na ton= Np and a path 
onn = Ny form — = En to & = Ep - The solution at the point Ep » Np will 

be independent of path if the data along the initial characteristics is 
analytic. In order to be able to solve for X’, Y', x?, Y? simultaneously 

we choose each initial path on Ny SO that it starts from — = Ea passes 
through — = &, and then through — = Eo: Similarly the initial path along 
E= EA begins at n = Na and passes through "B and Ne Thus the values for 
the right hand side of (3.4a) are available as they are needed. This arrange- 
ment is illustrated in Figure lb. 

The solution throughout the four dimensional domain can be computed by 
selecting appropriate pairs of paths. We are interested in obtaining the 
solution only in the real hodograph and would like to do so with as few pairs 
as possible. It has been shown that by choosing complex conjugate sets of 
paths we can arrive at one line of the solution in the real subsonic domain 
for each set of paths [9 ]. One set of paths can be found to find a portion 
of the solution in the supersonic region bounded by two characteristics and 
the sonic line. Generally one pair of paths can be chosen so that this region 


contains the full supersonic portion of the flow outside the blade. 
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Each set of initial paths is described by the vertices of a polygonal 
are and a set of grid points superimposed in such a way that each vertex is 
a grid point and that the distance between adjacent grid points is less than 
some prescribed amount. We use the notation Usy to denote the value of 
ulE(P;)s n(t,)) where r and t are real parametric representations of the 
initial paths. Note that r = 0 and t = O are the initial curves on € = En 
and n = uy respectively. Using this notation we can write the predictor 
step of the method of Massau for the solution of equation (2.5a) at the point 


j, k as the two linear equations 


(5.14) Fue 7 Yyeagee = Ae gaye Bye - Yaad 


(5.1b) v; 4k - Vy Jen = Ds ved (i, 4k = Hy Hew 


for the predicted values a, k 
3 


of u and v. The values of u and v along j = 0 and k = O are determined 


and v, k which are first order accurate estimates 
5 


explicitly from the hodograph transformation (2.6 - 2.8). 
For the corrector step, we use the values obtained in the predictor step 


to compute first order accurate midpoint values for i, given by 


(5.2a) 1, =F Op yegige H yl ge9%y 40 
(5.2b) X= FLODs yea + Cy ae F 1 
and solve 

(5.3a) if eas se tc a (Uy 7 Wy_1 i 
(5.30) v d 


J 5k i: "3 sk=1 - te (uy 9k 7 y ae) 


for u, and v. to obtain second order accuracy. 
J ok Jk 
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A similar set of equations is solved to obtain x!, Y', x?, Y? 
and X*, Y? to second order accuracy. Equations (3.la-b) are evaluated 
at points in the real hodograph to obtain x and y to second order accuracy. 
The value of the stream function is obtained by applying a second order 
accurate difference approximation to (3.8). Points on the blade are 
identified by the streamline passing through the stagnation point. They 
are computed by parabolic interpolation along the paths. 

The least square automation discussed in Section 4 is accomplished 
by repeatedly solving equations (3.4a-b) along the prescribed automation 
ares with varying initial data. For each automated parameter of function 
G given by (4.1) we assign initial data which has that particular parameter 
set to one and all the others set to zero. This set of solutions determines 
the coefficients of a linear system of equations for the automated parameters 
as described in [ 2 ]. The fumction G as defined by these autamated 


parameters is then used as initial data which determines the shape of the blade. 


The systems of two simultaneous linear equations (5.la-b) and (5.3a-b) 
become ill-conditioned whenever A, and A_ are nearly equal, i.e. near 
points on the complex sonic line; therefore, we expect our largest errors 
to occur in the vicinity of the sonic line. 

A third order accurate method can be obtained by doing the calculation 
with two different grid sizes and using Richardson extrapolation to the limit. 
We have incorporated this option into the computer program in order to obtain 
more accurate resolution in the vicinity of the leading and trailing edges. 
We cannot expect improvement near the sonic line because of the singularity 


of the hodograph transformation there. 
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6. Refinements to the Design Procedure arid Results 
A number of improvements have been made in the cascade design program 


that have not as yet been incorporated in the airfoil design procedure 
(cf. [ 2 ]). These modifications fall into three general overlapping 
categories. 

The first set of modifications was intended to decrease execution time. 
The cascade design program should take longer to execute than the airfoil 
design program since each path must connect the inlet and exit velocities 
in the initial plane. Since the gap to chord ratio varies inversely with 
the distance between these singularities we expect blades with low gap/chord 
to take longer to compute. Automation generally takes most of the computing 
time since the solution for X* and Y? along each automated path must be 
computed once for each automated parameter. Weighting for the least sauare 
automation was achieved by choosing additional mesh points in the corresponding 
interval for the airfoil design program. Since the execution time is 
proportional to the square of the number of points on the initial path, a 
highly weighted least square automation is very time consuming. In the 
present program the weighting is achieved by multiplying the value of w at 
the weighted point by the input weight factor. A significant reduction in 
execution time can be achieved if we compute X° and Y® for each automated 
parameter at the same point in the program. The right hand side of (3.4a) 
and ry would only be computed once. To achieve this we would require storage 
of the set of complex arrays for each X° and Y°?. For fine grids with many 
automated parameters the required storage would make this procedure prohibitive. 


A compromise, which allows for a fixed number of solutions of x? and Y? 
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at each pass has been incorporated. A typical run for the solution with 
ten automation parameters obtained at each pass has decreased execution 
time for the automation by a factor of about five using less than twice 
the storage. 

A second set of modifications was made in order to improve accuracy. 
The automation procedure was modified to ensure that it was second order 
accurate. The integration for least squares (cf. (4.2)) was modified so 
that only half the prescribed weight was used at the endpoint of each 
given interval. The automation near the leading edge was made second order 
accurate by adding an extra point on the path to the leading edge beyond the 
stagnation point so that central differences could be used for derivative 
approximations. 

An optional third order accurate scheme is provided which has been 
achieved by repeating the calculation at a finer grid and using Richardson 
extrapolation. Finally we note that the definition of A, (See equation 2.4) 
has been changed from reference [2] and so have the ensuing equations. 

This change avoids the inaccuracy which occurs when the denominator is small 
since for the flows considered u is greater than v throughout the supersonic 
zone. 

The final set of modifications was made to make the program easier to 
use. Some of these modifications have been achieved by more sophisticated 
computer coding such as the free-form input. Other modifications have been 
added so that more of the input parameters have a physical meaning. For 
example, instead of giving the values of x, and Vis at the stagnation point 


we now supply the leading edge radius of curvature and the angle between the 


central streamline and the inlet flow vector. Similarly the location of 
the trailing edge in the hodograph can be specified by prescribing the 
Mach number and the angle relative to the exit velocity vector. A number 
of changes have been incorporated in an effort to allow more general trans- 
formations for the hodograph plane as discussed in Section 2. An improved 
technique for the location of cuts of the logarithms in the initial plane 
removes the distinction between nose logs, tail logs, and side logs that 
was made in [ 2 ]. 

A computer program for cascade design has been written which incorporates 
these ideas. It has been used to generate a few designs. Each design 
required many runs in order to succeed but the time for each run is small 
so that computing time is not prohibitive. It takes about half a minute 
of central processer time on the CDC 6600 for a coarse mesh run. An inter- 
mediate grid and fine grid take about 2 and 8 minutes respectively. The third 
order option increases computing time by a factor of about 5. 

Figure 2 is a Calcomp plot of a blade which has been designed with 
this program in about 50 runs. All but the last few runs were made on a 
coarse grid. The inlet flow angles measured relative to the vertical stacking 
direction, is 45° and the turning angle is 29°. The gap chord ratio is 1.07. 
The Mach number distribution over the blade and the characteristics or 
Mach lines in the supersonic region have been plotted. With this Mach number 
distribution no separation is predicted for Reynold's number 1.6 million. 

The corresponding initial plane is illustrated in Figure 3. The input 
paths and sonic locus have been plotted. The plus signs represent the pre- 
images of points on the body and the arrows illustrate the location of 


logarithmic singularities in the initial function and the direction of their 


ay 


cuts. An asterisk is placed at the points corresponding to the inlet and 
exit velocities. We point out the similarity between this figure and those 
we have obtained previously with the airfoil design program (cf. [1, 2]). 

A number of experiments performed on airfoils designed by this procedure 
(cf.[5,6,7,8]) have proven the usefulness of this method. In the near 
future we expect a blade designed by this technique to be tested in a cascade 
wind tunnel at the Pratt and Whitney Division of United Aircraft. This test 
should verify the relevancy of this theory in this lower Reynolds number 
regime. 

Work to simplify the design procedure is in progress.A more detailed 
report on cascade design with an emphasis on turbine design should appear 


shortly [10]. 
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INPUT UtESCRIPTILOUN 


THE LNPUT FOR THE CASCADE DESIGN PRUGRAM CONSISTS Or TWO VATA FILES 
WHICd wt CALL TAPE7 ANO TAPCEG.)’ THESEs GCF COURKSEs WAY oe VISC FILES UR 
PUNCHED CARDS, VARIUUS FLOW PARAMETERS SUCH AS [HE INutl ANU EXLI CON} 


DITIONS ARe SUPPLIEV 10 THE PROGRAM 


DESCRIBING THE L2NITLAL FINCTION ANDO 
THe INTcCGRAILON 


BY MEANS OF IAPE7. 


AS LNPUT FROM TAPET7~ PAKAMcCTEcrS 
MAPPLNG FUNCTLON ARe ALSO ASdS1GNEV 


PATHS APPLAR ION TAPCE.6 


TO StT A VARIABLE LISTED IV THE GLOSSARY UF IAPE7 PARAMETERS we 
PUNCH A CARU WLIH THe VARIABLE NAME AND LIS VALUt OF THe FORMS 
PARAMETER = VALUE 


BLANKS ARE IGNORED 


ONE VALVE IS NEtVED 


COMMAS ARE YStyp TO SEPARATE THe VUMBERSe 


ARE DESLONATED BY ASSIGNING THemM 


PARAMETER CAN gpY SPECLFIED ON EACH CARD BY USING LITHteR A S&S UR § 


SEPARATUR, THE DATA FILE FOR TAPE7 


HAS ENUV= ON IT, 


THEIR OrFAULI VALUES. NOTE THAI 


CORRECT FORMAT 


TAPt& CUNTALNS THe PATHS OF INTEGRATION, 


THE STAGNATLON POINT tHe SECOND TO 


AUTOMATLON PATHS AND AS MANY INTEGRATION PATHS 


PATH WILL ALWAYS FE AUJUSTEU SO 


THE MACH NUMBER AT THt FRAILING 


SECOND PATH 


IN THE PLANE LIA=ETA(A). 


ANU MOVE CONVERSLON 15 AUTUMATIC. 


TU DESCRIBE A PARAMETER ¢ 


ANY 


ALL PARAMETERS WHICH ARE NOI 
THe OUTPUT LISTING OF 


fO Be JSe) AS INPUT TO REGENERATE 


THe 


THAI 
cO0Ge IS PRESCRKIBcy 
WILL BE ACCORLVINGLY AGUUSTLO. 


FOR SUSSONIC PATHS 


wdtN MORE THA 
Iece COMPLEX PARAMETERS s 
PARAMETERS TO ge AUIOMATED 

NON*=NUMcRIC VALULe MOKe THAN JUNE 
AS A 

MUST TERMLNATe WITH A CAKO WHICH 
PReSCRIscv WILL TAKE 
(APE? IS IN THe 


THE RUNe 


THe FIRS! PATH LeEAUS VU 


IRAILING cUuGeE FOLLOWED 3Y Ide NP 


AS D&ESIReu. THE -erRSit 
IT TERMINATES AT SIAGNATION, IF 
Or 


THE cNUPOLNI Tae 


ALL THe PATHS ARE PRESCRI38t9 


b4E PATH IN THE XI1-PLANe 
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OBTAINEV BY REFLECTIUN, FOR SUPERSONIC PATHS THe PATH LN THe X1=PLAVE 


IS O3TALNED BY CONJUGATING THE FULLOUWING PATH. 


EACH PATH Id DESCRIBED BY CAKDS OF IWU TYPES. [He FLRST TYPE oF 
CARD TELLS WHETHER THe PATH IS SUSSONIC UR SUPERSUNICe AOW MANY LLN 
SEGYENTS DeFINe The PATH, ANO THe NUMBER VF FURKS THE PATH wilt Avec. 
CARDS OF Trt SECOND TYPE CONTAIN THe ENUPOINT FOR CACH LLNtE sto%cNt LN 
THE ETA=PLANE AND A NUMBER wHICH IS USEL FOR LEAST SGUAXE WELGHIINS OR 


FOR LOCAL MESH REF INeMENT WITH NON-AUTOMATED PATHS, 


CARUS OF TYPE ONE CUNTAIN THE FULLOWLNG INFORMATIONS 


NCARO0S INTEGER, CULUMNS 123. THt NUMdbER OF LIVE StoYENTS WwHltd FURY 
THE PATH, THIS wiLl dt THE NUMBER UVF CAROS Or TYPt FwO wriltCH 
FULLUW. 

NTYPE INTEGER, CULUMNS 4-6. IF NIYPS<O A STEM FULLUNS AND <NETPc 


[S THE NUMBeR OF FORKS ISSULNG FROM IT. IF NIYPE=u A FURK OR 
SASE STEM FULLOWS, Jr NITYPESU A SUPERSONIC PATH FULLUAS; IF 
NTYPE=2% THIS PATH TS PLOTTEO?: IF NTYPE=e2 THE SINGYCARITIES OF 
THIS PATH ARE PLOTTED: IF ITYrP>2@ PAldA IS NOT PLOTTtLU, 


CARUS UF TYPE TWO CONTAIN THE FULLUWLNG LNFURYMATIONS 


POINT CUMPLEX, CULUMNS 1-20 AND 21-405. cD PUINT Ur A LINE déGM%e 
ON THE PATH UR FORK, 


WGHI Rte COLUMNS 41-50, FOR AUIOMAITIUN PAIHS Ag>(WoHI) LI * c 
WEIGHT GE THIS LINE SESSMENT JUStvu IV THE LEAST SQUARE PRO__SS, 
FUR NON=#AUTUMATED PATHS AdS(wOHl) 15 A FACTUR Voto J Su3= 


UiviIve THe FINITE OLFRFERENCE ORID UW THLS LiNce SceGMenvl, wodtl 
SHOULU -te A WHOLE NUMdcR, wWGnI<U “MARKS THE ficOLNNINGS OF Inte 
AUTOMATED StSMENT FOR AUTOMAILON PaldS ANDO dc BbeGANNING UP 
THE PART OF JHE PATH ALONG THE SONLE LIU rFUX A SUPCRSINIC 
PATH. 
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NAME 


RUN 


MRP 


NP 


NCR 


PEESLZ 


PPT 


JEFAULT 


1 


057 


GLUSSARKY OF IAPE 7 PARAMETERS 


ULFINITLON 


RUN NJIMBER 


MeSH REFINEMENT 
PARAMETER 


NU“3SER OF AUIO= 
MATION PATHS 


NU“WBerR Ur CUNE] 
STRAINTS 


PLUTP st Ze 


INYJEX FOR PLUI 
CUNTROIL 


CUMMENTS 


ABS(RUN) 2S Toe RUN NUMBER, NeGAe 
TIVE VALUES WILLE PRIUDUCE A wHITeE 
PAPER PLO! ON THE CYC 66uU0 AT tre 
AEC COMPUIING CENTEX AT NYU. 


DUUBLING THIS PARAMcTER CUIS Joe 
GRID LN HALF. A NeoATIVc VALUE LS 
UStU FOR IHIXI ORUCR ACCURACY, 


SINCe THede PATHS DeTeRYINe Trice 
RESULTING DESIGN THceY ARE PRINIED, 


Use NCR=5 FOR SINGLE CRITICAL PULNT 
AT THE TRAILING EDbe, 


USEU 10 SCALc THE PLOT OF teed 
PLANE. 


THLS IS A THree OLlotl NUMBER wHiCa 

COUNTRULS THE PLOTTING, fHe FiRSI 

VIGLU CUNTRULS THe eleTy INCARPEOT 

ANJ THe PLOT VF THe ReAL CHAARACTtro 

ISTIC PLAWE AS FOLLUWS? 

O NeITdeRK PLOUTTEV 

a CAARATIERLSTIC PLOT Owly 

2 LAKGE PLOT ONLY 

3-9 BOTH wIeke BE PLUITEU 

THE MLOOLeE VIGIT GUNTROLSL THE PEOT 

OF THe PAYSICAL PLAVE AS FULLUAD3 

0 NU PLOT GEVERATcCIU 

BY MACH NUYgc RX VISITXLBUTLON UNLY 

2 bLADE AND YACH DLSTRIBUTION 

3-9 SAME AS 2 PLUS CHARACTcRISTICS 

THE LEAST SIONIFICANT OIGLT IS USED 

TO CUNTRUL toe PLOT UF Tit AIT=9LAVE 

AS FOLLOWS? 

0 XL-PLANE NUT PLUITED 

1 SONIC LOCIS,. BODY POINTSs+s LIS» 
BRANCH POINTS, XLAe Xilbe XC 
POLE, NOISE ANDO TALL ARE PLOTTED, 

2 SA“c Ad 1 PLUS CLRCLES ORAWN AT 
AUIOMAIED POINTS. 

3 SsAMe— Ad 1 PLUS AYJTOMATtD PATHS 

4 SAME AS 1 PLUS NUNeAYTOMATED 
PATHS 

5 SAME AS BVJIH 2 ANU 4 

6 SAME AS BIUTH 3 ANU 4 wlIdnJsl 
SUPERSYNIC PATHS 

7 SAME AS BIUTH % AND Y 

B-9 SAME AS 7 PLUS CUNTUUR CURVES 

LF IPLT LS NeoOATIVce THEN tet rLUW 

COMPUTATIUN LS SKIiPPretOU. 
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NAME 


MAC 


RN 


TRANU 


TRANL 


KC 


ANGNOSE 


MTAIL 


ANGTAIL 


TR 


WLS 


GAMMA 


GRID 


MINLET 


ANGIN 


VEFAULT 


263 


210 


009 


20.00 


-15,G0 


0U3 


7260 


UtFILNITION 


MACH NUMBER 


REYVOLUS NUMBER 


TRANSITION ON 
UPPER SURFACE 


TRANSITION ON 
LUWER SURFACE 


KAULUS OF CUR-~ 
VATURE« 


ANoLe AT NOSE 


MACH NUMBER Al 
THe TRAILING 
LUSE 


ANoLeE AT THe 
TRAILING EDGE 


TRALLING EDvoe 
RAILO 


WreiGHTING FOR 
LeAST SQUARCS 
GAS CONSTANT 

GRLU SPACING 

PARAMETER 


MACH INLtT 


ANGLE AT INULI 


COMMENTS 


THE FLOW LS VORMALIZED 3Y ASSISNING 
THIS VALUt AT XI = Ue 


TO SKIP SBUUNJARY LATER IWIEORATLIUN 
SET RN = Ui 


THe VALUE OF X AT woAlCH we STARI 
THe FURBJLENT 3GUNVARY LATER 
INTEGRATIUN, 


PHe VALUt OF XK AT WALCH we START 
THE TURBULENT SUUNUDARY LAYER 
INTEORATIUN, 


THe RADIUS OF CURVATURE Al STAS- 
NATION S8ASED ION UNLI CHORD LENGTH, 


(He ANGLE BETWEEN IHe INecT FLOW 
VECTOR ANU THe CENTRAL STREAMLINE. 


If VALS LS OMELET VEU ole -VACUE 2S 
COMPUTED FROM THE VALUE OF XI wHiCH 
VEFINES .FHE cNOPOINE OF THE Se COND 
PATH PREStARI StU ON 'APCE. 


IHt ANGLE 3tIWwEeN Ite FRALLINO cUoEe 
ANDO THE cAIT FLOw veCror, 


THe RATLU OF THE 
TRAILING tJGe TO 


IHLCKNESS UF 
lot GAP 


Tre 


THIS CONTRKULS THE ReLATIVe welodile- 
ING speTweetV LEAST SaAUARES AND “MINI- 
MLZAILON, 


GAMMA IS 1.4 FOK ALX, 
1.07 FOR VRANIJM 


GAMMA 1> 


HAALFLOKIDc, 
GRLU/MRP 1S I[oae MAALMUM 


Mest Sidhe 


THIS VALJt ANY ANGLIN ARE USty Tu 
CUMPUTE X13, IF OMLTTE9S. LIS VALUE 


ANU ANGEXAT ARE Dt TcRMINEVD FROM XIA 
AND L1 WILL PRINT AS NeGatTive. 
THe ACTUAL INLET ANoLe WILL de Ute 


TERMINED FROY THE TURNING ANOLE, 
Ite IURNINVG ANOLE I5 THE DUIFRFcCRCNCE 
ANGIN|=ANGEXI I. 


NAME VEFAULT UcFLNITLON COMMENTS 


MEXIT MACA EXII THLS VALJe ANY ANGCALT Ake USeg 10 
COMPUTE Xin. IF OMLTTED, LIS Value 
ANU ANGIN Ake JETERYLNeD FROM ala 
ANO) 44% WEE PRINT AS NEGATIVES 


ANGEXIT -18.00 ANoLE AT EXIT UStU ALUNY wWLlIH ANGAN TO VelerMlNe 
THe TURNING ANGLE ANU EXIT ANoLE. 


XILA 0479203 XI AT A THE POLNT 1S Tot X1-PLANE CURRES= 
PONOING TU THe EXLI VELOCITY, iF 
MexL! LS S2ECIFIED IHe vALUt UF KIA 
WILL BE CUMPJTEU BY THE PROGRAM, 


X18 -.ls-.19 Xl Al 8B THe POINT 1S JHE XI-PLANE CORKES= 
PONDING TU THE INLeS VELOCITY, IF 
MINLtI IS :SPeCirFtieuv ihe: VAGUE DE 
Xig Wie ve COMPUTED. 


XIc -.12+.20 XL AT C THLS UeFINeS THe CHARACTERISTIC 
INITLAL PLANE FOR X5 AND Y3- 


POLE POLE LOCATION IF OMITT cu foe FRANSFORMATION FROM 
THe XIePLANE TO THE HUVOGRAPH PLANE 
WILL BE OETER“MINED SULELY 38Y Ide 
LOCATION UF IHE BRANCH PULNIS. 


BP 0054634 BRANCA POINT THE MAPPING FROM Toe XITePLANE Ty 
THE HUDOGKAPA PLANc LS DJeTeRYMiVNeEdD 
dY PRESCRIBING ITS 3SRANCH PULNTS 
ANU UVPTIONALLY A POLE. KS MANY A'S 
TeN SpRANCH POULNTS MAY BE SPEClFico, 


tT OevsG. CUcFFICIENT OF DETERMINES LIVEAR FeRM FOR THe INL 
EIA TIAL FUNCIION. THe REAL ANU/OR 
LMAGLNARY PARTS CAN ge AUIOMALZ Ue 
tETSQ Oer0e COLFRICILNT UF OLTERMINES QUAJRATLE TERM IN Tae 
t.TA«x«2 INTDILAL FUNCTILON, IHt REAL ANOsOR 


LMAGLNARY PARTS CAN 3st AUTOMATE. 


ETCu Ose. COEFFICIENT UF DETERMINES CUBSIC TERM FOR THE INle 
Li N«e*«s PLAL FUNCIIUN. THe ReAL ANUSUR 
IMAGINARY PARTS CAN 3te AUTOMATED, 


ETFO OorO, COLFFICTENT OF VETERMINES QJARTIC FERM IN Tee 
LI Ak#e4 INITLAL FUNCTLON, IHeo Qt AL ANDsJVUR 
IMAGINARY PARTS CAN Be AUTOMATEU. 


LOG LOULARITHM TERY tACH LOG LS JeFINevV SY PRESCRIZGING 
FOUR NUMBERS SEPARATLEU BY COMMAS. 
tHe FLRST TwoO SFECLPFY [THte ReAL ANU 
IMAGINARY PARTS OF IHE WELGHT ANY 
tlintR OR BOTA MAY ce AUTOMATED. 
THE LAST IWO SPECIFY [HE LOCATION 
OF THe LOv IN THE XL=PLANE, as 
MANY AS 12 LOGS MAY BE PRESCRIBED, 
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TRANSONIC CASCAJE VESIGON RUN 


RUN 
WP 
PLTS1L2Z 
MACH 
TRANU 
KC 
MTALL 
TR 
GAMMA 
MINLCT 

MeXII 
XIA = 0475 
SLC = ~-, 120 
BrP = 0,000 
el = 20%3 
ETGCU = 2106 
LOo = 2285 
SUG = c= 2976 
LOo = ~GU5 
LOlv = -=,010 
Los = eU55 
LO’ = 0,000 
LOS = 062 
LObv = -,085 


[APE 7 

-55 ; 
fs} ‘ 
1.000 ; 
2636 ; 
2050 ; 
-007 ; 
2986 ; 
2030 ; 
1.400 ; 
780 ; 
2480 ; 

» 2003 ; 
’ 206 ; 
een a, ; 
’ 0246 + ; 
’ e02U » H 
¢ 72966 ¢ ’ 
9 e724 *& ’ 
’ 0043 *& ’ 
’ eU6e * ’ 
t pa e0So ’ 
’ 0181 «* ' 
' 0171 »* ’ 
' -.064 ’ 
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PROGRAM CASCAUC(QOUIPUTe TAPE1Le TAPE2=OUTPUTs FAPEY=1LU00R, | APeES=1LU00 
1 3,TAPE7=1028+SF4PES) 

UIMENSION CARU(81) 

COMMON SUZ PI, AARP, NRN NX, NCR LOM, LPLT, “IDE NB gVCNIVZLL UI, KK 
LOMMON/SGZ NIisN2eNGeN7aN3eN4G."%1 

SECONU 15 a CUC 6600 CENTRAL PROCESSOR TIMING ROJIINE 

CALL SECUND(T1“E) 

Pe954.19159269357989 

CHANGE LBM WHEN CHANGING THE DIMENSION SIde 

LsmM = 200 

CHANGE KSBM wHEN CHANGING THE DIMENSION S12ZE OF KS AND YS 

K3M = 10 

M1 = NL 

REWIND NY 

‘IKPFAC = 1 

KREWLND Nid 

KeWlNU N6 

KEWIND N/ 

READ IN THE TAPEZT VATA 

READ (N7e120) (CAROQ(J) »« J = 1-80) 

Lr (LNCOVDE(CARU). VEU) GU TO 20 

MRP LS Nt Garivte wiatN CHOUSING {He IHIRU OXVER ACLURATE “Me TAdu 
MRPJULD = MRP 

“RP = JTABS(YRP) KMRPFAC 

tr \@MRPOLU SET.0) YRPFAG = MRPFAC+1 

Zr (NP.Geel) NY=N2 

VWP = ITABS(NP) 

CALE “ENE 

WRIJe (N2e60) NRN 

LEV CLP ET. 0)! 1GO" TO 1P0 

CALL AUTUL (N2eKBM) 

QV 40 J = 4,99 

“ODE = <J 

LALL MAIN(J) 

LE ACMODE «EQ ¢99))) GO TO ~4'5 

WRITe (N2eS0) J 
CUNTINUE 

ARTIE (N2@0130) JJ 
CAEL “SECOND “ClL) 

'iMe = T1-T1IMe 

WRITE (N2e110) FIME 
NT SSNS 

Lr (MRPFAC,EQ.2) GO TO 2 
1F (MRPOLOD.GT«0) V1 = Mi 
MRP = TABS(MRPOLD) 
CALE *TEILCE 

LPL T = IABS(EPLT) 

MX “= @ 

CALL BODYP7 (MX) 

CALL BLAVUE(™“X) 

IF (LPLT.EQ.0) CALL EXIT 


42 


CALL HOGRE (MX) 

CALL CPLUT((0610.) 9999) 

CAREC ESL! 

FORMAT (23H BEGIN CXECUTLIUN UF RUN ,14) 

FORMAT ( SH PATH 413, 9H FINISHtU) 

FORMAT (17H TOTAL CP TIME IsSeF6.1+8H StCONIS) 

FORMAL (80A1) 

FORMAT (23h OUI OF PATHS/18H LONGEST PATH HAS 414507H PUINTS) 
tNo 


FUNCILON INCUUE (JATA) 

UIMENSTON CATA(1) sVARNAY(7) yCHAR( S59) eVALUC (20) eVVAL(2) eVARLST (49) 
1 »sLUCV(40) eARRAY(1) eIARRAY(1) 

REAL MACHsLOGS 

LCUMMON SAS GAMMAyYMACHeROVUI( 4) 

CUMMUN /U/S PLleANeIRPeNRNe NX eNCReL OMe IPLT ese YUDE ge NGeNCeoNived bye Uig KX 

COMPLEX tTyU 

COMMON SES LUOS(4.12) eETIUIY) eXReNLUGS eo NK VP eNUe lL U1 352) 

COMPLEX XITAIL«+BPePOLE eXLAeXLBeXIC 

COMMON /H/ BP(1U) ePOLE se KALA ex be ALC eXLTAIL ePLISLewNePS 

REAL MACHA,YMACH3sMTIAIL 

CUMMUN SKS RNePCHJePCHL ss MACHA sANGLAMIALL se ANGLI se WACHRAYOLA RC, 

1 ANGLN,IR 

CUULVALENCE (ARRAY(1) IT ARRAY(1)) 

COUIVALENCE (tTsETIU) e(GPResP) »« (PePULE) 

VATA CHAR SLHAsLH3eLHCelLHI,LHZe LHP eo lLHGelLHAe LHI glide lHKeldbeltit, 

1 LHVeLHOsTHP eld de LHReIHSs LAT se LHUslLHVelHWelAxelLHVelet eye lHbGelalelH?s 
2 AHS elLH4GelLH5S el AS elLH7gLHBselLHGelH+esvlHeyIN,/ 

YATA VARLST /SHEND se 2HNP se SANRN « SHNCK,2HLD, SHYMRP,2ATR, ZHLUO ep eHEl, 
JYHE ISU seGHETCUsSY4HETFOs4HPULE g4HIPLI te 2HRNe SAWL Se GAHPLISI 2s 3HKL Ay SHK13 
2 9 SHXIC, 7HANGTALL «SHMTALL, QHRC, 7HANGNUSE » DHGAMMA,2HO2P,4nMACH , 

2 ZSHRUN,SHANGIN 7HANGEXIT e4HGRIi vr SAHIRANUsSHIRANL pOHMINGC CT eSAMCKLTZ 

VATA NCRU/O/ +» NLST /38/ « BLE /1H /eSTAR /1A*/S 9 K/1z7 2 LCXMP/S2/ 

LUCVINLSI4#1) = LOCK (LOCV(NLSI41)) 


NAS 
3) fa — A | 
MUDE = j 


1F (NCRUeNE.LO) GO 10 4 
LUCVI2) =-LOCF (NP) 


ie & 

LOCV(S) ==LOCF INR) 
Wau = 1 

LUCV(4) ==-LOCFINCR) 
VCR = 7 

LUCV(5) = LOCF(TR) 

(xX = .03 

LUCV ( 6) = =LUCrP (ARP) 
RP = ji 

LUCV( 7) LOLCF (IA) 


LUCV (8) = LOCF(LOGS) 
LUCV(9) = LOCFI(ET) 
LUCV(10) = LOCV( 3) + ICMP 


21) 


LOeVC11): = LOCV(10). + TEMP 
EOCV (12) =-EOCV (22). “+ ICMP 
UV 2JJS=1,4 

tlJid) = 0, 

VOCVlTSs) = -LOCrR(P) 

POLE = 0%. 

LOCVICTY) = -<LUCF(IPLT ) 
LPISE S="43:7 

LOUCV(15) = LOCF(RN) 

KiN= =" 705 

LUCV(16) = LOCF (XR) 

Ae =o 6 

LOCV (17) = LOCF(PLTSIZ) 
PERSPZ = 

KOE Va) t= 7 LOGV(1 3) "H+ ENCMP 
XLIA = (2«46,.-U02) 

LOEVCrTS) = LOCV(18) 4+ TCMP 
X1B = (tel,reld) 

EVEVC20 ) =~ LOUEV(1L9) -+- LEMP 
ALG = (-=e125220) 

LOcV(21) = LUCE (ANGET) 
ANGL! = -15. 

LUCV (22) = °KOCE (MTATL) 
“TAIL = Oe 

LUCV(23) = LOCF (RC) 

RC c=. 0079 

LOcV(24) = LOCF (ANGLN) 
ANGLN = ZO 

LOCV(25) = LOCF (GAMMA) 
CAMMA = 124 

LOCV(26) = LOCF (BPR) 

3BP(1) = (+6G5ee54) 

LOGY (27) -=—-LOCF (MACH) 

MACH = 265 


LUCV(28) =-LOCF (NRN) 


LOUCV(29) = 
ANGLB = 7. 
LUCV(30) = 


LOCF (ANGLB) 


LOCF (ANGLA) 


ANGLA = -16&,. 


LUCV(31) = 
AAS =. 08 

LUCV(352) = 
PCHJ = .05 
LUCVI4Z3) = 
PCHL = .10 
LOCV(34) = 
MAGHB = Ue 
LOCV(35) = 
“MACHA = O- 
Ae 
NA = 0 
NCRU = 
VV >L=1 
VARNAM(L) 

SAVE += 


LOCF (AA) 

LOCF (PCHU) 
LOCE (PCHE) 
LOCr (MACHB ) 


LOCF (MACHA) 


LOCF(LOGS) 


NCRO+1 


17 


= Un 
DATA(81) 


4y 


VAGA CSL t= LHS 

DOr SUG (be= 161 

LoNURE BLANKS 

LEE RCO ATA CE) SESH ") GO: SFO? 300 

ReSel ON cNCUUNTERING A VOLLAR SIGN 

IF (CDATA(L) ebU.1HS) OR (VATA(L) eEY1Hs)) SO TU 6U 
dr (NX.LS.2Q) GU TD 300 

Le SO ARAGL SEQ LAS 9 iGO! TNO woo 

FIND OUT WHICH CHAKACTER we HAVE 
BOMLORY = 71739 

tr (UVATA(L),Ede CHAR(J)) GO TO 2U 
CONTLNUc 

ga) AO 

tr (NX.EQ.G) GU TO 190 

LF (J.GIe26) Gv 10 270 

VARNAM(NAX) = JATA(L) 

NX = NX+41 

LF (NX.Ghee&) GY TO 260 

oJ TU 300 

CHECK FUR END UF DATA 

Dr (NX.LE.1) FOV TO 70 

PRINI 41UsiELeVARNAM 

CNCUVE (10.410eVNAME) VARVAY 

LF (CVARLS! (1) ¢bUeVNAME) GU FO 310 

SET VARNAM BACK 19 BLANKA 

UO 80 .J = Lev 

VARNAM(J) = 1H 

wo NX 

fv Geis 

eee QeNewG °C Aro esUo 

af Oc EO 0!) NK Sg 

LUCV(K) = ISLON(ITABS(LOCV(K) )+1eL0OCV(K4)) 
IF (KX,EU.1) GU 10 299 

LE MUBKXSGT F220)" GO? TO ero 

PIT NUMBER INIUY VVAL RIGHI JUSTIFIED 
ENCIVVE (20e41UeVVAL) (BeEeU = KXecU) oe (VALUE UH1) ed = DeKx) 
Wikre 1 

lt (MODOE.GT,.2) GU TO 110 

AN = LABS(LOCV(K))-LUCFI ARRAY) 

4F (MODE.-EG.2) GO TO 100 

INTEGER CONVERSILUN 

UeCIWWe (2Ue42UsVVAL) TARKAY( KA) 

tr (LOCV(IK).GIe0) ARKAY(AKA) = LARRAY (KA) 
vV [VU 3400 

CONVERT FTO REAL 

UeCUUE (204450eVVAL) ARRAY (KA) 

LF (LOCV(K),LIEe0) IARRAY(I4A) = AKRAY( KA) 4+.5 
MUHE = 1 

6Y FU 400 

AY = LOCV(K)=<KZ 

MODE! = i 

Ar (COCKY LEO) eURe (KY.GT.64%)) GU TU 120 
AUTOMATED PARAMETER PUT AYDRESS UN LISI 
VA = NKel 

LO(NK#1) = KY 


4s 


ae 


vV TU 300 
120 PRINT 470eVNAMtyVVAL 
vV TU 280 
150 NX = 0 
cCNCOVE (10041U0Ue*VNAME) VARNAM 
LUOK FOR THE VARTARLE ON THE LIST 
WO 160 K = LeNLSI 
LF (VARLSI(K) eb Qe VNAME) GQ TO 30G 
160 LCUNTINUE 
PRINT 4Y5U, VNAME 
6U TU 280 
190 LF( (ILE 025) eORe(JeEe40)) GO TO 200 
CHANGE THe 4YOUt UPON ENCOUNICRING A DECIMAL POINI 
Ir (DATA(L) tG@elHe) MODE = MOVEF1 
SAVe THe DIGLIS TO FURM A NUMBER 
VALVJECKX) = VALTAUL) 
Kix 4 KX 4+ 1 
bv TU 300 
CHECK FOR pf DECIMAL POINT 
CHECK FOR A CO4MA 
200 itr (UATA(L).-EWelHy) GO TY 99 


“JUDE = 3 
AX = KX+1 
00 TO 300 


270 PRINI 460 
280 PRINI G4Ue NCRUZ(DATA(U) oJ = 1480) 
LP = L+39 
PRIN' 410¢ (Bee J = 1lebLP)eSTAR 
SeT ERROR FLAG AND LOOK FUR NEXF DULLAR SION 
NAS a eed 
299 “ODE = 1 
400 CONTINUE 
310 VATA(@1) = SAVE 
INCOVE = Kel 
1F (LNCOUVE..NE U0) RETURN 
RESTORE LINE CUUNTER 
NCRU = 0 
VBPS= (LUCV(26) -LUCF(RPR))/2 
NLOGS = (LOCV(8)-LUCF (LOGS) )/4% 
ECG) = ECC) 
KETUKN 
410 FURMAT (8041) 
420 FURMAT (1206) 
440 FURMAT (t2e0.0) 
440 FORMAT (19HO****ERKROR IN CARD eL35e5H #€¥*# #2 IXe60A1) 
450 FURMAT (*0VARLASLE NAMED *,ABe * WAS NUT FIUUND*) 
460 FURMAT( *ONGCN©ALPHABETIC CHARACTERS NOT ALLIWED LN VARITAS3LC NAMES*) 
470 FURMAT ( *€OILLE GAL STATEMENT *A7ee@H= +2 Al1U) 
cNO 


46 


SUBROUTINE INIT 

VO LNITLALIZAITLON 

Re AL MACH 

COMMON /SA/S GAM“As4MACHsROVI»RUUTIL 

COMMON /C/ XARCAL yXAIMG e YAREAL se YALMGeAUI Moe YSIMGe TPTR 

CUMPLEX UeVeXleX2eX5eV1L eee SeLAMUAP eLAMLA%e XI acti 

COMMON XL1( 200) eWOHT (200) ec TAIL 20U) eU(20U) eV( 200) sA1l(20U),Y1(2U0), 
1 X2(260)*V¥2(200)¢X3(10e2UU) .¥Y35(10 9200) eLAMIAP(20U) eLAMUAM(Ze0U) 
COMMON JOS PLlreAhe ARPeNRNoNXeNCReLBOMe LPL Te MIDE Node VO a Mide hl eu e KX 
COMMUN SF/S XURIGsYURIGse SCALE ePSICPUOLI#VOLIeXOLVUIS50),10LU( 30)» 
PPP SLOL D'S) 

COMPLEX X1LTAI Le dPePOLE ae XLN ge Xl be KILL 

COMMON SHS BP( 10) ePOLE eXLAe XLS ee XLCeXLTAILePLTS1ZeNoPSsS 

REAL MACHA, MACHS,MTALL 

LOMMON SKS RNePCHIsePCHL se MACHA sANGLA,MIAIL se ANOLT s VACHS ,ANOL3ynKC, 
q ANGOLNe IR 

COMPLEX LeSy H eSOFXLALOFM,ROUT»RUUID 

ATA by £0. ele l/ 


SCALE «= Le 
XURLG = O. 
MOR LG: =. Ue 
eek Ge Sei 0.0 
RUOIL = (Oeele? 


RUOTE c= +61. Die) 

LAY DOWN PTH FROM XLA TO XI3 10 XIC IN EIA=PLANC 
LF (MACHA.GTe0O-e) XIA = XLIUFYMIMACHAsANGLA,XALA) 
dr (MACHBeCTeU0e) XIB = XLIUFYI(MACHB re ANGL Be ALS) 
MEO) a XT A 

COMPUTE THe NUMBER OF PULNTS FROM XI 10 ALZ 
MM = CABS(XLB-ALA)/SAAt+.,9 

AM = “AXU( 1, MMeMRP) 

Ne = MA+1 

MW = CALB=-XIJA)ZFLOAIT (M4) 

LAY VOWN GkKIU 19 x1 

YU 10 L = e,Na 

Atl) <= X12 wkd 

CUMPUTE THE NUM3tEX% OF PULNTS FROM XI8B TD KIC 
“4M = CABS(XICHEXALDB) /SAA+, 999999 

44 = MM&«&MRP 

VC = MM4+Nod 

rheo= (ALC@-XALBI/SFLUAT (MM) 

LAY VOWN GRIVU FROYW XIA TU XIC 


NSP = Nodtl1 
UO 2U Ll = NAPs NVC 
we OL) Ss XE 2 ee 


LAY VOWN REFLECTE) GrID IN THE KlePL AVE 
JU $U L = JeNC 

LiA(L) = CUNYO(XLI(L)) 

{AQ = Pl/1F Oe 

OAM = (LAMMAH1e)/2e 

LUSY = 1Lle/( MACH eMACH) +604" 

LF (MACHU eGTeVe) SU 10 4U 

FINU MACH NUMGeR AND ANSL«e AF INLeI 

> = SUFXI(AIB) 

YING?) = Le 


7 


C 
nN 


V(ING) = Ue 

LALL GETYUV (SeCONJG(S) eU( N53) eV (NBD) 

Jo = REAL(U(NG) ) 

VB = RKREAL(V(ND)) 

GSB= U3¥*UB4+VBFVS 

MACH = SORT (YSd/(CUSG"GAYMEUSD) ) 

ANGLB = ATAN2(V3,UB)/RAU 

bY TU 45 

ASK = COSQ+# YACH3*MACHB/ (16+ GAM EMACH GB EMAC) 
Js = SQRI(QSB)*«COS(RADKANOLSB) 

V8 = SORI (GSKB) €SILVIRADEANSLS3) 

U(NB)I = Us 

V(INd) = Ve 

LF (MACHAeGTe¢Ue) SO TO 5U 

FIND “MACH NUMbeR ANG ANGLEe AIT EXIT 

UY(1l) = le 
V(1) = Ue 
> = SOFXI(XIA) 

CALL FETUV (SeCUNIJG(S) eU(1) eV(1)) 

JA = REAL (UC(1)) 

VA RFEAL(V(1)) 

JSA JA*UA4VAKVA 

MACHA = =SGRT(USAZ(CUOSQ-bAMEQSA)) 

ANGLA = AILAN2(VA_gJA)/RAD 

y9U FU 55 

YSA = COSQ#YACAHAKYACHAS (1. + 6GAM*MACHAEYMACHA) 
UA = SQRT(CSA)*COS(RADKAWOLA) 

VA = SURI OSA) €SLNCRADAANGLA) 


UCL) = UA 
Viz) = VA 
RHO3 = RHO(GSp) 
RHOA = RHO( GSA) 


dr (MTAIL CTU.) XLTAIL = XIOFM(MTAILeNNGLA+ANGLI eX1A) 
CALL LAMBDA( UCL) oe V01) eLAMUAP(1) sLAMDAM(1)) 

ALASS = 1,7CABS(XIB]=-X1A) 

FACL = (Le-TR)*RHOAEVARRHUSkVB 

FAC2] = (leelTK) ®RHDAkUNERHIUB «UB 

TARLAL =\"FACA-FACL¥*REAL(LAMGAM(1)) 

JF (FACLTeNE Ue) YAREAL = YARCAL/S (FACIL *AIMAG(LAMDAM(1))) 


TL(1) = CYPLXCIVAREAL 1. )S5ORI (Let YAREAL®YAREAL) 
YI¢L) = (XTABS/CABS(LAYMUAM(1)))*(Y1(1)/CA43S(Y1(21))) 
A1(1) = =LAMUA%(1)#*Y1(1) 

LF (REAL (X1(1)).LTe-0-) GO TO 60 

Re) SP = X7 

WEG S =e) 

LCUNTINUE 

XAREAL = REAL(X1(1)) 

YARLCAL = REAL(Y1(1)) 

KAIMG = AIMAGO(X1(1)) 

YAI“MG = ALMAG(Y1(1)) 

ASIMG = =-XfIMG/(1.-TR) 

TSIMG = “$YAIMG/(lLe-Trn) 

IPTR = Seo¥*¥PI¥*IN 

KETJURN 

tNO 


48 


50 


40 
20 


100 


SUBRUUTINE AUTUL (N2eKB4) 

COMPUTE NOSE PATH, TAIL PAIH, AND K1 AJIOMAIEOD PAIAS 

Al AUIOMAIED PATHS FUR ZACK OF THe NK AUTIYMALEU PAKAYCIERS 
CUMPLEX Xl_l 

COMMUN X1(4C60) 

CUMMON SUZ PLleAAeARPeNRVe NX pNCReLOMe LPL Te 4IVDE ae NGe NC re Nie lLleuideKAK 
COMMON /E/S PARKAM(56) eXReNLOGSeNKeNPeNJOLC( 52) 

“WUGt=0 FUR VUNS AUTOMATED PATHS 

MODE =" '0 

fr (NK,EQ.C6) GI TI 50 

SET AUIOMAIED PARAMETERS 10 ZERU 

YO 10 Jeo = 1eNK 

Pee Et Ut 15) 

PRRAMCL) S Ce 


A2 = 2+MAXG (NPD) 

JO 20 MOVE = lekKe 

NX = MINO(KBMeNK4F1) 

No = 0 

CALL GTPATRK (NC+#1leNNol) 
Pees de 

lr (MOCE.GT,2) SO 10 30 
ee = NN@' 


1F (MOGE.EG.2) WRITE (N2,100) 
EXTEND THE PATH UE POI VI 

KL(NN41]) = XLONND F0OXLONV) @Xa ONNO1)) 
VN = NN41 

CONTINUE 

CALL MAIN (-1) 

REO NK. Lt eK SM)' 4GO TO 120 

VO 40 J = KOMeNKe KOM 

No As 

NX = MINO(KBMeNK41l“JU) 

CALL MAIN (-2) 

CUNTINUE 

CUNTINUc 

CALL AUTV2 (KB) 

SET THE “MOLE 10 .ZERO <SINGe AULOMALLON,.42S COMPLET. 
“UDE = O 

NU = CG 

aa | 

Xe TUKN 

FORMAL (25h) NUSE AUILOMATLUN CUMPLCIE ) 
tNo 


SUBRUUTINE AUIU2 (KB) 

SETS UP ANC SULVES THE EWUATIONS FUR [ME AJIOMAICU PAKAYMLIEXS 
LUMMUN SUS PLeANe ARP eVRVe NXP NCReLUMe IPL Te MODE eNGe UC oN Ved heute AK 
COMMUN /t/ PARAM(56) +e XReNLUGS eNKeNPyNIeLE( 52) 

LUMAMUN SPS KURLGS VURIGeSCALE sPSICeVOLUs VIL IeAUL UI SU), FULL S50) 
1 PSIULO(50) 

COMMUN/SG/S N1ieNVeeNSeN7eN35eNGLM41 

CUMPLEX ALITALL e SP ePOLE ce Ade ALB AIC 


4g 


Lo) 


10 


CUMMUN SHS SPCLU) ePOLE ea XLA CXL ye ALLeALTTAIL ea PLISIZeNyPrs 

REAL MACHA “ACHBsYTAIL 

CUMMON SKS RM ePCHUePCHL se MACHA sANGLA,MITATL es ANGLI 9 MACHR se ANGLSAyRLe 
1 ANGLNe IR 

LOMMON SKP( 2000) eARRAY (40041) eAKRKAY2(3U 080) eWGHIIKO) KALA BL) 
CUMPLEX X1L%7eXLPLaXINOSE eXT KX eUSZUWS UUW eUIL eUZ2UWeUULZS( 5U) eDIZU, OW 
ULMENSION PSINJUS(3C) ePSTIFAL (30) ePSLP1(90) ePSIM1(50) eXNOSe 136) eXPl1 
J (30) e¢XM1( 40) eXTAIL (30) sYVNOSE(30)¢YP1( 50) eY4M1( SUI es YTALL( SU) eKX( 30) 
WAT ALXU GAM O\e'e bh e:/ 

SAD = 186U.7P1 

UtLTA = AA/SFLUAT (MRP) 

LUMPUTE FHe NU“3ER UF RUWS ANU CULUMNS OF THE MAIRLX 

(HE RIGH! HANU SIVE UF THe SYSTcM VF ceUATLONS Ld STOKtevD IV 
COLUMN NSIZP 


VSIZ = NK+tNCRE?) 
NST ZPs= NSIZ+L 
NSPZRO = NSTZP +1 


SET IHE LEFT AND RIGHT HAND SIUt UVF JHE MAIRLX IY ZERU 

OY, 20-Jo = TeNsleZ 

UO 10 K = 1,NS1ZP9 

ARRAY (UeK) = O. 

ARRAY(NKt+2,NSIZP) = 1. 

ARRAY (NAF2],NSL12P0) = 1, 

LP RCNGR SEE. 3) "OO Pho 7S 

ANGL = 2e*(ANGLNFANGLB) /RAD+PI 

AU = COS(ANGL) 

AV = SIN(ANGL) 

ARRAY(NK4+7,NSIZP) = 1. 

MeAU 1N THE VATA FUR NOSE AUIOMAILUN 

KewiNo 2 

VAP = NK+1 

VO 20 JJ = LeNKP2%XBM 

VX = MINO(K8MeN<+2-JuU) 

NJ = JJ4tNXe] 

XeEAJ (N1¢260) 4K1e%2 

READ (N1<2G60) YUMLeVM1L eC XML (gy) sY¥M1(4U) ePSL4L OU) o> UWS BPN ex 
READ (N1*2GC) UsVe(XNOSE( I) e YNOSE(U) ePSINUS (JU) y UY = JUeNU) eXiNOSE 
READ! (NL: 260) *UPLeVP Le CxPLlCd) «eYPECd) «PSiPl (J) +s J = Lend) eATP1 
JO 2U Ko SA eNX 

Y= JJ+K-1 

ARRAYV(NK+2,.U) = 1, 

CUMPUTE XU ANU XV 

IZOW = CMPLXCXPLUK)-XM1L(K) ep YPLOKI MR YM1(K) )SCMPLX(UPLleUM1lyVMlevPl1) 
ARRAY(NK+6_,J) = XV¥REAL(UCOW) -XUeALMAGIOZIA) 

ARRAY (NK+7,JN) = XJ#REAL(UCOW) F+XV¥FALMAG(OZIA) 

VUZIVJ) HSCMPLX(XPL(K)4XM1°4) -2e*®XNUSE (JU) ec YPL(K)+Y%4L(K)-2.¥*YNIDE( J) ) 
CUNT INUE 

Jw = CMPLX(UPL=UM1LeVN1leVPL) 

UJW = CMPLX (UP1+JM1,VP1l+V%1) 

XEAD IN DOAIAN FUR TAIL AUYIUMALION 

VO 3U JI = 1eNKP X34 

NX = MINU(KAMsN<K42—-JJU) 

NJ = SJtNX=-1 

READ (N1e26G) KX1sK2 

READ (N1is200) UML eVALe (XML CODY VYMECU)« ePSLAL Gol ls JS pax) oe XDSL 


50 


50 


50 


60 


65 


5n 


IS 


70 


95 


96 


40 


ReAD (N1e¢206) JeVel(XTAIL( IJ) -VYIAIL (CU) ePSITAL(U)s Y = Ide 
READ. (NGe206.). URLeVPLe CXR EG). eoV¥PECGS) wPSLP?CI) « J= 
OU SO. IK a= «1:9 NX 
Y= JISeKel 
ARRAY(NK+5,J) 
ARRAY (NK+4,J5) 
ARRAY INK+t+54¥0) 
ARRAY (NK+8,J) 
ARRAY (NK+9,J) 
CONTINUE 

NPTS - 0 

A CNP SLE ese) GO TO 95 

ReAU IN LEAST SQUARE RESULTS 

UO (55 KP = deNP 

UU 5U JJ = 1sNAPeK3M 

WX = “ING (KAM se NK4+2"JU0) 

KEAD (N1e2CG0) 41e%2 

VU 50 K = 1,K2 

READ (N1+2C0) (XX(IU)e J = LeNX) eWGHT (KA) eX X(K) 
UO 5U JX = LleNX 

J = JSN4tdXHj 

ARRAY2(JU9K) = AK( JX) =PSINOS (J) 

CONTINUE 

NORMALIZE 2£NU WEIGHT ROWS 

VO 55 L = 1,K2 

SUM = 0. 

YJ 60 J = 1,NKP 

SUM = SUM + ARRAY2]( Jel) FARRAY2] (Jeb) 

WOHT(L) = KWGHI(L)/SUM 


PSILTAL( J) -PSINOS( J) 
(XP1L(K)-AM1L(4))/UELTIA 
(YP1(K)*¥M1(K))/UELIA 
(XPLIK)FAML (K-20 *XTAIL (JI IS(DELIA®UELT 
CYP1(K)4Y41(K)-2.*VYTAIL (UN) D/(OLCLIA*®VELI 


ion ou wou 


CUNILNUE 

JU 80 K = ieNKP 
UU BU J = KeNKP 
v0.60 L = 2,.h2 


ARRAY (OeK) = ARRAY (UeK) tWwoH{ (lL) FARRAV20 Se LD) eARRATAZIA,L) 
CUNTINUE 

VWPTS = NPIS4K2 

LONTINUc 

UU 70 K = 1,NAP 

FACT = 4e*xXR/FLOAT (MRP) 

vO 70 J = KeNAP 

AKRAT(JeK) = FACT eARRAY(JUeK) 
LCONTINJE 

AXITic (N2e250) NPTS 

V0 Jo J = 1sNSLZ2 

OU 26 K = JUANSIZS’ 


WARAY(JeK) = ARRAY(Ke J) 

UJ 40 J = 1,NK 

UNE = le 

ARRAY ( UtleJUt1) = AKRAY(Ut1leUt1l) + UNE 


LALL LEQLARKRRAT se ARRAY (1 eNOLZP) eNSILSYNeSe4UeGOe DET) 
WRITE (N20249) VET 

Ait = O- 

TIL = O« 

UUZ = (Uesl.? 

[HE SUFFIX O REPRESENTS IHE VALUCS FOR ABSIUZUW=06 


VJ) eXa FAIL 
LoeNX)eXfP1 


Q) 
N) 


XURGU = Ue. 
YURi0 = U, 
AILJ ma Oe 
TBE OY c= OFe 
PSrGo -==0). 
VOZO -= oC Oeg God 
Le ANCReST. 63), GO. TO. 7120 
AU = O% 
OVE Sie 
129 VV 100 J = 1eNKP 
PSIC = PSICHARKAY(JUygNSIZP)*PSINUS (J) 
XURLO = KOKIGHARRAY(UeNSL2P) *XNUSE (YU) 
TURIG = YURTGOHARKAY (Je NSLCZP)*¥YNOSE(U) 
XIL = XTL +¢ARRAY(SeNSIZ7) *XTALL(IS) 
YIL = YIL +ARRAV(SeNSIZPI*YTAIL (J) 
JUZ = DUZFEARRAY (Se NSLZP) *#00ZN(U) 
JF (NCR.GT.3) SU TU 105 
XJ = KRUMARRAY (IU eNSIZP) €ARRAYOINK+6 00) 
AV = XVOEARRPY (JeNSIZP) KAKRAY(NK+7 eJ) 
40 TO 100 
105 PSICU = PSICOFARRAY(UeNSIZPO) *PSTNOS( JU) 
KORGU = KORGUFARRAY (GeNSLCPO) *XNOSE (JU) 
= YORGUFARRAY (Ge NSLCPO) *¥YNUSE (J) 
= XTLO + ARRAY (UeNSLZPO)*XTALL (J) 
YILO = YILU + ARRAY (JUeNSLEZPO)*YIALIL( J) 
= DUZEC+ARRAY (UeNSIZPU) &€VD0ZN (JU) 
100 CONTINUE 
COMPUTE SLUPE VF CENTRAL STREAMLINE AT STAGNATIUN 
ANGLE = eD¥*¥(ATAN2(-XVe-XU) tPi) 
ANGLN = RAG + ANGL-ANGLB-39U. 
UNS 2 E=LosCOs (ANGLE) 
VN = 1,€-4645IN(ANGL) 
CNX = COS(3,.*ANGL) 
ONX = SIN(Z,¥*ANSL) 
UZDw = CMPLX(XIUeXV) 
USZOWS = 4,4 (U0Z2-IZ0W*00W) /(UW*DW) 
do = REAL(DSZUWS) *SNX-AIMAG(USZOUWS) &CNX 
CHORD = XKXTL-XVURIG 
1F (NCR.GT.3) GO TO 115 
XC = 3.*(XUKKUFKVEXV) Z(30*CHORD) 
ASDZUW = SGRICXUeKXUtEXVEXV) 
90 TO 125 
FINI MAGNITUTE OF U2ZUW TU SATISFY RIGH!I wIUSE CURVATURL 
115 CUO = XTLU-xXORGYU 
VELL = CHORD=-CU 
USZIWS = 4,40UUZ0/(UWeUW) 
BO = REAL(DSZUN5) «*SNX-AIMAG(USZOWS) €CNX 
UcLs = BuHYvo 
YVISCR =( RC (DEL C#K3ZU-“UELB*CO) )¥* 2412, * KL KC 0430 
LF (VISCReGEeUe) GO TO 12d 
KC = =126*HO*CU/(JELC*¥BIU-JELB*CU) €*2 
WRITE (N29190) RC 
UISCR =- U0 
125 ABDZUW = (RC*K(DELC#BO+GEL dD *#CO) +SQRI (DISTR) IS (60e72e¥*RO*KUECLB4IELC) 
KRORILG = KOKGO+H3UZDWe(XIKILG-XOROU) 


52 


110 


190 
200 
240 
250 


10 


40 


TORIG = YORGU4A3UZUW4 (YORKLG-YORGU) 

YIL = YTLO¢ABUZOWe(YTL-VILO) 

PSIC = PSICUFAB0ZOW*e(PSIC=PSICO) 

CHORVU = CO¢AKBUZUW*EDELC 

SCALE = 1-/CHURV 

PSTLe =" PSIC/CHORD 

TIL = (YIL-“YORKLG) CHORD 

STORt THe VYALULS OF THE AUTOMATED PARAMEI CRS 

VO 110 J = 1¢eNA 

CC I+1 ) 

PARAM(L) = ABUZOWEARKAY (Jt1leNSLIZP)+(1.7ABUZUW) FAXRAY( Ut1eVSIZPO) 
XEWwLND Ni 

WRITE (N1e260) (PARAM(L)« L = 1¢56) 

AXITE (N10200) Gov eYTL 

WKITte (N1e20U) UNeVNeXITVUSE 

Ke TURN 

FORMAT (SOHO**IHAILS VALUE JF RC CANNOT SE ACHIEVLI#*® YolNG <C=r6,3) 
FURMAT (4020) 

FORMAT (14+ DETERMINANT = 4612.5) 

FORMAT (10H THeke ARe 614e30H POINTS USEY FOR LZEAST sQUARCS ) 
tNo 


SUBRUUTINE MALINC(ITYP) 

CUMPLEX UsVeXlLeX2eXSeV1leVer VY Se LAMUAP sLAMDA% XLTeclA 

LCOUMMON X1L( 200) eWGHT (200) ec TAC 20U) eU(20U) pVI 200) eAL( enue V1 eu), 
1 X2(200)+Y¥2(20U) 4X53(10,2U0) ,¥3(1U2200) eLAMIAP(20U) eLAYMDAN(20U) 
COMMUN SUZ PI eANse ARP eNRNeNXeNCReL OMe LPLT se MIJEeNSeNCeNve lb etter KK 
CUMYMUN/SGS NleN2eNGeNJeN3e V4 eM1 

COMMUN /SL/ XARLAL Ae XALMGe YAREAL eS YALMG 

CUMMUN FS XURLGSYURIGe SCALE sPSICsUODLD+eVULIeXOLUI30),YULUI50)» 
1 PSIOLO( 50) 

CUMPLEX UNNe VNNeXITNNG VINNG TEMP 

LOGICAL [Sw 

JATA ISW /,FALSE/ 

tr (1SW) GO TO 10 

ESW4 =F) TRUE Ss 

UA = REAL (U(1)) 

VA = REAL(V(1)) 

CALL SAVE(1,1+1) 

Xt AU IN THE NEXT PATH 

Ar (LIYP.G6E.0) CALL GTPAITH(NC4#1eNNeol) 

IF (MODE-EG.99) RETURN 

bET UsevVeXlLeX2eL AMJA4eLAMUAne, ANU CIA UN XL = CUNIJO(XTA) 

Jizz) = UA 

V(1) = VA 

Ret) = XAREAL 

T1(1) = TAREAL 

CALL GETAYG teeNN) 

Ar (ATY¥P.6t.0) Ll = JJ 

JO 4U LX = JeNX 

PSITOLO(LA) = PSIC 

CHECK FOR SUPLRSOVIC PATH 


60 


30 


50 


20 


70 


lr (KK.61.2C) FY TU 50 

SUBSUNIC PATH, COMPUTE SULUTLON IN THE TRL ONGLE 
lr (KK.EQ.0) GU TO 60 

SAVE UeVeXte AND Y1 AT END OF PATH SO JHALT WE MAY AOD FORKS 
UNN = U(NN) 

VNN = VONN) 

AINN = X1(NN) 

YTINN = YICNN) 

CALL TRIANG(2eNN) 

CHECK TO SEE LF THERE ARE ANY FORKS 

Her o(KK ees’ Os the TORIN 

NFORKS = -KK 

SAVE REAL UeVeXeY ANU PSL AT ETACNN) 


URNNV = UULE 
VRNN = VULG 
ARNNW = XULD(1) 
YRNN = YULOG1) 


PSINN = PSJOLU(1) 

CUMPUTE SOLUITIUN ALONG EACH FORK 

VU 30 K = 1eNFURKS 

VNP = NN4+1 

RtAJD IN THE FORA 

CALL GTPAIH (NNP» VF el) 

dr (MODE-EQ.99) RETURN 

OeT UeVeXleX2eL AMIA+e¢LAYMUA=, ANU CIA ON XL = CONULIXTA) 
U(NN) = UNN 

VIN) = VAN 

X1(\N) X1NN 

YT1(NN) = YINN 

JULD = URNN 

VULJ = VINN 

XOLU(1) = XRNN 

YULJU(1) = YRNN 

PSTOLD(1) = PSINN 

CALL GEIXYG (NNPe NF) 

CALL RECIAN (2eNNeNNPoeNF) 

“ = NF+¢1-NKN 

WRITE (N1*SC) KKeM 

WRITE (N1*SG6) JGULIsVULDeXOLDO(1) eVYOLU( 1) ePSLOLD( 1) XL (iN) 
CALL TRIANG (NN41—-NF) 

CUNT INUE 

RETURN 

SUPERSONIC PAIN, REAU IN THE OTHER PATH 

LF (LTYP.GE.71) CALL GIPATH(NC4+1¢MM,1) 

le (MOOL-EG,I9) KETUKN 

CALL RECTAN (2eeJIdeeeNN) 

COMPUTE FHE COMPLEX STREAM FUNCTION ON XI (JJ) 
DO 20" = Nee lt 

CALL GETPSI (LeJJ) 

CONTINUE 

POLI =" PSICGEUCL) 

cS ac es 

NX = NN 

“ = NNX+1l-I1 

LF ((N12EQ.3) eAND-(MOUD(Me2),.EQ.0))60 TU 75 
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WRITE (N1050) KAeM 

HEMP = “EMPEX (REAL X EGU SREARCEPAGCLE IFES LON CL OO oe AI MAG GALES) 2) 
WRITE (N1e9C) JULDeVULDeXILD(1) pYULO(1) ePSLOLU(1) 9 ter 
1F (JeEQeMM) GO TO BO 

CALL SAVE(CII+1eNNXsJ) 

WNX = NNX-j 

oar ee oan b 

CALL RECTAN (Jede2eNNX) 

HSPOLOELY c= PSLL 

CALL GETIPS1 (IieJ) 

PSTiies: PSIGLU(1) 

SOMEO 7.6 

CONTINUE 

FURMAIT (%020) 

RETURN 

aT) 


SUBRUUTINE TRLANG (KleK2) 

COMPUTE FHE SOLUTION LN IHE TRIANGLE xl XI (Kl) Al clAs 
cIiA(K1) TIA cTA(K2) 

COMMUN /SA/ SKP(2)eKOUT 

LCUMPLEX Use Xl eXSeXZeVlsVeeVYSeLAMUAPeLAMDAMsATsclA 

COMMON XL1(20U) ewG4T (200) ec TAC 2UU) eU(20U) ,V (200) sA1l(2Gu) +L 2uuU), 


1 X2(200) © ¥e(20U) 6X3(10e2UU) eVSC1U e200) eLAMIAP(20U) sLAMUKM(eOU) 


CUMYON S/O PLleAWeARPeNRNeNX se NCRe LBM LPL TeMIDE eNSe NC oNivelb edd eAXKX 
COMPLEX LEMP °¢ ROD 

VU 30 J = K1eK2 

SET THE GRANCH FOR THE SUUARE ROOI IN LAMJA+t ANU LAMDA~ 
ROOT «= (Ocele) 

GET VALUES GF SeVeX1L ae X2eLAMDPAteLAMUA] AT AL (U)ectIA(J-1) oY 
KEFLECTION OF THESE QUANTLT1¢€S AT XI(J"1l) et TA) 

Ulyel) = CONJol(JU(J)) 

Videl) = CONJG(V(J)) 

Al(J-1) = CONJo(xX1(J)) 

Yl(J71) = CONJO(Y1(J)) 

teMP = CUNJGILAMVAP (J) ) 

LAMUAP(J-1) = CUNJG(LAMVA(J)) 

LAMJAM( Jel) = TEMP 

CHECK TU SEE IF WE ARE PAST XKI(ND) 

LF (JeLEeNo) GY TO 1U 

GET X2 ANDO Y2 SY REFLECTIUN 

A2(J-1) = CONJo(X2tJ)) 

f2(jJ-1) = CONJG(Y2tJ)) 

CHECK TO ‘SEE LF WE ARE “PAST XDUNC) 

LF (JeLEeNC) GU 19 20 

Get AS AND ¥5 SY VWEPLECTFIUON 


UU 25 EX = 1LeNX 

KRS(LAed=-1l) = CUNUGIXSILXeU)) 
TS5(LXeda1) = CUNUGIYS(LXed)) 
LUNTINUE 


beT SOLUILON AI XI(U) FRU ETA(U) 10 clatAe2) 
CALL LINtJ(UedeK2) 


40 


30 


60 
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20 


30 


Ar (SeGEe NC) CALL SAVE (JUeJded) 

CHECK TO S€éeE LF We ARE UN A PATH UK FORK 
DE (OL NE) "GOP TOs Se 

WE AKE ON A PATA. SAVE VALUES AF X1( U0) SETAIKE)® SU THATS HES ean 
LATER EXIEND IAIS PATH wllTH FORKS 

Nga) os Wa 2) 

WiGgiad). ‘sc Vike) 

A= a=) xa CKe) 

Ya 1) > = YI we) 

LAMUAP(JW1l) = LAMIAP (Ke) 

LAMIJAM( Jel) = LAMDAM( Ke) 

1F (JeGTeN5B) GU TO 4u 

LF(JeLTeNB) GU TU 36 

CALL GETX2(JU+1leK2) 

RK2(J-1) = x2(K2) 

Ke(u=1) =: yetKe) 

LeCsgok Tt. NE) GO TO: 30 


UO 60 LX = JeNX 

XSCEKe Jel) = KS(LXeK?S) 
VSCEXeJ-1) = VY3QLXeKe) 
CUNTINUt 

CONTLNUE 

tCND 


SUBRYVUTINE RECTAN (K1leKe043eK4) 

SULvE IN THE RECTANGLE xl(K1) xl xXI(K2) ¢« ETA(K3S) ETA EDACKE) 
CUMYUN SAS SKP(2),RO0UT 

COMPLEX UeveXleX2eXSeV¥leV2eVSeLAMDAP yLAMUAY  xXIetia 

COMMUN XL° 200) ewGHT (200) ec TAC OOU) eU( 200) VI 200) eA1( 200) ¢Y1(ouUU), 
1 X2(200) © ¥e (200) 5xX5(10520U) .Y3(1U1200) eLANIAP(20U) sLAMDAM(2e0L) 
COMMON SUZ PLleAheAKPeNRNe NX eNCReLBMe IPLT se VIDE eNB ONC eNivedilLedJdeAK 
CUMPLEX ROUUT2 » ROOT 

YO TO J = KieK2 

GET JHE VALUE OF THE SOLUTION A AIL(J)cETAUK3S) 

PREY CKRKSCESG) GU TO 20 

N SUPERSONIC PATH GET UeVexyeX2eLAMDA+eLAMIA= ON EIA = AIA 
CALE GETXY2? «J) 

3U TU 30 

UCKS=1) = U(d=-1) 

ViGK S21) <= v(d=2)) 

X1(KS-1) = xX1l(Je-1) 

WL(KS=2) = Vie) 

LAMJAP(K5-1) = LAYMODAP (Jel) 

LAMDAM(K3-21) = LAYMDAM(J-1) 

X2(“5=1) = x2(Je-1) 

Y¥2(%3-1) = Y2(Jel) 

UO 25 (LX Sy Nx 

RS(LXeKS91) = AS(LXeIH1) 

YS5(LAeK3—-1) Y5(LXeIU-1) 

CUNT 2NUE 

KEEP TRACK OF THE BRANCH OF THE SQUARE ROUT FOR LAMUA4sLAYIA= 
RVOTe2 = ROOT 
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40 
10 


CALL LINtCU (UGeK4e KY) 

LF (JeEQeNB) CALL GETX2(MAXG(NB41043) eK4) 

KUOI = ROOT? 

SUPEKSONIC PATH, COMPUTE PSI ON EIA = LCONJUG(XIA) 
LE ((CAK LEO) URS (UeLTeNGC) eOReltebE.JJ)) GU 10 1U 
IF (J.GT.NC) GU 10 4G 

CALL GETPSI] (=NCeJ) 

6U TU 10 

CALL GETPSI(NCeJ) 

CONTINJE 

RETURN 

CNbD 


SUBRUUTINE LINE J(J1sK1eK2) 

SOLVE THe EGUAIIONS ALONG THE CHARACFERISILC X1L1(4) FRUYM CIA(AlL) 
10 EIA(K2) 

COMMON SLES XAREAL se XAIMG se YARENL se YALMGe XbIMoe YEIMGHIPIR 

COMPLEX UsVaeXlL er KX2eX3eVLle Vere VYSeLAMUAPeLAMDAYM AT ectla 

COMMUN XKL( 260) ewGHT (200) ec TAI200) eU(20U) ev(200) eA1(20U).71(2UU), 
X2(200) es ¥2( 200) oXS(10e2UU) eV 3(1U2200) sLAMIAP(2UU) sLAMUAM(eOU) 

COMMUN JUS PleAN se IRF aNRNeNXeNCReLOMeLPLT, “IDE Node NC oN Ve Ll ye Jte KX 

COMPLEX UX] setEF ecb e RH] RAS eRHS te RHY se L AMS eLANAIMeLAIBPH, LAMB A4s USe V3 
© LerATHeETAHeXSNT( SU) eIJLAMT * OY eE IGE 

DAT ACL: £10.01 67 

Y= Jil 

LE UCU NE CNC) cORS(K26GTNC)) GO. TO 5 

WU 2 LX= 1_NX 

CALL INITFN(ETAINC) eX 3NC(LX) LX) 

YS(LAeNC) = (06,02) 

XRS5(LAeNC) = ASNC(LX) 

4 (KK.GTeG) GU JO 2 

AS(LAeNC) = REALUXSNC(LX)) 

CONTANUE 

UXT = xXI1(JU)-AL(J-1) 

KALH = .D*(XI(U) 4x1 (J-1)) 

U0) BOnkL = Ki «Ke 

SOLVE FOR U ANU V TO FIRSI ORVER ACCURACY 

MHz = V(L)-LAMUAP(L) *€U(L) 

RH2 = V(bL=-1])*LAMDAM(L=-1)*J(L=-1) 

JS = (RH1L=-RH2)/7(LAMDAM(L=1)-LAMUAP(L)) 

VS = RH1+LAMUAP(L) *US5 

LCUMPUTE LAMUA+e+LAMUA]e TO FIRST URUER ALLCURACY AI XL (GdectAtL) 

CALL LAMBUA(US+VSeL AMBP +L A444) 

LCUMPUTE “LUPOLNI VALUES Ur LAMBGUA+ ANU LAYSIA- 

LAM3SPH = ,Se¢(LAYMBP+LAMOAPIL)) 

LAMSMH = eS#(LAMGYF#LAMDAM(Lel)) 

SULVt FUR U ANU V TO SECUND URUCR ACCURACY 

MH] = VILI-LAMSPHeU(L) 

KH2 = Vll-jJ)“LAMBYreU(La=1) 

JIL) = (HY -KA2)/( LAMHMH=-LAMSPH) 

V(L) = RHL4+LAMSPHeU(L) 

FINY LAMSGOA+ ANU LAMBUAe AT KIC dd et TAL) 


- 


20 


16 


50 


40 


60 


30 


CALL LAMBDA (UCL) eV(L) LAY BP eLAMYWM) 

CUMPUTE MIDPOINT VALUES UF LAMUA+ AN) LAYMSIA}- 
LAMSPH = e5#(LA4MBP+LAMUDAP(L=1)) 

LAMSMH = oS (LAYMBY4+LAMOAM(L)) 

S1Oxkt THe NEW VALJES OF LAY30A4 ANU LAMBDA= 
LAMVAP(L) = LA“2BP 

LAMUAM(L) = LANBM 

COMPUTE X1L-AND Y¥2 AT. XT(J), ETAL) 

RH] = X1(L)4¢LA%34%H*YL(L) 

RKRH2 = X1l(L-1)+LAM%3PH*Y1(L-1) 

T1(L) = (RH1"RH2)/(LAMBYMH-LAMBPH) 

X1(L) = RH} -LAY3AMY*Y1(L) 

CHECK -TO SEG IF We HAVE KEACHED Xls 

LF ((UeLE NS) CORD (LeLTeNS)?) GO TO YU 

CUMPUTE X2 ANU Ye AT XI(J) ETAL) 

RH3 = X2(L)¢lLAMSMHA*YVS (L) 

LF (L.GTeN&B) GO TD 20 

CUMPUTE X2.'¥2 ON EWA = 01S Xl = xXECY) 


KHy = I*(,314G+LAMDAP(L)*YBIMG) 
Y2(L) = (RHZ"-RA4) /(LAMBMHELAMDAPI(L) ) 
K2(L) = RHSZ-LAMSZMH*Y2(L) 

GUO TO 99 


RAY = X2(L=-1)+LAMSPH*Y2(L-1) 

Y2(L) = CRHZ"RHYG)/(LAMBYH=-LAMBPH) 
A2(L) = RHZe-LA“BMA*Y?S(L) 

LF ( (SLT oNC) eOR(LeLIeNC)) GO TO 9D 
CUMPUTce THE LNHOMOGENEOUS TERM 


Cr = (RHI +L RE CKALMG+LAMBYMARYAIMG) )AZOUXT (1) HALA) + CRASS tLe (X81 454+L AM 3M4 


1 *Y31MG))I/( x1 (N83) =-XIH) 

iF AE eG .NC) GO TD Go 

GE c= (0.00) 

CHECK FOR SUBSOVIC PATH 

LPF -OKK.GI.0) GU TO 16 

PATH IS SUBSONLC, GET SYMMETRIC EF AND GE 

SR 6= 5 EF 

EPA = sSECETACL) FETAL =1)) 

Se = ,5¥*((KH2-Le(XAIMG+tLAY3BPHEYAIMG) DVV ETAL) -EIAH) + (RHy-T* 
1 (Xx3L4G+LAM38PH¥*YBIMG) )/Z( CI AC NB) -E TAH) ) 

tr (J.NEeNC) GUY TI 30 

SET iXde VS GN LNT REAL CHARACTER TS TLC Xb = XLC 
tlGe = (cTA(L)I-ETA(L=1))¥*GE 


UOU 50) xX’ = 1x 

CALL INITFN(ETA(L) eXS5(LXeL) LX) 

dF (KK LEQ) XS(LXeL) = eBS¥(XSCLXoL) tHCUNUOCKX SINE CEX))) 
YS(LXeL) = YSCLXeLH1) + (LZTGER(XS (LX eL I AHeXS(ILXgLOL)) ) LLAYSPA 
CUNT INUE 

60 FU 90 

TF ( (AK LEO) eURe(JeEULNC)) GO TO 90 

UY = (OXL*EF)/LAM3BMH 

UO 60 LX = JeNX 

TSCLA¢ EC) HS -VYSCLXeE) +. OY 

CONTINUE 

vVU 1U 9D 


LOMPUTE X3 ANU Y3 AT X1(J) FETAL) 
UY = UXI*EF 
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cCIGE = (tTA(L)=-cTA(Lel) ) &OE 

ULAYL = 1./(LAMSMH=-LAMBPH) 

DOVE OD MEX = T ONX 

RAL = XS5(LXeL) FLAYMBMHFYS(LXeL) + UY 
RHQ = X35(LXyeL=1) +L AMBPHEYS(LXeLaH1) + CIGE 
YS5(LAeL) = (RH1L=RdA2?)40L 441 

XS5(LAeL) = RHL@-LAYMBMH* YSZILXeL) 
CONTINUE 

CONTINUE 

CONTINUE 

Me TURN 

-NO 


SUBRYVUUTINE SAV (41¢K2eK5) 

SAve THe SOLUTLON IN THE REAL HODOGRAPH PLANE FRUM ETA(KAL) Tu 
tTA(K2) ALUNG KI(K3) ON TAPE 

COMMUN /L/ YAREAL se XAIMG s YAREALs YALMG se XOIMGe YKIMGe IPTR 

COMPLEX UeWVeXlLeXSeXSeVle Vers Y¥SeLAMVAPeLAMUAW XT ec la 

COMMON x1(200) ewGHT(200) se TA( 200) eU(20U) eV( 200) sALl(2NU)d se VY1L(2uU), 
1 X2(200) eTe(2U0) oX5(10420U) ,Y3(1U9200) eLAYMIAP(20U) eLAMUAM( eOU) 
COMMON SUZ PL eAAeAMRPeNRNeNXeNCReLBEMe IPL TeMODE NB e NC eo NNe ll yeild gKK 
CUMMON ZFS KURL Se YORIG se SCALE sPSILeUOLIeVCLUyXOLUI(30),YOLU( SO)» 
1 PSIUL0D( 30) 

CUMMON/SG/S NIW*N2eNGeN7eNSe V4 QM1 

CUMPLEX ALTALL eBPePOLE se XdHeX1BeKIC 

LUMMUN /H/ BP(10) ePOLE se XLAVALBA ALC ae XL TAIL eae PLISL2eNBPS 

REAL LUG1sLO0G2 

COMPLEX [FeMP 

LUGILCAL 1Ldw 

YATN ISW/eFALSE./ « TOL /1,.E°6/ 

FE GLSW) 460 8049 

ES We=". PRUE. 

Us = REAL(UI(NS)) 

VS = REAL(V (ND) ) 

YPSI = RHO(UB*eUB4VB*V8) €(UBeYRBLAGHVeeKNIMY) 

ReTUKN 

LFe(CCABS UXT CKSeXECALE) ek be LOL) GOchY 

IEMP = CLUGCCUXL(K3) -XITALL)IsS(XIT (KS) -XIA)) 

ALPH = ALMAG(ITCLMP) 

CUNT LNUE 

femMP = CLUG CUKLIAS) eXIBIS(XLIKS) MALA)) 

LUG1 = REAL (TEMP) 

(HEITA = ALMAGC ICMP) 

LF (K1.EUeNC) IHETAR = Oo 

1eMP = CLUGUIALIKS) =XITA) ©&CMPLXICUSITHLIA2) eeSINGINeTA2))) 

LUG2 = REAL( TEP) 

IHETAD = AIMAG(TEMP)4+THEIAe2 

VU 1U L = AleAKe 

UNEwW REAL (U(L)) 

VNEwW REAL(V(L)) 

AXT = SCALCECXCALIX2(L) PFLOGLFE( REAL (KL(L) P#REALUCA2(L) DD FLUE! 
1 =“XSLMGeIHETA = (XATMGeXdIMGI#IHETAQ-KURT So) 


98 


15 


50 


30 


40 


TYT = SCALEKIREALCY?S(L) ) *LUGIF+(REALCY1L(L) DtREALIVY2(0L))) *LUS2 
“Y3LMGaATHETA = (YAIMG+YBIMG)*THETAC-YURIG) 

YU 50 LX = LeNA 

XNEW = XXT + SCALE*REAL(AS(LXeL)) 

YNEW = YYT+SCALE*REAL(Y3(LXeL)) 

Le (eb seo. NC) GO. T9° 20 

UH = e-5*(UNEWFUOL,) 


VH = .5*(VNEWFVOL)) 
UX = XNEWeAULU(LX) 
UY = YNEW=-YOLUDILX) 


RHOd = RAVU( UH*JH*tVH*VH) 

PSITJULU(LA) = PSITULO(LX) +RHOA* (UN*®DY-VH*)DX) 

LF (ABSCALPH@ALPAP).LT.PL) GO TO 15 

Wwe HAVE JUST CROSSED LIVe JOINING EXLT VELOCITY '0 FRALLINo LUGE 
PRINI 964X1(K5) 

FURMAT (* CRUSSING CUT AL *2E15.4) 

FAC = TPI Ra SLON( SCALE» ALPHP) 

GO TU 5 

LF (ABS( THE TA=TAHETAP).LTePI) GO 10 20 

we HAVE YUST CROSSED THe CUTs ADYUST Pol 

FAC = SIGN(2.*PI*SCALEsTHeTAP) 

PSIUOLU(LX) =PSLILO(LX) +F AC ®(OPSI-KHOH* (UH* Y3IMG-VAFEXBTMHG) ) 
XOLJ(LX) = XNtwW 

YOLU(LX) = YNew 


CUNTLNUE 
YULU = UNEW 
VULD = VNtu 


fe (Est teLL) GU 59 10 

TEMP var XL4K3) 

LF (KK.GI.G) TEMP = CMPLXA(REAL(UXI(K3) ) eREALIMETAICL J) 4¢5LTGN( 10008 
AIMAG(xX1L(K3)))) 

“ = NN¢l-l] 

LF (LZEGeII) WRITE (N194U) KKeM 

Lr (MODE«GT.2) GO TO 30 

WRITE (N1946) JOLIeVULDe (AOLU(U) ep YOLI(Y) sPSITOLD(Y) sJ=1eNX) os Top 

60 TU 10 

WRITE (N194O) (PSTOLO( UJ) eI = LeNX) ewGHl (tL) oe TeMP 


CONTINUE 

ALPHP = ALPH 
IHETAP = THETA 
METURK\ 

FURYAT (4020) 
tND 


SUBRUUTINtE GETPSLILSIGNedJ) 

REAL “ACH 

COMMUN /SAZ GAMMAsMACH 

CUMMON /C/ XAREALeXAIMG e TAREAL se YALMG se XBIMGe YBIMG 

COMPLEX UeVeXl ae K2eXSeVlLeV2eVSeL AMUAP,LAMVDAY XT ec lA 

COMMON X1L( 209) eWGHT (200) ee TA( 200) eU( 200) eV( 200) eK1( 200) +Y1( 200), 


1 X2(200) *¥e( 200) +X35(10e2UU) eYS(1U0200) sLAMIAP(20U) eLAMUAN( EDU) 


COMMON /DZ PLleAAsMRPeNRNeNXeNCReL OMe IPLTeMIDE NB eNCoNivet le ude KK 


60 


10 


20 


COMMUN /F/ XURLGeYURLGySLALE ePSILeUOLIJ se VOLIsXOLUI30), YULUIS0)» 
1 PSIULD(3G6) 

CUMPLEX KLTAlL s+ S3PePOLE se XLAgKLBvAIL 

COMMUN /H/ BP(10) ePOLE se XLAe XL BeXAC eo KL TAIL se PLISIZeNoOPS : 
CUMPLEX CUNEWeCVNEWeCXNEW Se CYNEW UX DY UA V4eCSeCXHVeCLUGl oleceCLVU 
1 sCUULDsCVCLDeCXOLUeCYOLU+eLOG1+LUG2 

DATA tLe AG0 el thks ce FSW 707 

COMPUTE IHE COMPLEX LOG UVF Zz AND CHOOSE THe BRANCH wWHOSt ARKSUMENT 
KS .CLOSEST TO TH . 

CLOG2Z(Ze!H) = CMPLX(Oee IH) +CLOG(Z*¥CMPLA(CUS( TH) s=SLN(TH))) 

Ley CLSWeNE. 50); (GO TO sO 

4sw =1 

GAM = (GAMMAr1.)/2-6 

COSY =1e/(MACH*MACH) + GAY 

VAMIL = 1-4/(GAMA-1.) 

RHOINF = (MACH*MACH) €*GAMIL 

LF (LSIGN.GTeO) GO TU 15 

LUG1 = (UeeDe) 

LOGS = (UeeDe) 

L = IABS(LSIGN) 

CUNEW = U(L) 

CVNEW = VIL) 

LUG1 = CLOGC2(X1L(J)-X1(1) sAIMAG(LOG1)) 

LUG2 = CLUGZIXL(J)-XLINS) sAIMAL(LUGE) ) 

CXNEwW= SCALE*((X1(L)+1#xALMG) €LUG1+(X2(L) tL ¥eXBIMGO) ELOGStKS(1eL) 
1 -XORIG) 

CYNEWS SCALE (CYL (L) +1 YALYG) *LUGI+(Y2(L) tLe YBIMo) *LOLA+T3 (19) 
1 “YORKIG) 

LF (LSIGN.eLT«O) GO 10 20 

JH e5*(CUOLUFCUNEW) 

VH o5*(CVOLUtFCV NEW) 

UX CXNCw-CXULO 

uT CYNew-CYULQ 

Ls COSY = GAY ( JH*UH4+VH*VH) 

CRHO = RHOINF *e€CEXP(GAMI*+CLOGI(CS)) 

PSIUJLU(1) = PSLOLD(1)+CRHI*(UH*®DT VHX) 

CUOLU CUNEW 

CVOLU CVNEW 

LXOLU CXNtw 

CYOLU LYNEW 

LF (CULT oe Jd ce UR (LeLTe14)) RETURN 

JOLY = REAL (CUULDL) 

VULU = REALICVOLO) 

XULU(1) = REALICXOLD) 

YULD(1) REAL(CYILO) 

Xe TURN 

tO 


FUNCTION RHO(QS) 

CUMPUTE VENSIIY KHO AS A FUNCTIUN UF THE SaJARE YF THe SPecOyeS 
KeAL MACH 

CUMMON SAS GAM4A, MACH 
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Ce Sot earl ae a 


10 


100 


LUGICAL LSw 

JATA ISW /,FALSE,/ 

tF (LSwW) GO TO 10 

LSWie = sURUES 

GAM = (GAMYA-1e)/26 

CUSQ =1./7(4%ACH*4MACH) + GAY 
GAML = 1e/(GAMMA=1.-) 
RHOINF = (MACH*MACH) «*GAML 


CS = COSU=GAM*US 

RHO = RHOI\F*CS*«*5AML 
MC TURN 

CVU 


SUBRUUTINE GETXYO (K1eK2) 

GET UeVeAl ANU ¥1 ON XI=XIA FROM EIA(K1) 10 ETA(AK2) 
COMPLEX UeVeXlLeX2eXZeVYL eV2eV¥SeLAMUAP eLAMDATsXITecla 
COMMUN X1L( 200) ewGHl (200) ec TAI(200) sU(20U) ev (200) +AL(2NU) +¥1(200), 
1 X2(200) +V¥2(200) 6X3(10.2UU0) ,YS3(10 e200) sLAYMVAP(20U) sLAMUAM(cOU) 
COMMON SCY YXYAREAL s+ XAIMGs YAREALe YALMG 

COMPLEX LeSe Te LAM3PHeRH1eRkH2eSU0F XL 

DATA T° 70.9 le) 7 

COMPUTE SAT XL t= xia 

> = SUFXL(XI(1)) 

CALL LAMBDA (U(K1-1)eV(K1-1) sLAMDAP(K1-1),LAMDAM(A1l=1) ) 
JO 100 L = K1leK2 

COMPULE ICETACEy) 

1 = CONJUG(SOFAL(XI(L))) 

EFACL) = CONJCGCXI(L)) 

JC) =Ss"uUCk =) 

V(E) = vit=2) 

CALL GETYUV(SeTeJd(L) eVI(L)) 

CALL LAMBUA(U(L) eV(L) eLAMUAP(L) eLAMDAM(L)) 

FINU X17 AND Y1 ALONG XI=KIA AT ETAL) 

CUMPUTE LAMODA= AT MIULPOILNI 

LAM3SPH = eSe(LAMDARP(L)+LA“%DAP(L=1)) 

RH1 = X1l(bL-1)+LAM3ZPH*Y1(L-1) 

XHOQ =-1*(XAIMG+LAYDAM(L) ®YAIMG) 


Y1(L) = (RHI RAZ) /(LAMBPH=LAMUDAM(L) ) 
X1(L) = KH1 = LAM3BPH*Y1(L) 

CUNT LNUE 

RETURN 

LNO 


SUBRUUTINE GEIXY1 (K1) 

bET UrVeLAMOA+e LAMDA=eX1 AND X2@ AL XI(K1) e«CONJGIALA) AND STORE 
IHe RESULTS IN THe FIRST ELeMeNT OF Toe RESPECTIVE ARRAYS 

IHE FIRST ELEMENT OF THe RESPECTIVE ARRAYS WILL CONTAIN IHcIr 

VALUES AT xJI](K1l-1)+CONJG(XIA) 

COMPLEX VUaVeXlL oe X2eXSeV1Le Vee VSeLAMUAPeLAMDAT XLT ectlA 
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LCUMMUN X1( 200) swoAl (200) ec TAC 200) eU(20U) .V( 200) sA1L1(2NU)eT1(2U0), 
1 X2(200) eY2(20U) eX5(10,2UU) -Y3(10+20U) sLAMIAP(2UU) sLAMDAM( eCU) 
LUMMUN /LC/ XARLCAL se XALMGe YAREALs YALMG 

COMPLEX LeSeleXHdleRKRH2+LAM3M,LAMBMHe SOFXI 

VATA te™ 7° CO! Nel Lier) 7 

CUMPUTE 3S ANU XI(41) 

S > 0F x POx TOKE) 

COMPUTE T AT tfiA = CONJGIAILA) 

§ = CONJGO(SOFXL(XI(1))) 

CALL GETUV(SeTeU(1)eV(1)) 

CALL LAMBUA(U(1) +V(1) eLAMOAP(1) «+L AMBM) 

PING XL ANDO YIsAT ETA =. CONUGUXEAN At IXTGCSD) 
COMPUTE LA*“O0OA+ AT MIODOPOLNI 

LAMSMH = 2«&4(LAMUA(1)+4LA%8™) 

RH1 = X1(1)4LA%BM4A*Y1(1) 

HHO = L#(XAIMG+LAMUVAP(1)*1TA1MG) 

Y¥1(1) = (RH1"RH2)/(LAMBYH=-LAMOAP(1)) 

X1(1) = RHI = LAMBYH*EY1(1) 

a1Ort NEw LAMUA+ 

LAMUVAM(1) = LA“3SM 

RETURN 

LNO 


SUBRUUTINE GELTX2(41eK2) 

FINY x2 ANG Y2 ALONG XI = XIdB FRUM CTAIKi) [0 EIAIKA2) 

LUMPLEX UeVeXLeX2eX5eVlieVeeVSesLAMVUAP eLAMUAM ALS CIA 

COMMON XL( 200) eWGHT (200) sc TAC20U) eU(20U) eV 200) eAL(2NU) se ¥1(2U0), 
1 X2(200) eYZ(20U) oX3(10e20U) eY3(10e200) sLAYMIAP(20U) eLAMUAN(20U) 
COMMON /U/ PL eANeMRPeNRNeNx eNCReLBEMeIPLT MODE NG e NC eNNe LL ede KK 
COMMUN /C/ XAREAL se XALMO se YAREAL se YALMGe XBIMor YBIMGsITPTR 

COMPLE X KH] eKH2eLAMBPH,I 

OATA T° SC ee Lael7 

CHECK TO SEE LF WE ARE AT xi 

te (K1.G). (NB+1)) GO TO 5 

FIND x2 AND Ye AT ETA = CONJG(XI 8) 

XH] = 1.-/7LAMVUAM(NS) 


K2(NKH) =-(YBLIMG+EREAL( RHI) *X53LMG)/SALMAG( RHI) 
YT2( NH) = “KHL*(X2S(NB) F1LexXSIVMG) e-LeYolMo 
CONTINUE 


JU 1U L = Ki efKe 
CUMPUTE LAYOA]= AT MIUPUILNI 


LAMBRPH = eEe(LAAMUAP(L)+4LAMUAP(L=1)) 
MH] = X2lL-1)+LAM3PH*Y2(L-1) 

442 =I *#(XBIMG+LAMVAM(LI*TGIMG) 
Y2(L) = (Ree RN2)/( LAMBPHH-LAMDAM(L) ) 
X2(L) = RH] = LAMSPHEeYV2(L) 

CUNT LNUE 

KReTJUNW 

cNO 
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SUBRYVUUTINE GEIIV (Sel eUev) 

bilyEN S ANOQ T FIND U AND v 

LCUMMONSGZ NieN2@eNGeN7TeNSe NG 41 

CUMPLEX SeTeUeVeGSeWeWST Re] eHSdeHSUPR,YUSNicWeSTe COUR 
CATERNAL CSORI 

VATA EF sAcC0se hel 7 9 TOL Leb =107 


ws = U*UtVEV 
We-S8 

WSTR = T 

oh = (Swe 


Ke (CASS (ST)..6E .TOL) ‘GO-T0 (35 
CALL GETHSQ((Uee00.) tHSQeHSuUPR) 
w= W/SURKT (REAL (HSQPR) ) 
WSTR = WSTR/SART(REAL(HSUPR) ) 
60 [0 100 
VO AT MUST 20 NEWTON ITcRATIUNS 
25 vO 8U L = 1,20 
LOMPUTE H(GS)**2 AND ITS JVERLVATIVE wllH XcSPECI TO QSe HSSPR 
CALL GEIHSG(USr+HSQeHSUPR) 
ASNEW = US-(HOUHKST)/HSQPR 
LF (CABS(USNEW-QS) LIeTUL) GU TO 90 
80 YS = ASNEW 
WRITE (N2e110) SI 
CALL EXII 
90 as = VSNtw 
W = CSQRI(GS*(S/T) oW) 
ASTR = wel 1/5) 
100 YV = e5*(W+twSIR) 


Vv = eO*]4#(WeWSTR) 
Xe TURN 

110 FYURMAT (42H NEWTON LTERAILUN DID NYUI CUNVERGE AT SI =,2F12.4) 
cNO 


SUBRXYUUTINE GETHSU( US eHSeeAHSGPR) 

CUM-UTE H(GS)**2 AND ITS DERIVATIVE WIId RESPECI TO 9S HSAPR 
Xe AL ‘ACH 

CUMMUN SNS GAMMAeMACHeRRAROGI 

COMPLEX US eHSdeHSAPReyBPHI 1, BPH] SeCSGRT se ROUT 2, TEMP eS ¥ eCSeRIOI1LeRR 
cXTcRVAL CSQRI 

LUGICAL 1Sw 

VATA ILSw /,FALSE./ « CONSI /1./ 

He CLS Wi) eGo VEO: 10 

GAM = (GAMMNAH1e)/26 

LUSQ =1./7(MACH*4MACH) + GA 

ASTRSOG =(2,./(GAMMA+1.))*CUSY 


CMU = SORT ( (GAM4A=1.)/(GAYMMAt+1.)) 
FAC = EMU/QGSTRSYQ 
ISX = GS 


WS =) (14-0/0'5,) 

10 CS = COSUSGAM¥eUS 
SPHI1l = (CS#CS)/NUSTRSQel1. 
PHILS = (C5+#C9)-05 
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noQarn 


20 


10 


XUOT2 = CSGERI((CS#CS)*(3PHI2-QS) e«ROOTL) 
KRUOTL = FAC#RUUT2) 

(eMP = CEXP(CLYG(3PH1L1-RUUT]L)/EMU) 

1EMP =CONSI/(IECMP*e(BPH12+RO00!2)) 

IeMP 1S n*€ae/US 

HS@ = wS* TEMP 

LE CLS ede-COx TO20 

4SW = TRUE. 

Sem GCONS! SO THAT HEL) = 2 

CONSTI = 1./REAL(HASQ) 

as = OSX 

6U TU 10 

COMPUItE ITHE UtRIVATIVE OF HSQ@ wlIH REsPECI TO Qs 
HSQ@PR = IEMP*RUOT2/(CS4+C3) 

RETURN 

tND 


SUBRUUTINE LAMBDA (ULeVEeLAMUAP s+ LAMDAM) 
COMPUTES LAMBUA+ AND LAMbUA]- GIVEN U AND V 
KReAL MACH 

LCOUMMUN /A/ GAMMAs“ACHeROUI 

CUMPLEX ULyVL+eLAMDAP sLAMDAMsVeCSsQSyUel SURI KOU! 
EXTCRKNAL CSGRI 

LOGICAL 154 

VATA 1SW /.FALSE./ 

LEVI ISW) GO? TOeg0 

iSwW = , TRUE, 

OAM = (GAMMAele)/26 

LUSO =1./(MACH*MACH) + GAM 

Bis OL 

Veise VL 

dS = J#U+VOV 

CS: =) COSQ: = GAN*«QS 

RUOT = COURTICS#€(aS-CS)eRkU0T) 

LAMOAP = (USV4#R0U0T)/(CS-VEV) 

LAMUAM = (U*V=ROOT)I/S(CS-VEV) 

RETURN 

tNo 


SUBRUUTIWE GIPATH (K1eK2943) 

KEAJS IN THE PATHS ON THe COMPLEA CHARACTERISTICS AND ScIS YP THE 
ANITIAL ORILD 

PATHS Akt PRKESLCALSZEU ON JHE PLANE KISLONJO(ALA) ANU TAce 

bRID IS SET UP IN THE PLANE ETA=ALA 

COMPLEX UeVeXlLe XoeXSeVle VereVYSeLAMUAP,LAMUAMe XLT ecia 

COMMUN XL (200) eWGAT (200) ee TAC 20U) eU(20U) eVI 200) eAL( 200) eTL(2uU)s 
1 X2( 200) eV¥E(20U) eX3(10e20U) pV S( LU 200) eLAMIAP(20U) sLAMUAM(eDU) 
COMMUN /U/ PI eANe ARP sNRNeNK eNCReL OMe IPLT se MIDE Noe VC eNVe bbe Je AK 
LUMMUN/SG/S Ni *N2eNGeN7eN3Se VGe M1 
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10 


20 


30 


60 


1 


LCUMPLEX K1I1WALLegPePOLE se XLAeXL Be KLE 

COMMON SHS BP(1G) ePOLE eXtbeXLBeXLCeXITAIL ePLISIZsNBPS 

RKRrAL MACHA, MACH se ATALL 

LCUMMON /K/ RNePCHUePCHL e MACHA sc ANGLA,MIFAIL, ANGLT,» YWACHR,ANGLS, RC, 
ANGLNe IR 

CUMPLEX PUINIyH,yXILOFS 

JATA KAT 4.017 9 K2MAX SUS 

READ (N6e1N0) NCAROS KK 

**ee4CHECK FOR cNUO OF FILE KKe* 

lr (tUF(NG).NE20) GO TO 8O 

ICARUS = IABS(NCAROS) 

LF ((K3Z.CO,2) ¢-ANDe(MUDE.GI.LO)) WRITE (4-140) NCARUS,AK 

KSwITCH = KK 

A = Klel 

OU ves) KK 

ASwWICH = O 

LSOw = MOvUe 

Lr L(MOVDERMEG SS) eAND ec (MTALL CLE. 0. )) ESwW -= 0 

VO 60 L = 1,NCARVS 

READ (N6e110) POINT swGHTA 

LL = WGHIX 

JSeEe XJIALL IF MIAIL HAS BEEN ReAU IN 

LF ((LSW eEGe2) eANUe(L2-EXeNCARDS)) POINT = XITAILL 

tr ((K3.08,2) eANUe (MUDEGI1.L0)) WRETeE (NG4,150) PULNTeLL 

AUJJUSI PATH TU THE NOSE SO THAI IT ENDS AI A STAONATION PUINI 

LF ( (MODE eEG 21) eANDe (LeEGe NCARDS)) POINI = XIOFS( (Use. ePUINI) 

Le MKSWICH Gilke.0)316.0> TOs20 

Xen) 

CALL ADUYUST( IXeIXePOILNT) 

Le -CEX cee) GO TO 26 

POINT = POINI + RATH(XI(K)-POINI) 


60 [UV 10 
4“ = CABS(POLNI=-XI(K))/sAAt+,999 
AN = Kl 


Le (CEE SEQs6 0) EL S12 

MM = MMaAMRP 

LTE (lsSwiCH.t les) 4M = AM*eLASS( LE) 
LF (LL.GI.G) GU TO 56 

JJ = KOM 

LSwTCH = MGDE 

KSwllH = ej 

Le. set 

IF (KA.LGI.C) CALL ADJUST (1e1ePO0INT) 
Ho = (POINT=X 20K) ) FLOAT (44) 

A = K+MM 

VU 60 J = KNek 

WOHI(J) = FLOAT(LL) *FLOAI(LL) 
AL(J) = XI(JU-1)4H 

LE’ CLSWECHCET eS) GO 10 65 

ir CAV EE. dd) oJ TO 695 

AN = JJUt2 

1F (KN,.GtEeK) GUY TU 64 

VU 62 J = KNetk 

IF (WGHT (JU) Ede wGHI1(J-1)) GO TO 62 
WOHT(J-1) = ¢3*(WOHT(J=1)+WwGHI (J) ) 
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62 CUNTL2NUc 


64 WGHI(JJ) = ,S*€WGOHT(JUt1) 
WOHT(K) = ,5*#WOHT(n) 
65 K2 = A 


CALL NOIWYUST (JdeK2eX1) 
K2MAX = MAXO(KEMAXKeK2) 
LF (LBMebE,KeEMAX) KETURN 
WRITIE (N20129) K2eLBM 
CALL EXI1I 
80 “ODE = 99 
JJ = K2“MAX 
KETJURN 
100 FORMAT (313) 
110 FORMAT (26€20.12.F10.5) 
120 rtUKMAI (/735H Ide GRIO IS TOU LARGE FUR PRIGRAMe43903H K=l4) 
PEO “RORMAT -U/S1 KX eA) 
150 rURKMAT (51% e2r4.3014) 
tNO 


SUBROUTINE AUVIYST (K1leK2eXK]) 

Xe AL MACH 

LUMMON /AZ GAM“Ay4ACH 

COMMON /U/ PI eAAeMRPeNRKNeNXeNCR ALAM ITPLIT eAMODE NG eNCoNive Ll edd enk 
CUMPLEX AL(1L)eSeSCRLTeXLUFS,SVOFXS 

LUGICAL ILdsw 

VATA LSW /.FALSEe/ 

fF (1SW) Go TO § 

L5w = TRUE, 


YERLI = (GAMMAH1L.)S(GAMMAt1.) + 2e/( (GAMMA41 ©) *MACHEYALCH) 
LALL GETHSG (CMPLX(QCRIT*#U.) eSCRLTeS) 
CONS! = REALCSCRIT) 


VO Aid = KiVKe 
S = SOFXL(xItU)) 
AdSSU = S*CONUG(S) 
IF (KK,LE.0) GU TO 20 
SUPEKSONIC PATH, PLALE PULNT OF THe SOWIC LOCUS 


SCRII = S¥*¥SURTI(CUONST/ABSSa) 
A1(J) = XLOFSC(SCRIT,AI(J)) 
oY TU 10 


“AKe SUR THAI THE PATH KeMAILNS SUBSONIC 
20 LF (ABSSHYeGTeCUNST) GO TU 3¢ 
10 CUNTILNUE 

AcTURY 
30 “2 = Jel 

Ke TUKN 

CNO 


OaNnNAN 


20 


10 


40 


399 
30 


SUBROUTINE INLTFN (ETAgZeL¥) 

CUMPUTES IRE LNITLAL CHARACFERISTIC DAIA UN XL = ALC 
IHIS ROUIILNE CAN 3t REPLACED BY ANY OTHER ANALYILC FUWCILUN 
ETA 0S. THE VAGUE OF ETA 

2 IS THe VALUe OF THe ANZ2ITIAL FUNCTIONe G( LTA) 
COMPLCLX ETAeZeXIT Pele TEMet TU 

Xr AL LOGS 

LUMMON /t/ LOGS (4o12) eETIU( 4) eXReNLUGS NAKA pNP NUL (32) 
VIMtNSION PARAM (56) 

CQULVALENCE (PARAY(1) eLOGS(491)) 

DATA To FCO ne hie) 7 

LX = LU+NJ 

PRE EXC EUS) 60 [2 5 

KA = ESCEX) 

PARAM(KX) = le 

XLP = CONUG(ETA) 

SUM UP THe PULYNOYILAL TER“S 


ai 

f= ESC!) 

Js Jel 

2 = iZ¥XIP#tETUCJ) 


Le iCu.Gilel) GUO TO 26 

22s 2*xXIP 

LOG TERMS 

JU 10 J = 1,NLU0G6S 

LF ( (LOGS (1,0) EG Oe) -ANDe (LUGS(204).E0-0.)) GO TY 10 
lemMP = Le = XLPSCMPLX(LUGS(SeU) eLUGS(4eJ)) 
Zs Z+CMPLyx(LObS(1LeJ) eLUGS(2eJU)) *CLOG( TEMP) 
CUNTINJE 

Z = CONJG(Z) 

dF (LX,EFe1) RETURN 

PARAM(KX) =*U. 

Re TUKN 

CNO 


CUMPLEX FUNCTLUN XLOFS(SeXIPAST) 

CUMPLEX XIPAST ee XITOLUe XL NEW SUF XL ea SUFXIP eS 
VATA TOS Jlee=207 

SLOLU 1S XTPASA 

VO AT MOST 20 NewlON IIERATIUNS 

JO 4U K = 1,420 

AlWtW = XLOLUV-(SOFXI(XLOILI)=-S)SSOFALIP(ALOLD) 
LF (CABS(XINEW -XIOLD).tt.fOl) GO 10 50 
XIOLU = XINEW 
CONTINUE 
PRINI 99 
FURMATL (* AO 
XLOFS = XINEW 
Re TURN 

tNy 


CUNVEKGCENCE*) 
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30 


40 


50 
39 


60 
10 


20 


10 


20 


LOUMPLEX FUNCTLUN SUFXL (Al) 
CUMPLEX A1,TEMPyPOLI,C(12) 
COMPLEX XITALL oe dPePOLE se XLNe XI Be ALC 


CUMMUN SHS BPI10) ePOLE se XLAe ALB eKALUeXALFALL se PLISI 2+ NoPS 


LUGILAL LSweNPJULc 

VATA ISW /,FALSE./ « NPULC/,FALSEo/ 
1F (LSwW) GO FU 10 

1Sw = TRUE. 


LF (( REAL (POLL) NE COe) eORe (ALMAG(PULE) NED.) ) 


POET ="0¢6 
NDEG = NBPS + 1 
dr (eNOT.-NPOLt) GO TU 30 


NOEG = NBPS 

POLIT =91e/P0LE 
C(1) — it Gs 

C2)" =" BP (L)*EP1) 
&(3) = -26 


tr (NBPS.EG,1) GU TO 55 
VO 50 J = @yNBPS 


VEMP >=  =1,78P (J) 

C(JU+2) = TEMP*U(J41) 

sa CY 

C(K+tL) = C(KI*#TEMP+C(K41) 
ACs hej 

1F (K.2GTe1) GU TO 40 
CUNTINUE 


Clo) = Cl2)=POLE 
UU 60 J = 2@_NULG 


COS+1) = (C(CU+1L)+FLOAT(J-2) «C (JU) «POLL)/FLIUAT(J) 


1EMP = C(NGEGF41) 

VU 20 J = 1ieNUtES 

1eMP = XL*¥TEMP+C(NDEG+1-U) 

SUFKL = IEMP 

Ler (NPOLE) SOFAL = SOPXIACL =POLLexXs) 
KETURN 

tNo 


LCUMPLEX FUNCTIUY SUFAIP(AL) 
CUMPLEX ALITAI Le SPePOLE eXALAeXI Bre ALC 


COMMUN S/H BPC LU) ePOLE sXLAcXL Be XLCeXLTAIL ePLISI2Ze NSPS 


CUMPLEX KI,TEMP,POLI 

LUGICAL 1SweNPULe 

YATA ISW /.FALSEe/ « NPOLC/S. FALSE e/ 
tr (1Sw) GO FU 10 

lsw = eITRUE, 


AF (CREAL( POLE) NEO.) ceORe (AIMAG( PULL) - VEU.) ) 


Le (NPOLt) POLA = 1,./POLte 
VUES = NBPS 

feMP = devl(BP(1)-XI) 

LF (NOEG.LE.1) GU 10 30 
V9 20 J = 2@eNULS 

1emMP = ITEMPalle-KI/BPIJ)) 
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NPOLE 


VPUOLE 


elRUE. 


el RUE. 


30 SOFXIP = TEMP 
tr (NPOLE) SUFXIP= TEMPs(1,-POLI*XL#¥(2e-PULI*X1)) 
Ke TURN 
CND 


CUMPLEX FUNCILUN XKLOFM( “Me 4NGeX1LULU) 
CUMPUTE XI CLUSEST TU XLULO WHICH HAS MACH VUMBeX M AND ANOLE ANG 
LUMMON SAS CAMMAd4MACH 
Xe Au MeMACH 
LCOMPLEX HSU,HSUPR,S,AIOFSeXIOLU 
LUGICAL ISw 
VATA PI £43.14199259358979/ .« ISW /.FALSZ,/ 
4F (LSw) GG TV 10 
iSw = .INUE, 
RAD = PI/1EGe 
SAM = (GAMMA=1.)/26 
LUSY =1e/(MACH*EMACH) + GA 
10 EMSQ = Mem 
ANGL = ANbs#RAVU 
WS = COSQU*EYSI/(1.4+GAM*EMSQ) 
CALL GETHSG(CYPLX (QS 00.) *ASQ+HSUPR) 
> = SURTUIREAL (ASQ) )¥CMPLX(COS(ANGL) »=-SIN(ANGL) ) 
XLOFM = XIUFS(SeALOLU) 
RETURN 
ND 


SUBRUUTINE TIIce 

PRINUS tHe I2L1LeE PAGE OF AHE VUTPUT 

REAL MACHyLUGS 

CUMMON SAS GAA’ “ACH 

CUMPLEX UeVeXleX2eXZeV1leVeeVYSeLAMUAP,LAYMVDAM, XI eLIA 

COMMON XL0 ECU) eWGHI (200) ec TAI 200) eU( 200) eVI 200) eAL(2NU) pe VY1(2UU), 
1 X2(200) eY2(2UU) eX3(10e2UU) Y3(1Ue200) eLAYIAP(20U) eLAMDAM( 2G6uU) 
CUMMON JUS PLeAAs4RPeNRNVONXeNCReL GM LPLTe VWIIEeNGeNCeoNNe LL eS O KK 
COMMON /t/ PARAYM( 56) eXReNLOGS eNKeNPoeNIelLC( 52) 

LCOMMUNSGS NieN2eNGeN7eNZe V4 QML 

CUMPLEX XAXLTALL eBPePOLE se XLACALB eA XILC 

COMYUN /H/ BP(10) eFPOLE se XLAeXLBeAXALCe XLT TAIL ePLTSLZsNBPS 

XtAL “MACHA, “MACH3YIAILL 

LCUMMUN SKS RNePCHIU + PCHL MACHA se ANGLAMTALL eANGLT se MACHR eANGL3ekCe 
1 ANGLNe TR 

VIMENSION P1 (356) eLX(56) eLOGS(4e12) 

CUYUIVALENCE (PARAM(1)eL065(101)) 

RtwlNO N1 

Le CrP 0). SO TO 7 

READ (N1¢5G) (PARAM(L), & = 1456) 

Jr (N1.EU.M1) GO TO 7 

RewLNo M1 

Kr AD (N1+560) (P1(L)-. L = 1456) 
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on 


10 


12 


14 


vO 5 LC = 
PARAM(L) 


156 


= (4e*PARAM(L)=P1L(L)I)/S6 
VATELLS ACCOCUROUTING WHICH (GLVES ChHE 


CALL CAILC(IGATE) 


VMLeNRNeN“M2eMRP 
NMLeNPeNM2eNCR 


NMLePLTSIZeNM2eIPLi 


NRUN = IABS(NRN) 
WRITE (N4eE0) NWR 
VM1 = GH RUN = 
NWO = BH MKP = 
ARTTe (N4e60) 
Ny = 8H NP =s 
VWM2 = BH NCR = 
WRITE (N4e60) 
Ny = BHPLISIZ 
N“2 = 6H TRI 
WRITE (N4e1G0) 
N“W1 = 8H MACH 
YWM2 = BH RN 
Ber="2nes 

tr (KNX.EG.00¢) EC 


KVX = 1etegekN 


YM1leMACH ea NM2eRNXet6 


NMLePCHUe VMS ePCHL 


VMlLeRCeNM2sANGLN 


NYLeMTAILeN“%2eANGLI 


NMleTReNM2eXKR 


NML se GAMMAeNY2¢AA 


NM1leMACHB «N42, ANGLE 


VM1leMACHAeNMS +e ANGLA 


WRITE (N4eS0) 

tr (RMN.tY.U.) 90 
NM] = BH IRANU = 
NM2 = BH TRANL = 
WRITE (N4%*906) 
V“4y = BH RE aes 
VM2 = @HANGNOSt= 
WAI Te (N4e90) 
NAte=(6H APALL &s 
VM2 = SHANGTALL= 
WRIT Te (N4e90) 
V4M1 = BH TR = 
NMO = 6H wks is 
WRITte (N4e90) 
NM1 = BH GAMMA = 
NM2 = BH GRIU = 
WRIFE (N4e9G0) 

LF (MACHSBeLTeU) 
NM = BHMINLEF = 
WM2 = @H ANGIN = 
ANXITTe (N4Ye 90) 

LF (MACHA.LT+0) 
V4] = BHMEXTT = 
NM2 = BHANGCEXIT= 
WKI Te (N4e90) 
WRITe (N4970) 
N41 = BH xfA = 
WM2 = BH Kip += 
WAITE (N4e913U) N 
NM] = 8H xIiC.<s 
W442 = &H BP = 
Avece 

are <i 

A 1 


BH 


POL 


UNe IVATE 


esa th 


TO 10 


69 TU 12 


GY TO 14 


MLekXlAr NVM e X1U 


MONITAeOAYs 


REAL (POLE) «ee Qe Ne) eANDe (ALMAGI POLE) eEde Oe)? 


29VU 


ANU 


| n 


fcAR 


15 


WRYETe (N4e130) NMLeXl1Ce NMS ePULE 
vO TU 18 
16 WRIT (N4e130) NMLeX1CeVM2e8P (1) 
1A LF (NBPS.GE.K) WRITE (N40130) (NM2eBP(U) ed = Ke NSPS) 
WRITe (N4e70) 
VY 2U L = 1.56 
20 UX(L) = 4H ’ 
VO 3U L = 1,NK 
J SALE 5) 
30 LX(J) = 4H # ¢ 
WRITt OUI POLYNJUMIAL COEFFICIENTS 
NM) = 6H tl = 
N“2 = BH ETo@ = 
WRI Te (N4e110) NMLePARAYM( 49) eLxX0 49) ePAKAM(30) eLX090) eNMerPARANI51) 
1 #£LX(51) ePARAM(52) eLal5e) 
NM] = @H ETCULS 
N“o = BH ETFU = 
WRITE (N4e110) VM1ePARAYM(53) eb X(95) ePAKAM( D4) gL X034) ep NMZePARAYM(59) 
1 9LX(55) ePARAM(56)eLX(56) 
PRIVEE OUI LOG CUEFFICIENTS AND LOCATIONS 
dr (NLOGS.~EQ-0) GI TU 45 
VU 40 J = 1eNLUGS 
WRITE (N4e140) LUGS( Led) eh X(4# 593) eLOGS (295) oe L XC F#IHF) eLUGS (Std) 0 
a LX(G*eI-1) LUGS(4eJd) 
4O CONTINUE 
45 WRITE (N4_1S50U) 
CALL PLPATH(2+4AX0(OeNP) 02) 
LF (LPLTeNE.U0) CALL GOPLUT(NRN) 
Xe TJKN 
60 FURYMAT(S1IT7XsAB re L4 el 3ZX1LHi IAA T4) 
70 FURMAT (1H) 
BO FURMAT (1H1 /15xXe28HITRANSUNIC CASCADE UVESISN RUN 1L5917XA10///38X_ 
1 6HIAPE 7 7) 
JO FORMAT (SLTX «ABV FT eS elLOXIHit Dx eABre Tere) 
100 FURMAT(/17X AB FT. SeLOXLHIEIXe AB,lY ) 
110 rOR“AT ( EXABeF Te SaeNMbe Fle 3eAh2e OXLHs 4YXABeF Te SeA4eFT7e3eA2s) 
150 FORMAT ( EXABeE7 a SedXlHe EF (SSe 6X2 4XeABsPh7e3eSK1LH, FILS) 
14M FURMAT (/11X5AL0G =F7eSrAteF Te SehleGX1lNeleXF7 eS eA4eF 7.5) 
150 FORMAT ( //3UX25H* AUTOMATED PARAMETCRS ) 
50 FURMAIL (40206) 
tNp 


SUBRUUTINE BUUY?T (MX) 

KEAUS TAPE DATA FINOS BUNDY POUINIS AND CONSTRUCIS PeN DIAGRAWS 
CUMPLEX X1,X180DY+CHAR 

Kr AL MACHN, KAPPA 

COMMON X1(2UU) sUREAL(20U) »VREAL( 200) eAREAL( 200) *YREAL( 200) PSI ( 
1 2UU) eXLBODY (200) eUSODY (200) eVBUUY(20U) sXASO0UY(2U0) eYSUUT(20U) 
2 CHAK (41010) sLLINE( 200) +5( 200) ee MACHN(2U0) eX APPA( 200) ,ANGL(20U) 

3 eTHeTA(200) sH(2U0) 

COMMON /LC/ KAREAL se XAIMG se YAREAL se YALMGe XBIMGe YBIMGe FPTR 

COMMON /US PleANeARPeNRVe NX eNCReL OMe LPL Te VIDE CNB eNCeNNe Li geJSy AK 


T2 


LCUMMUN/SG/S NiwsN2yeN5eN7eNZe NG QML 

CUMMON /JS POI405) eXMINe IMINO XFACe YTFACONCHe*LSZeNLINe TWUeUL Ve ACH 
COMPLEX ROTATE 

INTEGER HP,HN 

VATA HPeHNeNRUW/1APeLHNeL/ 

ANGROT = ATAN2(-X31MG+-YB1IMG) 

RUTAITE = CMPLAICUS(ANGROI)»,SIN(ANGROT)) 
SET UP AXES FOR THE PeN UIAGRAM 

CALL PNAXIS((U6ee0-)) 

SET UP PLOT PAGE 

1PL = MOU(IPLI/10.10) 

Sr (2PL2690,6) GO TO 5 

CALL ORIGIN (345015) 

ARITTe (N4e20U0) 

1F (NY.NECN2) WRITE (N2e2U0) 


XSOJY(1) = le 
AX300T(2) = De 
TSIUIT(E) = Oe 


REAU(N1T +210) USOUY(1) eV300Y(1)-.1TBODY(1) 
dr (ECFIN1) .NteeO) RETURN 

REAU (N19210) USUDY(2) eV8UDY(2) exleODr(2) 
Lr (NLsNEeN3) GO. TO 7 

READY (M1-¢210) URS eVR2¢YR2!) 


YVSOUV(1) = (4.*UBIQUY(1)-UR2)/3. 
VBODY(1) = (4e*VBOUY(1)"°VX2)/5.~ 
TSOUT(1) = (4e*YHOUY(1)-YX2)/3- 
REAVU (M1s210) JUReevR2eXI BUDDY (2]) 
UBOIT(2) = (4e*UKBIUY(2)-UR2)/5. 
VBLUDT(2) = (4e*VBIDY(2)-VR2)/3. 
At ae 3 


TIL = YROOY(1)*REAL(KOTAITC) +ALMAGIROTATE) 

ir CVTLekEDe Oe) CALL CPLOI (CMPLA(Oe VIL) e=3) 
READ DATA FOR A PQOTH OR FURA 

ReAVU (N1s9210) KAXeNPTS 

*eeCHECK FOR END OF FILE Kee 

lr (LOF(N1).NWEe0) GO TO 160 

O00: £5 b=  2eNPits 

Ke AU (N1e¢2)0) UREAL(L) eVREAL(L) eXAKEALIL) eYREAL( HE) oe PSLIL) eAT IL) 
CONIINJE 

fF (N1.NE.-N3) GO FO 20 

ReAU (M1210) KAeNPTS 

J = 1 

VO 29 L = 1eNPIS 

Re AU (410210) JURSeVROeXKRere YRSePSIRSeXI(L) 
UREALIL) = (4e *®JREAL( J) -UR2)/3. 


VREAL(L) = (4e#VAZAL (JU) “VN2)/3. 
AREAL(L) = (4e ®XREAL( IU) MARZ) /S5e 
TREAL(L) = (4e *®YREAL (JU) -¥n2)735-6 
PSI(L) = (4.8PSL( J) ePSIR2Z) 43 
we wr? 

CONTINUE 

LN = NPTSe¢ 


CHECK FUR ¢ SUPERSONIC PAIH 
lr (KRA.GIeL) ©Y IDO 40 
at AKE UN A SUDSONIC PAIN UR FURK 
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PUT SYMgUL P AI POINTS OW PATH UR FORK WwHerte POL 15 PUSLIIVE 
PIT SYMBUL V AL POINTS ON PaTH UR FORK WHERE PSL IS NEGATIVE 
vO 30 J = 1eNPIS 
Ll = He 
BA MPSEGd)eiie Ue) Lb = HN 
30 CALL PNC XI (JU) eL) 
LF(XKeLT2¢0) GU TO 10 
we AKt ON fg FORK,LOOK FUR BODY PUINTS 
GO LF (ULN.LT21) GU TQ 100 
VU 90 J = 14LN 
CALL INTERP(Ue4Xe I CHK) 
HE. CLCHKES.0) GO TO 309 
WA BUDY POINT HAS BEEN FOUND IN THe INTERVAL 
CALL BAUPTS(MX) 
CHECK FOR ERKUNEOIS POINIS 
LE CAR. GieG) GU FO 110 
90 CONTINUE 
100 £F (KK.LEe0) GU TD 10 
SUPERSONIC PAIH wIIH NO BODY POINIS 
130 tr (NROWeNE.1) GO 10 60 
MRQ = NPISs65+41 
Lr (N2,.teeN4) FSO TO 60 
PREPARE FOR TTY PAGE OUTPJT 
MRT = NPIS/36+1 
ING =e1 
LF(MRGSNEeMRT) INC = 2 
60 LF (MOD(NRGWe1le4Ra) NEO) GO TO 150 
CONSTRUCT PeN TRIANGLE FOR SUPERSONIC PATH 
NCOUNT = 0 
UV 140 J = LeNPTS eR 


NCOJNI = NCOUNT+¢]1 
L = Hp 
de (PST(JU).LI.G) FC = HN 


1490 iLINeE(NCYUNT) = L 
WRITE (N4e1i80) (ILLNE (JU) ed = Le NCUOUNT) 
LF (COND LUNG) OR. (MUD(NKRUweleMRI) NEU) )GO9 TO 190 
WRITE (N20180) (ILINE( Ue = Le NCOUNT+1INC) 
150 IF (LPL.61.2) CALL. PLTML (NROW.NPTS) 
NROw = NRKOwt+] 
Lee NP TS6 GT.) GO TO 10 
TRIANGLE IS COMPLETE 
WRI ft (N4e.2600) 
LF (NY NteeNe) WRITE (NO@e200) 
NROw = 1 
60 10 10 
PLOT EVERY “MRP CHARACTERISTIC FROM THe 380UY TO THe SONLC GL NeE 
110 12F (CIPLeLE 2) eOR.IMUDINROW=1le MRP) NEU) ) OU FTO 150 
CALL CPLUI (ROTATE *€CMPLX(X309V(MX=-1),YBUDY(4X=-1))¢3) 
JO 120 L = JetN 
120 CALL CPLUIT (ROTATE *&®CMPLX(XAREAL(L+2) eYREAL(L+2) Dee) 
SU FU 130 
160 “xX = Mxel 
WRITE (N4—e290) PG 
{IF (1PL.NE.G) CALL ORIGIN (3251-95) 
Tr (N20EQG6N4Y) KETIRN 
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310 
320 


180 
200 


210 
290 


10 


20 


450 


VU 310 J = 2¢403e13 

JP = J4+6 

WRITE (N2eZ20) (PS(L)eL = JedP) 
FORMAL (1X,8A10) 

RETURN 

FORMAT (65(1X9A1)) 

FORMAT (1H1) 

*#e*eCHANGE (4020) TO (20A4) UN 15M 360*¥** 
FORMAT (4020) 

FURMAL ( TEFPEPFFFF FACILE OLSALOY) 
tNo 


SUBRUUTINE PLIML(41,K2) 

PLOIS CHARACTERISTICS IN IHE SUPERSONIC ReGION 

CUMMUN /L/ XARCAL se XALMGe VAREAL +e YALMGeXOIMoeYbIMG 

COMPLEX XKI,XIdUD0Y~+CHAR 

REAL MACHN, KAPPA 

CUMMUN X1( 200) pUREAL(20U) pVREAL( 200) eAREAL (200) te YREAL(20U) pePsl ( 
20U) eX1LBOOY (200) eUBS0DY (200) sVS0VUY(20U) +X300Y(2U0) sYBUDT(20U)»s 


2 CHAR (41910) eLLINE( 200) +S( 200) sYACHN( 2U0) eX APPA(20U) ,ANGL (200) 


THe TA( 200) eH 200) 
COMMON /U/ PL eAhseARPeNRNeNXeNCReL OMe IPLT MIDE NG eNCeNiie Ll eddy KK 
COMMUN/SG/S NleN2eNGeN7eN5e VG, 41 
COMPLEX ROTATE 
JLIMENSTON MPOINT(41)e¢PSIIMP(41) 
LOGILAL 1Sh 
UATA 1SW /.FALSEe/ « NAR 4417 ¢ MAR /10/ 
LF. ¢LSW) GO TO § 
4SW = TRUE. 
ANGROT = ATAN2(-X31MG,eYB1MG) 
RUTATE = CMPLXICOUS(ANGROI) «SINC ANOROT)) 
es. 0 
CHECK FOR FIRSI ROW OF TRLANGLE 
AF ¢K1.Gl.1) GU 10 30 
FAC = =PS1(1)/485(PS1(1)) 


“KQ = MRP 

VU 20 J = 1eK2eMRQ 

LE = t+¢{i 

LCHAR(Le1) = RUITATE*CMPLX(AXREAL( JU) se YREAL (JU) ) 
PST IMPIL) = PSi¢tJ) 

"POINI(L) = 2 


LoS toa 

CHECK TO SEE IF WE HAVE CXCEEDED THE JIMENSION OF CHAK 
dr (L.LEeNAR) RETJRN 

WAITE (N4.100) 


MRQ = MKU+4RV 

& = UV 

oY TU 10 

JU 70 J = 1+K2e4%Ra 
. = &+j 


Ir (FAC#® PSI(J),Le.0.) GU TO 60 
PUINI 1S UUTSIUe THE ROUT 
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NANNY 


MPOINTC(L) = MPOJNT(L) +1 

LF (MPOINT(L)-etQe2) GO TU 50 
LHECK FOR CIMENSLON UF CHAR 

LF (MPOINT(L)eLE»YAR) GO FO 60 
MPOINT(L) = 2 

LALL CPLUT(CHAR(Lel1) +3) 

UO 40 K = 2,MAR 

CALL CPLUT (LHAR(L+K) 22) 
CHAR(Le1) = CHAR(LeMAR) 

SO TU 60 

FIND A POINT ON THE BODY 3Y LINCAR INTERPIULATION 


K2-= PSIUd)=-PSITMP(L) 
Rl = PSI(J)/R2 
K2 = =-PSLIMP(LI/R2 


CHAR(Le1) = RL*ECHAR(Le1) FROTATE*®R2*CMPLX(AREAL(U) Pe YREAL (IJ) ) 
STORE THe POINT FOR PLOTTING 

A = MPOINI(L) 

CHAR(IL eK) = RUTATE*CMPLX(XREAL (J) e YREAL (UJ) ) 
PSTINP (CL) s= ePsiLod) 

AF (K2.GTe1) RETURN 

We HAVE FINISHED THIS SUPCRSUNIC PATH 
FINISH PLOTTING THe CHARACTERISTICS 

VO 9U J = jJelt 

E> MPOLING Cup 

De Cotes) GO sTO .90 

CALL CPLOT(CHAR(Jel1) ¢3) 

WO 80 K = 2,L 

LALL CPLOT(CHAR(JUeK) +2) 

CONTINUE 

RETURN 

FURMAT(2@/7H TOO MANY SUPERSONIC POINTS ) 
tNo 


SUBRUUTINE LNICRP(K1+K2eK3) 

CUMPUTES POINIS UN THE B8UU0Y BY AUAURATIC LNIERPOLATION 

'dE BOOY IS THe STREAMLINete PSI = Oe 

K1,K2 ANU KZ ARE 3 CONSECUTIVE POINTS 

KS IS ZERO wHeN THRE IS NO BODY PULNE LN THAIS INTERVAL 

COMPLEX XI,XLBUDYe CHAR 

REAL MACHN,KAPPA 

LCUMYON XL( 200) eUREAL(20U) eVREAL( 200) sXREAL( 200) *YREAL( 200) gPS1( 
1 20U) -X1LBCOY( 200) sUBODDY( 200) eVBUVUY(20U) -X30UY(2UU) eYBODY(20U) + 
2 CHAK( 41010) sLLINE (2060) 06S 200) eMACHN(2U0) eK APPA( 200) ,ANGL( 200) 
3 eTHETA( 200) eH( 200) 

COMPLEX X113+eX123 


site jee et 8 
X85 = PSI(K1741) 
AC = PSI(K1+2) 


CHECK FOR A SlvoN CHANGE BeTwEEN THE SttONU AND TALRD INTcRVAL 
Lr (Xs*XC.LT°¢O02) GO TO 20 

IF (X8*XC.GTe¢Ue) RETURN 

CHECK TO Stet LF BOTH XB AND XC ARe -ZERVU 
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Lr (XB+XC.EGe00) RETURN 

A3 = 1 

XA = (PSEK1) 

EA = -xXC=XB 

Es ta. XC=KA 

EG: = X6=XA 

YU LINEAR INTERPOLATION BETWEEN PUINTS 2 AND 3 
U23 = (UREAL(K141)*XC=URLAL(K14+2)*XB)/04N 
V253 = (VREAL(AKA141)*XC-VRCAL(K1+2) *XB)/tA 
K25 = (XREALIK1L4+1)¥XC-XREAL(K1+4+2)*XBI/C0A 
Y23 = (YREAL(K141)*XCHYRCAL(K142)*XB)/CA 


Al25 = (XI (K1+1)"XCH-xXI(K1+2)*XGb)/tA 
UU LINEAK INTtRPOLATION scTweteNn PULNTS 1 AND 3 


Y1Z = (UREAL(OK1) *XCHURLEAL(K14+2)*XA) 763 
V13 = (VREAL(K1) *XCHVREAL(K142) *XA)/03 
X14 = (XREAL(IAL) ¥*¥CHXREAL(K142)*XA)/03 
Y1Z = (YREAL(IK1) *XCHYREAL(K14+2)*xXA)/63 


K115 = (AL (K1)*XCHXT(K1+2) XA) /e33 
1F PST IS NOT MUNOTONIC USE LINEAR INILRPJULATIUN 
AFL (CEAXEC*GT¢0%) GO TO 30 


WS 4 wi: 

UBOUY(K2) = U2s 
VBOUY(K2) = V23 
XBOUDY(K2) = X25 
YTSOVUY(K2) = Y25 
XLBOUVY(K2) = X123 
Ke TUK 


JO QUADRATIC LNTeRPOLATIUN 


USODY(K2) = (UL3*XB=-U23*XA) SEC 
VSOUT(K2) = (V13*X8=V23"XA)seC 
RBOIV (Ke) = (X15*#XB-KX2354#KAN) seC 
YHOUT(K2) = (Y13*Xdb-Y23#KA) SEC 
KLBOIUY(K2) = (X1L13%Xb-X125*xA) EC 
KE TUKN 

tCNO 


SUBRUUTINE BLAUC(A1) 

KeORVER IHE K1 BODY POINISsPRINI AND PLOT FHte ReSULIS 

Xt AL MACH 

CUMMON /A/ GAMMAs4ACH 

COMMUN /C/ XAREAL e XAIMGs YAREAL se YALMGeKUIMGe YUI MG 

CUMMUN /F/ KORLG,YURIG,SLALE 

CUMPLEX AI,XI3U0Y«CHAR 

McAL MACHN, KAPPA 

LCUMMUN X1( 200) sUREAL(2Z20U) eVREAL(20U) eXREAL (200) se YREAL( 2UU) PSI, 
1 2U0UU) eXLS8COY( 200) sUBOULY (200) eVG0UY(20U) sASO0UY(2U0) +e Y¥nUUT(c0U)» 
2 CHAR (41610) eLLINE( 200) 451200) sMACHN( 2U0) eX OPPA(CUU) ANGEL ( eCuU) 
3 eTHETA(200)+4(200) 

CUMYUN /U/ PLeANeAKP eNRNO NX eNCReLoMe LPL IE eMUDE se NG ae NC eNive ll euie AX 
LCOMMUN/SG/S AleN2eN6eN7eN3Se V4, 41 

CUMPLEX AITAIL + BPePOLE eALAe KIBVAIC 

CUMMON /HZ BPILO) ePOLE eXAAeXLBe ALCO XLPALL es PLTSLdendrs 


Real MaCHaAeYACHBAYATALL 

CUMMUON SKS RNOePCHU + PCHL MACHA sANGLAMIPAIL sAQGLI eMACH3ZyANGLSy RL 
1 ANGLNe Ik 

COMMUN S1L1/ SFeSIZEseANG ae XMAX se YMAXeXURGVURePSSIZ 
COMPLEX LINT(5) sBASE(3) 

CUMMON ZIT, VITLECE) 

CMACH(QS) = SURT(IASS(COSU=-((GAMMAZ1.)/2e) * aS) ) 
GAM = (GAMMN-1.)/26 

CUSQ = GAM + Lle/(MACH*MACH) 

KAD = 180./7P1L 

1PL = MOD(CIPLT/10.10) 

onl =: SiaASk 

ANGROT = ATAN2(-X31MGy-YKIMG) 

CUSA = CUS(ANGROI) 

SINA = SILN(ANGROT) 

SUL = 2e¥*¥PI*SUNT (XB IMG AXGLMG+YBIMG*YBIMG) 

DPE CRS Les i); GO’ 12) 25 

CREATE A POINTER LIST 

VV 10 J = 14h1 

AL INE(J) = J 

URDER THE POINTS AROUND THE BODY AND FLND THE IALCKNESS/CHORY 
TEu= 30 

LMIAXesS 7 

REBODY CL) = XD VALE 

VV 2U J = 2eK1 

K = ILINte (Jel) 

LMAXP = LMAX 

KMX = XBUDY(K) 

LMAX = J=1 

UO $U L = JeKl 

Ce = re TNE ( Ee) 

AF (XBODY(LL)ILE.XMX) GO TO 50 

XX = XBUDY(LL) 

LMAX = L 

CONTINUE 

IO = AMAX1(TCrHABS(YBUDY(LYAK)-YB00Y(LMAXP))) 
LF (LMAX-EQG,(J=-1)) GO TO 2U 

LOEPNECd<=1) = LEINECEMAX) 

LLINE(LMAX) = K 

CONTINUE 

SEZeE = «07 

SUL = SOL¥*SCALE 

WRITE (N4e11U0) TCeSOL 

[THA = 90--ANGLA = KAU#ANGXOT 

'HB = 9Q."-ANGLB = RAU*XANGXOUT 

TURN = IHA-THB 

WRITVO (N4e120) MACHB se TH3e MACHAs THAe TURN 
LEeCKL SCE ST) VRETURN 

LL = ILINE(K1) 

BASE(1) = xIB 


LINT(1) = CONJG(XIBODY(LL)-x1B) 
SASE(2) = xIA 

CLINI (2) = CONJo(XIB-XIA) 
BASE(3) = xITALL = 1.E=6 

CINE (3) = CONJO(XIN-XITALL) 
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RAGT-= Gls 


ANG = 45-6 
re Le: 
tl = -1 


SO) VOr PU) we, = tei 
ee BETINE: CLs) 
Cige = 
tr (REAL (X1B0UY(J)).G6T, REAL (XI) ) L2 =,2 
PFA CREAACKIBODY (J) ) GTOREAL (XDA). £2. =. 3 
LF (FACT¥AIMAG((XLBODY (JU) -B3ASE(L2))*ETVI(L2))-LIeUe) SO TU 40 
BaP Sets 11 
KREAL(L1) = XSUDY(JU) *COSA-YBUDY(U)*SINA 
YREAL(L1) = AXSUDY(J)4*SINA+Y S001 (4U)*CUSA 
WS = JBOUY (J) *UB00Y (VU) +VbUDY (JU) #VBODY (J) 
tMN = EMACH (QS) 
MAGHNGE1L) = JAIN 
ANGLIL1) = ATAN2Z(JBODY (J) *SINA+VB0UUY (J) *CUSAeUBUUY (J) *CUSA= 
1 VBOUY(J)*SINA) 
wee GLP sits) 40:): 16,0; of.0¢ 0 
CALL CSYMBL (C4PLX(XREAL(L1) «SM*(-5SO0FEMIN) ) 0359-1) 
4Q CUNTINUE 
EPRAUF AGT 6 BG) .G0) i710) 160 


ANG = 0-6 
Eee: (0 

Pe LS ae 
FACT = -1. 
9V TU 50 


60 “P = K1+1 
PLOT MACH NUMBER AT THE FAIL POINI 


ANG = 45e 

tMN = SM*(MACHN(1)+.50) 

CALL CSYMBLICMPLX(XREAL( 1) ¢EMN) 030-1) 
ANG = 0. 

KREAL(KP) = XREAL(1) 

TREAL(KP) = YREAL(1)+TR*SUL 

ANGL(KP) = ANGL(1) 


“MACHN(KP) = MACHN(1) 

[Pos GOLRGALUL ) als he Oia) OGAL LS CPLOT « CCAPLX (Us y =YREAL CL Pe ad} 
CALL OUTRPI (KP) 

AF (LPL eEG.0) RETURN 

CALL ORIGIN (5051-5) 


ANG = 0. 
Si2te sue ko 
Se = il], 


CALL CSYMBL ((*,.80+6820) e1lHMel) 

cMA = ABSIMACHA) 

tMB = ABS(4MACHY) 

cNCUUE (60,14UeTILILE) E%S et MA eT URNe SOL 

CALL CSYMBL((Uee=eS)eTIFLe 60) 

CALL KYAKLS (CH LeOe Pe Dd) eVUe eGo Ue ecUs GU.) 

CALL FRAMER 

METURN 
110 FURMAT (CIAL SSSA LBATHICKNESS/SLCHURUAeF6e4//3IKeLONVPAP/SCHORVEa Foe 5/) 
120 FURMAT (15x slLYHLVNLET MACH NUMBER=eFO,S5e1DXeL7HINLET FLOW ANGLER» 

1 F7e2//15X— ISH CXIT MACH NUMBERS eh 6e.35915Xe L7H CALT FLOW ANOLEX-e 
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140 


20 


10 


2 FT ee//3OXeLHAHTURNING ANOLE=eFTec//) 
FORMAT (4H MIEF4 35 ,5A3H“M2=F4.3,95X7THOEL TH=F 6.2, 9A4HG/C=F 4,2) 
cNO 


SUBRUUTINE BAUPTS(K1) 

cCLIMINATE MOST ERONEOUS POINIS 

COMPLEX X1I,xLlBUIY,CHAR 

MtAL MACHN, KAPPA 

COMMON XI(2UU) sUREAL(20U) eVREAL( 200) eXREAL (200) *#YREAL(2UU) PSI ( 
1 200) eXIBCDY( 200) eUBODY (200) sVBUUY(20U).x30UY(2U0) sYBUUDY(20U). 
2 CHAR( 4110) eLLINE( 200) 45(200) »MACHN(2U0) »XAPPA(2CU0) ,ANGL(20u) 
3 eTHETA(200) eH( 200) 

CUMYMUN SCS XARLAL eXAIMG,YAREAL se YAIMG se xXbIMGe YBIMG 

COMMON SUS PLeAAe4IRPeNKNe NX eNCReLoMeIPLT se “WIDE eNO tNCeVINe tLe Ute KK 
XT = XRBOVY(1)+¢0YBOUY(K1)°YB0U0Y(1))*(XB1L4MG/Y3IMG) 

Lr (UBODY(K1) eLTemo2) RETURN 

LUOK FOR POINIS ON THE CENTRAL STREAMLINE 

LF (XBODY(K1)eLT-0-) GO FY 10 

LF (XBODTY(K1) ¢STeXIT) RETURN 

SAVE THIS POINI 

Al = MINU(K141+¢LB4) 

RETURN 

LF (ABS (UBODY(KA1)) LE«ABS(VBO0UY(K1))) GU TO 20 

RETURN 

tNO 


SUBRYVUTINE HUORF(KA1) 

COMPLEX KO] 

CUMPLEX XL,~XLBUDYeCHAR 

REAL MACHN, KAPPA 

COMMON KL(20U) sUREAL(20U) sVREAL( CUD) eAREAL (200) eYREAL(20U) PSI ( 
1 200)4xX1B80DY(200) sUBUDY (200) eV3UDY(20U) eX300Y(2UU) eYEVUDY (200). 
2 CHAR( 41910) eLLINE(C2U0) eS(2CU) eMACHN(2uU0) eX APPA(CUU) -ANGL(20U) 
3 eTHETA( 200) eH 200) 

COMMON /U/ PLleANeARPeNRNONXeNCReLBMeIPLTs MODE eNSeNCoNNeLLe Je KK 
COMPLEX KXITAIL ese 3PePOLE se XLAeXl be ALC 

COMMUN /H/ BP(1O) ePOLE se XLhe X18 eXALCeXLTAIL ye PLTSLdeNBPS 

COMMON /I/ SFaSILZEs ANG ae XMAX se YMAXeXURC VYURePSSIZ 

EP = MOD. CEPLT.10) 

tr (LPL.LQ.G) RETURN 

XI1(NC) = XIC 

X LUNGS) = KIC 

CALL LOCUS (IPL) 

LE CKD SG. 1) RETURN 

her U 

DUS TO J =) Pek a 

CHECK FOR SUPERSONIC POINTS 

Ir (ABSVAIMAG(XIBOVY(J)))-GT.10-) GO FY 10 

A = K+l 
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10 


15 


20 


30 


40 


50 


CALL CSYMBL(ALSOUY( JU) 539-1) 
CONTINUE 

PPR =4IPET./ 100 

LF (2PL.EQ.0) RETUKN 

LE (LPE«EG. 2) GO. TiO SO 

PAL = SPL a7 

Lr (KeGE.~K1,) GO TO 15 

PLOT IHE SUPERSONIC I[RIANGLE 

CALL FRAMEK 

CALL URIVIN (7230725) 

SF) =) 10% 

VU 20 J = 1eKl 

LF (ABSCVAIMAGCXIBODY(J)))-+LT¢10~.) GO TY eu 
TEMP = AIMAG(XL30IY(U)) 

(EMP = TEMP] SIGN(100eeTEMP) 
XLBOVY(J) = CMPLX(REAL(XLSO0Y(U)) oe TEMP) 
CALL’ CSYMBL (XIS00Y (J) eSe-1) 
CONTINUE 

CALL KYAXKES (UO eeOe DV eVoeedetelteOe) 
GAL UY XKYVAXKES(( U5 0's) tee SeveleAOe) 
CALL CPLUT (Cee 7eme7) 93) 

CALE CPLOUT (063253) 92) 

LF (LPL.EQ,1) RETURN 

CALL FRAMEK 

CALL ORIGIN (220Ue5e0) 

XMAX = 6U, 

Seas’ SOs 

KP = Kitt 

RUT = CMPLX(XREAL( KP) e=YRLALC KP) ) 
KUT = ROI/CABS(ROT) 

1PEN = 3 

DOr ie =p, he 

CALL CPLUT (RUF eCYPLX(XREAL(L) se YREAL(L) ) oe LPEN) 
IPEV = 2 

JPEN = 3 

VO 50 L = 14KP 

VELS = THETA &aHCL) 

XS = XREAL(L) + VELS#*SIN( ANGLE (L)) 
YS = YREAL(L) + DELS*COS(ANGL(L)) 
CALL CPLUT (ROT#CYMPLX(XS*YS)eIPeEN) 
IPEN = 2 

RETURN 

tND 


SUBRUUTINE LOCUS(IPL) 

COMPLEX SeXTOFSe XT eXlLOeCSART eo XIMINe XILMAXKeASOAPRec IY 

ReAL MACHs+LOGS 

LUMMON /A/ GAMMA, MACH 

COMMUN /U/ PL eAe ARP eNRNe NX eNCReL OMe IPLET,VODE Noe NCaeNN eli edie MAK 
LOMAON /t/ LUGS (412) EI UIG4) XR eNLUGS eA NK e NP eNJe iv (32) 

CUMPLEX XITAIL +s SP ePOLE es XLAVALB AIC 

LCUMMON /H/ BP(10) POLE XtAe KL ALCeXLPATL e PL ISL de NorPs 
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FINU POINT ON BOUNDRY ANU DRAW LINE TO IT 
TP ENE t= <2 
LE (CLPENSEG, 3S) GO 10. 20 
XP = XMAX 
BE CXR se Lb es) AP -= O6 
YO YP = YYF¢CYLA*YV)/S(XL“XX))*(XPHXKA) 
dr (ABSCYP+YP=YMAX).GT.YMAX) GO TO 110 
§50 CALL PLOI (xXPeYPeIABS(IPENL) ) 
LF (LPENLeEGe3) GO TO 10 
bU 1U 20 
60 IPENL = LABS(LPENL) 
JF (LPENLeEQe2) GO TO 20 
FIND POLNT ON SUUNURY ANU ORAW LINt TO IT 
LPENL = 2 
1F (LPENeEQ.5) GO TO 20 
YP >= YMAX 
REMI ete.) VP = Oe 
70 XP = XXF((XLOEXAVS(YLAVY) IV *(YPHYY) 
60 TU 50 
80 IF (LABS(IPEN).tu.3) GO 10 10 
IPENE. = TPENE+1L 
1F (LPENL.EGe5) GO TO 70 
LPENE -= 5 
GU 10 4O 
90 CALL FRAMEK 
XL = 0. 
YL = 0. 
LPENE = 0 
KETUKN 
100 XOR = XP 
YOR = YP 
RETURN 
GO 10 60 
cNO 


SUBROUTINE CSYM3L (X+NeL) 

CUMMUN SL/S Ske SIZE eANG A XMAXeYMAXe XURVVURe PSSIZ 
UVIMENSION x(2) 

CHANGE RELATIVE MOVEMENTS TO ABSULUTE JNCHES 
XX = KOR+SF*¥X(1) 

VY J=) VORtSFE#X (2) 

CHECK TO SEE IF we ARE WITHIN THe PAGE 

LF (CXXeL Fee) eORe (VV eL T0004) eORe (XXe Gl e XMAX) -ORe (VV e Gl eo YMAX) ) 
TRETURN 

LALL SYMBOL (XXeYYe+SLZEyNeANGeL) 

RETURN 

tO 
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ANNAN 


2) 


30 


29 


30 


10 


20 


30 


PI(2) = YMIN+FLOAT(NLIN]=J) €}y1 
CALE PNAC PT HP) 

CUNTLNUE 

Re TUR 

ND 


SUBROUTINE GOUPLUT (NR) 
LNITLATE PLOT 
RROD RD ROR RK ERD RE EKA KR EERE EEE K ED EERE ERE K EERE RE ERE KEE 
IHIS SUBROUTINte SHUULD BE REPLACED BY ANY ROUTING WHICH LvVSIRUCIS 
PHE “SYSTEM: TO LNLTLATEs A PLOT 
OR KOR ORR RE ROR OR RR RR RRR RR ERA RE REE REE EKER KEKE ERK EEE 
UDIMceNSION 103) 
14U(3) = 0 
NRUN = IABS(NRN) 
tCNCUVE (20,3Us1)) NRUN 
LF (NRNeLT.O) GO TO 25 
GAEE ‘PEGS (120 -LDe«LLeylle) 
RETIN 
CALL PLOTSB6BL (120+¢10e11-ee11,.) 
RETIRN 
FURMAIL ( Sxe3HRUN «© 15) 
cNo 


SUBRUUTING CPLUP (XeIPEN) 

MUVE PEN TO THe POSITION VEFIweD BY X WHEN X IS Jw FHt PAor 
COMMUN JZLZY Sk eSIZEe ANG se KMAXe YMAXe KURA VYURePSS1LZ 
UVIMENSION x(2) 

VATA XLoeVYL,IPeENe S0e+0e9U/ 

CHANGE RtELATIVe MOVEMENTS TO ABSOLUTE LNCHes 
KX = XOR+SF4X(1) 

YY = YOR+SFeX(2) 

CHECK TU SEE LF we ARE wiIHIN THe PAGE 

Lr (ABSUXX4+XX=-X4AK) -GTeXMAX) GO JU 30 

tr (ABSITY4YY=YMAX).GT.Y“MAX) GU TU 60 

dr (LPENLE.NEeO) GO TO 80 

CALL PLOI (xXXeVY¥eIABS(IPLN) ) 


APENL = 0 
Ak = AR 
VOI S) TY 


AF (L4PENSGT.9) RETURN 
dr (LPENL.NE*O) GO TU 100 


XOK = XX 
TUR = YY 
RETURN 


IPENL =-<LARSC(APENL) 

Le (UXXe OT XMAK) CAND. (IPENLT2U0)) GO FO 90 
CHECK TO SEE IF PEN IS ALREQUY OUT UF BOUWS 
LF (APENLeNE*O) GO TO 20 
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op Ja corte op I ge) 


10 


10 


20 


CUMPUTE THe LCAARACIER NUMS3ER 
lWO = (P(1)=-X“M LN) «XFACH1~-9 
COMPUITtC IHE LINE NUMBER —- OWE 
JLN = FLUAT(IUNLINW@IFIX( (P12) °YMINI *¥YFAL) II &°.5 
CHECK TO SEE LF we ARE WLIHIN PRESCRIBED LIMITS 
Je (CIWUeLE 0) OR (ULN LI 0) eURe (lw GIeLSZ) OR (JLN,GI.NLIN) ) 
1URN 

CUMPUTE THE WORD NUMBER 
J = 1+(IWU+JULCN*LSZ)/NCH 
COMPUTE THE CHARACTER NUMBER 
ACH = 14+™U0U(LWUeNCH) 
REPLACE IHE CHARATIER BY L 
**x ENON ANS] = SEE WRITCUP AT LCNU*¥4* 
VeCWUe (NCH e10ePG(J)) (LYK) eK = leNCH) 
LUJ(KCH) = L 

**€kANONEANST = SEE WRITEUP Al ECNOK*¥** 
CNCUVDeE (NCHe1l10ePG(I)) (LYK) eK = le NCH) 
RETURN 

FURMAT (10401) 

tNpD 


SUBRUUTINE PNAKILS (P) 

URKAW AXES THRUUGH THE POLVT P 

THE X-AXLS wltl CONSIST UF eet - WHERE THe + APPEARS EVteRY i0U 
CHARACTERS ON tHE LINE 

IHE YeAXLS CONSISTS OF + EVERY FLiVeE LINES FROM THE POINT P 
ANTebeR HP,H" 

UIMENSION P(2)0 PT(2) 

COMMON SUS POI4H0S) eXMINOYMINGOXFACOYFAC NCAP LSZeNELNe TwWOe UL Ne KCH 
>EE SUBRUUTINE PN FOR THte MEANING OF THe cleeMENIS LN THIS sLUCK 
StT THE GRAPH TO ALL BLANKS 

VATA HBeHP,yHM/1LH glHteln-/ 

11 = NLIN*LSZ/NCH 

UO 20 Je= 1 eit 

P6O(J) = HB 

AL = 1. /XF AE 

CALL PN (P,HP) 

Lt =" 1Wb 

BPibGeie=s (Ce) 

UU X-AXIS 

ES) = £S2=1 

UO 20° J =12LtS 

PT(1) = XMINt+FLOAT(J=1) «HL 

L = HY 

dr (MOO(J-ITI+1U).EU.U) L = HP 

CALL PN (PT,L) 

VU Y-AXIS 

Hl = 1./7YFAC 

11 = L4+NLIN-JLYN 

Phen = 707) 

VO 3U J = 1eNLIN 

Ir (MOO(J=ITI1+5).NE-G) GU 10 450 


84 


RET 


30 


110 
146 
130 


1 


IHeETAH= RTHO/RIH 

IHT = THETAH 

IHETA(L) = THeTAH 

FC = 1,.0+.066*“HSQ-,U08#MNeMHSQ 

FR = 1,-6134*4HSU+.027% 4HSQ «MH 

VO AT MOST 200 ITERATIONS 

DOs 29 00 et 1200 

RIAU= LeJS(FCK(2e4711*ALOG(FRERTH* THETAN) +4475) 4+10¢34GF4+1/24./ 
(Gt+*GE+2(10+)-16.87) 

1AU = RTAUsRTAU 

HB = le/l1l.-GE¥RTAU) 

4H = (Hd4+10)*(1.+.178*4HS0) 1. 

deP = -TALTAN*UQUS 

PLE= HH*ostEPsIAJ 

PIE= AMAX1(PIMLNeAMINI(PLMAXsPIE) ) 

6b = 6.1*d0RT(P1E4+1.81)-1e7 


12 = ABS(G-GE)/Ge 

Se S'S 

Ulo = OT 

wit = (HH +267 4HSQ)*SEP+IAU 


Ee COED. 1). GO fOr 110 

§1 = ABS((OT*VUI2)/vUT) 

Cee Vbieel, TeS2)cANDe(CT2ek1.7TEL)) GOTO >k50 
IHEITAH= THT+e5*0T*US 

CONTINUE 

IHETAC(LP) = UVI*FQStTHT 

CHE-PAWS = THETA EP) 

AT = Fre lA Cee) 

4(L) = (NL) *US+H440SO0OLU)7(0S+USOLU) 

Ha(LeP) = HA 
JSOLU = US 

Ca) ER 

dr (lLeNEeK2) GY TO 10 


H(K2) = 2e#H(K2)-H(KE-INC) 
BCKi2 SAGs. 

XE TURN 

tNO 


SUBRUUTINE PN (Pel) 

PUT CHARACTER L I1VIO ROW AND CULUMN DEIERMINED BY PUSIILUN P 
JIMENSION P(2)0 145410) 

CUMMUN /IJ/ PG(4035) eXMINe YMINeXFAC Se YFACeNCAeLSZeNLANe TWOeUL Ne ACA 
AMINe YMIN IS THe LOWER LEFT HAND CORNER OF THE GXAPH 

AFACsYFAC FRE Ide SCALING INCREMtNIS FUR X ANDO Y 

NCH 14S IME NUMBER OF CHARACTERS LN A MACHINE WORY.10 FOR CIC 6590 
LSZ 1S IHE NUMBER OF CHARACTERS WHICH wILL COMPRAdt A PRINI LINE 
NLIN IS THE NUMdgER OF LINcS [THe GRAPH WILL AAVE 

IHE CHARACTER & wILL APPEAR IN [He IWvU AWURD UN LINE JLN 

ACH 1S THE NUMBER OF THE CHARACICR WITHIN THe WORU WHICH wlLe Be 
MREPLACEOD sy IHe lot CHARACTER CONTAINeV IV L 

IHE VIMENSION UF PG MUST 3E NLIN@eSZ/NCH 

StE WLOCK QATA FOR THE DEFINITION UF CLEMENTS IN FTHLS BL UCM 


SeP = =$THETIA(L) e0QUS 
WRITE (N%4+710) KREAL(L) se YREAL(L) eANGe KAPPA(L) eMALAN(L) 
LL STHETACL) «SEP axXSsYS 
120 (EGC) = kGe2 
ar (CC FE. 55). G0, 10° 130 
Lc = | 
WRITE (N4.200) 
WRITE (N4e700) 
130 CONTINUE 
Re TURN 
200 FURMAI (1H1//) 
690 FORMAT(13X,4¥9HCOORDINATES FROM LOWER SURFACE TALL TO UPPCR 
Jel2qSURFACE TALL /) 
700 FURMAL € DXelHAeBXe LH +e bXe SHANG eOXeLHK se 7TXReLHMeOXsDHTHL Tae dA, 
1 SHSEP se 7Xe@HXSe7Xe2HYS/S/) 
710 FURMAT (F11,SeF9e5eFS 1 eFleoeF Be 4 e4F ILD) 
720 FORMAT (FID Der Dede FFel eh Teer F Bo tet Xe ALOe 4+Xe2F9-9) 
tND 


SUBROUTINE NASHYMC (K1leKe2) 

CUMPUTE THE BUUNDRY LAYER FRUM PULNT Kl TJ) K2 

CUMMUN /A/S GAMMA, 4ACH 

CUMPLEX X1L,XIBUDY,CHAR 

XtAL MACHN, KAPPA 

LUMMON X1( 200) sUREAL( 200) eVREAL( 200) eXREAL( 200) *YREAL (200) ePSI ( 
1 2UU) eXIBCOY (200) sUds0OY( 200) eVACUY(20U) eX30UY(2UU) sYBOUY(20U) 
2 CHAK(41 910) eLLIWE (200) e5( 200) e4ACHN(2U0) eX APPA( 200) ,ANGL( 200) 
3 *THEIA( 2006) +H( 200) 

REAL MACHA,MACH YS eMITAIL 

LOMMON SKS RNePCHIePCHL se MACHA ANGLAMTALL +e ANGLT s4ACHQeANOLSyRL oe 
1 ANGLUNe IK 

REAL YMHeMHSGeNUe MACH 

VATA IR 4e3424/ 4 RTHOsS320.7 « IEL/5.E°3/ »« TE2/79eEH=5/ « PIMIN 
L 42-165/ +6 PIMAK/1.b4/ 

GAM = (GAMMA=-1.)/2-6 

GAML = 1./(GAM“A-1.) 

CSIINF = 1,+GA4%«MACH*MACH 

INC = ISIGN(1lsKe-K1) 

bE = 6,5 

ime 3 9 Sh f 

VSOLU = ABS(S(L)=-S(L=INC) ) 

10 LP = L+INC 

“H = «5S*(4ACHN(L) +MACHN(LP) ) 

MHSQ = MAHeMH 

CSTH = Let AM*MHSYQ 

US = ABSIS(LP)=S(L)) 


VIDS = (MACHN(LP) =YMACHN(L))/(0S*MH*CSIN) 
f = CSIINF/CSIA 

KHOH = T**«GAM1 

NU = [Te(1.4+TRI/S(RAGH*A( THEIR) ) 

KtH = RN*MH/(CMACH*ENU) 


fF CE SNE SKI) GO’ TO 3g 


70 


80 


90 


100 


105 


IHETACL) = Oc 

LF (MACHN(L) eGT.Q4MAX) KGMAX 
dr (MACHN(L) LI,GYIN) KoMLN 
GMIN = AMINICMACHN(L) sQYIN) 
IMAX = AMAXL(IMACHN(L) «QYMAXK) 
UX = XREAL(LI“AREAL(LA-1) 

OY = YREAL(L)=“YREAL(t=1) 

LF (ANGL(L)-ANGL(L-1)eGTePI72.) ANGL(L) = ANGL(LI=Pl 
JIH = ANGL(L)-ANOL(L=1) 

PAC = TT. 

Lr (UTHeNE,Ue) FAC = oS*UIH/SIN(29D*UTH) 

US = FAC*eSOURT(UX*OX4UY*«0Y) 

SCL) = S(L=7)!+US 

KAPPA(L=1) ="U¥IA/0S 


“ou 
ro 


CONTINUE 
LRU = KP+2 
ty) f° | Leda 


DE CRNSEE.O,) GO FO LOS 

L = KUMAK 

Ee CRGHU elt. Ue) 8G "ro s0 

VO 79 L = KQMIN,KP 

fr (XREAE(L)  -GtresPCHU) GO 10 80 

CONTINUE 

CALL ABORT 

KXIRANS = XREALOL) 

LFRY = L 

CALL NASHMC (LeKP) 

Te (PCHLeGT,0.) XTRANS = PCHL 

UU 90 L = 1eKOMIN 

LF (XREAL(L)eLESXTRANS) GU TO 100 

CONTINUE 

LF (KN,NE.G.) CALL ABORI 

SrRO =": 

CALL NASHMC(LIRLel1) 

WRITE (N4e690) 

WRITt (N4.7CU) 

CC =°15 

eNG o> 4 

V0 150 L = 1¢KP 

ANG = RAUV#ANGL(L) 

4F (LPL.6T.1) CALL CSYMBLICMPLX(XREAL(L) ce YREAL(L)I ) 015-1) 
VELS = THETA(L) eH (L) 

RS = KXREAL(L)F+UCLS*SIN(ANoL(L)) 

TS = YREAL(LIFVELS#*CUS(ANOL(L)) 

CF CCE SE PeL IRL PeIRel(EsGlekE TRU) ) 60 FO 110 

“SG = 1H 

LF (MACHN(L).LI.TOL) “SG = LOHSTAGNATION 

LF (CL eEMe LITRE #1) eORe (LetO LIRU -1)) M3sS = LOAIRANSITLUN 
ARTTO (N44720) XREAL(L) eYREAL(L) PANG KAPPAIL) pMALHN(L) eS0eKOeTS 
INC = -] 

bY TU 12U 

CUNT IANUE 

CrP = ERNE 

LSI = Le tGAMe4ACHN(L) &4ACHNIL) 

UVUOS = (MACHN(L) -4ACHN(LP)) SI CST e®MACHNIL) ©ABS(S(LIHS(LP))) 


10 


20 


SUBRUUTINE LINeS (K1sK2) 
LOMMON ARRAY(1) 

COMPLEX ARKAY 

CALL CPLOT (ARRAY(K1) 43) 
Wess 1K] 

JO 10 J = KPeKe2 

CALL CPLUT (ARRAY (JU) 9?) 
Ke TURN 

cCND 


SUBROUTINE POLNIS (K1eK2) 
LOMMON ARRAY(1) 

CUMPLLEX ARRAY 

OO. 20) J = Ki dK] 

CALL CSYMBL (AKRAY(JU) 0260-1) 
RETURN 

CNO 


SUBRUUTINE GUIPT (KP) 

ReAL MACH : 

LOMYON /A/ GAMMA, 4ACH 

COMPLEX XI,X1LBUUY»+CHAR 

REAL MACON, KAPPA 

COMMON XI (200) eUREAL (200) eVREAL(2UC) sAREAL (200) eYREAL(2U0) »PSI( 
1 2UU) -XLBU0Y (200) sUSUDY (200) eVBUVY(20U) eX 30UY(2U0) eYBOUDY(E0U)s 
2 CHAR(41¢106) eLL INE (200) +5¢4200) «eMACHN(2U0) eXKAPPA(2U0) ,ANGL(20U) 
32 ee THETA 200) eH( 200) 

COMMON SUS PI eANeYMRPeNRNeNxXeNCReLBMeLPLT, VIDE yNBeNCeoNive ile Ue KK 
CUMMON/SGS NI*N2eN5eN7eNZeNGe M1 

COMPLEX KITAILeBPePOLE es XLAeXIBAAKLC 

COMMON /H/Z BP( 10) ePOLE se XLbeXLBe ALC eXLTALL ePLISI Ze NBPS 

COMMON SL/S SFaedSIZEeANG eo XMAXe YMAXeXURA VURPSSIZ 

Re AL MACHA.MACHB«¢MTAIL 

LCUMMUN /K/ RNePCHUePCHL eMACHAsANGLA,MPAILeANGLI eACHRy~ANOLS RKC, 
1 ANGLNeIR 


LPL = MOU(IPLT/10.10) 
KAD = 180./Pl 

10L = 1.£#5 

GAM = ,5*(GAMMA-1,.) 
KYMIN = 1 

AQMAX = 1 

Q™MIN = MACHN(1) 
W4AX = QMIN 

H¢(¢7) = Oe 

(IHETA(1) = O- 

o(1)) =. Oe 

KAPPA(KP) = O- 

VO 70 L = 2eKP 

H(L) = Oc 


88 


20 
25 


40 


50 


30 
60 
70 


1n 


CONTINUE 

MODt = 0 

LE (Ke sb eS) hETURN 

VOSS 0nd) =~ 99 

CALL GTPATH(NC#+#1eNNeK2) 

LF (MODEeEQ. 99) RETURN 

1F (KK GI eG) GU TO 50 

SUBSUNIC PATH 

CALL LINES (NCeNN) 

lr (KK.EQ8.0) GU TO 30 

'HIS PATH HAS KK FORKS 

NFORKS = =<KK 

VI 4U K = JeNFURKS 

LALL GTPATH(NN+1 NF eK2) 

CALL LINtS(NNe VF) 

CUNTINUE 

bU TU 30 

CUNT LNUE 

LF (KK,EQ.2) CALL XXPTS (NCE,NN) 
Le. (K2,NE.&) CALL LINES (NC,NN) 
CALL GTPATH (NC+H1lLeNNeK2) 

LF. (KK,EQ.2) CALL. XXPIS, (NE,NN) 
LF (K2,NE,E) CALL LINES (NC,NN) 
CUNTINUE 

FURMAT(LHL/S9XeL4ATAPE 62 PATH UO /) 
FURMAT (LHO/S/S9X15SHAUIOMALION PAIHNSs) 
KETURN 

cNO 


SUBRUUTINE XAPIS (K1leKe2) 

MEAL MACH 

CUMMON /A/S GAM“4As “MACH 

CUMPLEX AKRAY se XLOFS es SUF XI so seSCrRII 

COMMUN ARRAY (1) 

LOGICAL ILSw 

YJATA LSW /,FALSC./ 

bE. (LS) 4G0 tO-5 

45wW = TRUE. 

WCRLI = (GAYMA=1.)/(GAMYMA41,) + Se/( (GAYMMA4L, ) eMACHEYACH) 
CALL GETHSG (C4MPLX(QCRIT+6.)eSCRII9S) 

CONS! = REAL(SCRIT) 

ARRAT(K24+1) = ARRAY(K2) 

tj =s Ke 

AKRAY(L) = XLUFS(CONUG(CUNST/SSOFALIARRAY(L) DP eARXAV(L41)) 
Gs Let 

Le (beGieK1) GY TO 10 

Me TURN 

LNO 
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KL = CMPLX(LUGS (Sed) eLOGS(4,U)) 

ANG = RAU*ATAN2(AIMAG( XL) +REAL(AIL)) 
CALL CSY™“BL (XlLe4enml1) 

XL = KI¥(1,+-e0S5/CAB8BS(XI)) 

SLZE t= 4659 

CALL CSYMBL (X1.20e-1) 

SEZE = 614 


CONTINUE 
ANG = 0. 
SCF = SF 
OF S21. 


Ir (NRN.GE.G) GO [0 55 

CALL XYAXES ((O0ceOe) eSeSeSedeSFleVe) 
CALL XYAXES ((UeeOe) ce 2eUr5Se5eSFI2 90, ) 
LALL CSYMBL ((-°3.00°3.5) e1 1 TLE +60) 

GU TO 586 

CALL XYAXLS((UceOe) e605 Se 5eSFLe0e?) 
CALL XYAXES ( (Use) e4oetoesSFle Ie?) 
CALL CSYMBL ((#4e095e¢5) e1LTLb +60) 

SF = SCF 

ofzZze = .07 

TF (2PL.6E,1) CALL PLPATH(24MAX0(OeNP) se IPt+2) 
stZe =e le 

Ke TURN 

tNO 


SUBRUUTINE PLPATH(K1+¢K2) 

LUMMON /US PLlaAWeIRPaNKVe NX eNCRebLoMe LPL Te WIDE se Noe NC eNivelLledSeX%A 
COMYMONSGS NI weN2eNGOeNT7TeNSe NG eM1 

AcwINu N6 

dr (K2.eGee) wrRITe (NY_,6U) 

MODe = 1 

LALL GTPATH (NC4H1eNNeKe) 

fe (MODE.ES 99) RETURN 

LF (Ke.Gt.e5) CALL LINES (NCANN) 

NUD Ease 

CALL GTFATH (NC41eNNeK2) 

LF (MODE-EG.99) RETURN 

LF (K2,.Gt.5) CALL LINES (NC,NN) 

4F (K2.61T.2) CALL POINTS (NNe NN) 

1F (K2.E0.5) RETURN 

FE AKL Le. 2) SOTO 25 

Zr (K2,EU.2) WRITE (N4e7U) 

VO 2U MOVE = 3441 

CALL GTPATH (NCHLyNNeKeQ) 

tr (MODE.EG,9I) RETURN 

LF ((K2,.b0,4) UR (K20EQ.7)) CALL PULNTS (JdeNN) 
dr ((K20E0,5) eURe(K206E.8)) CALL LIVES (NCeNN) 
tr (KA,LE.0) GU FO 20 

LALL GCTPATH (NC+1eVNeK2) 

TF ((K2.E0.4) eURe (K20EQ.7)) CALL POINTS (Jd_yNN) 
LF ((0K€2.00,5) -URe(K2eGE.4%)) CALL LINES (VC eNN) 
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CUMMUN /1/ SFeSIZEsANGe XMAXe YMAXeXURVVYURIPGSIZ 
COMMON /ITy TLlite(8) 

JLMtLNSTON CONL(20)+TEMP( 61) 

tXTERNAL CSQKT+SOFXI 

GAM = (GAMMN=1le)S2e 

CUSQ = GAM+ 1¢6/7(MACH*#MACH) 

CUMX = AMIN] (CUSUs3.) 

RAD = 160./PI 


SFI = PLTSIZ/5. 

LALL ORIGIN (Ocee5S) 
KIMIN = (=7.Us74e0)*5F1 
KIMAA = (4.095¢U)*5SF1 


1x = 44 
ek) iar a6 
ASZ = 11.0 
YSZ = 9.0 


fF (NRN.GE,G) GO TO 6 
CALL ORIGIN (20e25e4.0) 
SFI = SFL+t+SFI 

XIMIN = (-2,/7d9-2.0)4SFI 
MLAAXK =H (5, 75.95.05) *Shl 


JX = 4Yy 
LS 60 
XSZ = 7.50 
YSZ = 5.29 


UD LO!) =) LiwiNu 

tMx = FLUAT(INU+leJ)/FLOAI( IND) 

KL = COMAPEMX*CMX/ (1 6+GAM*EEMX¥EMA) 

CALL GETHSG(XLeSeHSOPR) 

CONLCINO+1-U) = REALS) 

CALL CONIUR (IXeUXe INO eXSfYeYS!feSOFXIeCUNL eo FEMPeXlLMIN, KLYAX) 
SF = 1./SF]I 

CALL ORIGIN (7460445) 

1lF (NRN.Y%ELO) FO TO 5 

CALL ORIGIN (620450) 

LCALL CSYMBL (KLAglle-1) 

CALL CSYMBL (ALBelle-l) 

CALL CSYMBL(XLL.129-1) 

LALL CPLUT( X15) 

CALL CPLUT(Xide2) 

CALL CPLUIT (X1LA_g2) 

X10 = XIOFS((UesOe) el -1.90.)) 

CALL CSYMBL (XlJ063e2"1) 

JU S3U J = 1eNdPS 

CALL CSYMBL(BP(J)ele-l ) 

LF (CUCREAL (POLE) -EQe Oe) eANDe (AIMAGIPOLE) eEdeGe)) CPV ID 355 
KRPOLE = REALIPULE) 

AF (CRPOLE,LEsKEAL(XIMIN)) ORS (RPULESGESREAL(XIMAX))) ©O Tu 59 
PLO! THE LOCAILON UF THE POLE 

CALL SYMBOL (PULte del) 

IF (NLUGS.EGe0) GU TY 42 

VI 4¥U J = JyNLUGS 
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i) 


SUBRUYUUTINE ORLGOINOXeY) 
COMMUN SL SFeSIZEeNNGe XMAKe YMAXe XURG YUN PGSIZ 
xXOR* =" X 


TOUR = Y 

Re TURN 

CNTY FRAMEA 

XUR = O- 

YUR = O. 

CALL PLO! (PGS12+00e4-35) 
RETURN 

cND 


SUBROUTINE XYAXES (XeBOUTe1OPe SCF e ANGL) 

X IS THE LCOCAILIN OF THe JRIGIN UN THE AXLS 

BUT IS THE LteNoTH OF THE AXxIS 10 THE LeFT JF THE ORIGIN 
1UP LS THE CeNoTH 10 THe RIGHT OF THE YRIOLN 

CUMPLEX 28becl ene COR 

CUMMUN SL SFeSIZEeANnGeXMAX se YMAXeXUR VA VURIPSSIZ 
UIMENSTION x(2)0 Ye) 

ANGO = ANG 


SFO = SF 

SLZO0 = Si Ze 

VC) =F X0R 

Y.2) = YOR 

ANG: = O'¢ 

SF = 1e 

SU ZES = el 

XUR = X(1)+XOR 

TOR = X(2)+YUR 

é3 CMPLX(-BOT,G.) 


rae CMPLX(TUP+¢0.) 

CUR = (e257 5) 

NC = 16 

tr (ABSCANGL) .NtE.90.) GO 10 10 
VERTICAL YeAXLS 

23 = (O4«e1,)*2Zs5 

210 = 0s 61. eZ 

CUR = (-e6,0e) 

VG 5 

URAW LINt FOR THE AXIS 
CALE CPLOT (2143) 

CALC CPLUT (Z82) 


ES UP 

a = SOT 

N = 14#K+4Ll 

S = eFLOAT(K) *SCF 
n= 207 T0P 


28 = #-FLUAT(K) *H 

ZT = 2B+LOR 

UV 20 I = 1,N 

URAW HATCH MARK 

CALL CSYMBL (fdeNCeH-l) 
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S$ = S4tFLUAT(L=1)*«SCF 

*#*e*INON@ ANSI = SEC WRITEUP AI ENUKH Oe 
ENCUVE (10,230eA) 3 

LABEL AXLS 

CALL CSY%0L (ZI eAs4) 


Z6i= 2Bth 

foal | tear AS ee 

SF = SFO 

SLZE = SLZG 
ANG = ANGU 
XJR = YO1) 
TUR = YC2) 

Ke TURN 

FURKMAT (F4,1) 
ENO 


SUBRUUTINE CONTOR (MeNeNECOXStVoeFNeCUNL eI cWPeXLMLVeATMAA) 

M IS THE NUMBER UF SURDIVISIUNS ION THE X JLRECTILUy 

V IS THe NUMBEX OF SUBDIVISIONS IN Toe Y JLRECI LUN 

Roers THE *IUIH OF THE PLEOTTtO AREA 

Po oke THE HIGHT OF THE PLUTTED AREA 

LABS(MC) 1S JHe NJIMBER JF CUNTOUR CURVES, VNC NeoAILVE UnfAWws A SUX 

AROUND THe PLOTTtU ArcA 

PN IS THe COMPL-X FUNCTIUN WHOSE ABSOLUTE SQARtU 1S ctu 3st rPLUITEDU 
XI“IN AND XIMAX ARE VALUCS OF XI ALT CUAER LEP ANU URPCR ALonT 

CUNL LS AN ARRAY VF COUTUIR LEVELS IN LNCRASTNG VRUEK 

teMP 1S A TEMPURARY ARRY UF LENGIAH 

CUMPLEX FNewoXLMLVeXLMAK 

UIMENSTON CONL(1)+TEMP(1) 

MUXL = (REALOXLMAX) -REALCAIYMIN) )ZFLOAT IM) 


AUXL = (AIMAG(AIMAX)-ALYMAG( KIMIN) D/FLOAT(N) 
JX = XS/FLOAT(M) 

UY = Y¥S/FLUAT(CN) 

NCL = JARS (NC) 

Lats as 8 

VV 10 I = 31,"Pr 

Wt FN(XIMIN4EFLUAT (1-1) *#RuUXxl) 

TeMPC1L) = REAL a) eREAL (wa) + ATMAGLIW) &ALMAO(W) 
Views 8s 

UY ¥UU J = Jen 

Aree” Oe 


W = FN(XLMIN+EC%IPLA(06. eFLUAI (J) KL AL)) 

IM = REAL (CW) eRe AL (ke) + OLYMAGIW) FALMAG(A) 
YPP = YRP#UY 

VV 360 1 = 1eM 

APP = Xb+D0A 

A = FN(XLIMINFC4PLXIFLOAT( 1) eR UXL ge PFLOATI IU) © A DX1)) 
IP = REAL) e Kec AL IW) + ALMAG(W) CALMAG(W) 
1UP = AMAX1(TPe IT 4e TEMP (1) ee TEMP(141)) 

SUT = AMIN] (FPe IT Me TEMP (J) e TEMP (141) ) 

dr (CONL(NCL) eLI.e3UI) GUY IU $50 

YOU 4¥U L1 = LeNCL 


93 


Jr (CONL(L1).61e30T) GO TU 5y 
4o CUNTLNUE 
60. FO 350 
50 YUL S#0ick =-L «Nee 
tre CCUNL(L) -GTefOP) GO TU 350 
N=. LOWERALE ET s=LOWER RIGHI C=UPPEK LEFT D=UPPEK R1LOAT 
TeMP(1L)-CONL(L) 
(EMP CL4+1)"-CONL(L) 


ion 
iowoun 


C TM=eCONL (L) 
wv) PP-CUNL(L) 
FURY ITYP EY LUIKING AT JHE SIGNS AT THE FUUR CURWERS 
LA = ISION(1+A) 
ig = ISIGN(2@ed) 
1C = LSION( Yet) 
LU = ISIGN(8eU) 


LTYP = &t(LBbtldtLC+Iv+1)72 
bU 10 (34U,22U1250 92409290, ZOU e270 2BN se 28UeeTI0 ve eOUI CEU ea C40 ye dUe 
1 22U.340) iIfYr 
SLGN OF A 1S UVIFFERENT FROM BeC ANU VU 

220 Xl = XP 

Tl = YP#VUY4A/S(A-C) 


X2 XP+UX4¥ A/S (4-5) 
Y2r= YP 
bO TU 300 


SIGN OF Jd Is DIFFERENT FRUM AC AND DU 


230 Xl = KPHUXK#e SJ (Nes) 

Vek yp 

X2 = XPP 

T2 = YP+VUTa3/ (2) 

bd TU 300 

SL1GV OF A ANU 8 LS OIFFERENT THAN SLGN OF C ANU JV 
240 Al = XP 

Yl = YP#UY#A/(A-C) 

K2 = APP 

T2 = YP+UTaB/ (3-0) 

60 TU S40U 

SLGN UF C€ 1S UVAFFERENT FRUY Aes ANU UO 
250 Xl = XP 

Yl = YP+UYah/(AN-C) 

Xe = XP4LUXrS CAC C=—D) 

| reson fH OD 


vbU FU 300 
SIGN OF A ANO C ls CIFFERENT THAT SIGN OF 3 AND JU 


260 AL = XP+UK4h/(N-B) 
Lear 
K2 = XP+UK4#eC/(C-0) 
12 = YPP 
bY TU 300 

270 l%p = 1 
GU 1U 340 


SLGN UF VU IS UIFFERENT FRUM Agts ANU C 


280 Al = XP4UX#eC/ (Cov) 
Yl = Yer 
Ke =: APP 
Y2 = PtUY43/ (=U) 
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oe x i | 


AanNaaaann 


300 


340 


350 


360 


400 


CARE IGPEOT COMPLEX UXKLs ¥1.)'\03) 
CALE VERPEOT (CCMPEX(X20V2) 02) 
CONT LNUe 

lEMPC(l) = TM 

XP = XPP 

itiess te 

CONTINUE 

REM CALL) s= TP 

ViPies ¥PP 

CONTINUE 

Re TURN 

CNO 


CUMPLEX FUNCTLUN CSGRT(2s3RANCH) 

CUMPUTES THE CUMPLEX SWUARE RUUT UVF Z ANDO CHUUSES THE BrRANCA 
IHE CUT UF THe SQUARE ROUS IS A STRAIGHT LINE FROM THe 
URIGIN PASSING THROUGH C=#dRANCHD 

UVIMENSION 2(2)+3RANCH(2) 

REAL LMAGZ 

MEA =) 2 bt } 

LMAGZ = Z(2) 

KR = RKREALZ*KREALZ+I4MAGZ¥*IMAGZ 

Bae (Oe 

LF (RZEG.O.) GU 10 4U 

Q = SORT(IS*(SART(R) +FABS(REALZ))) 

AF (REALZ.GE-0.) GU TO 1U 


a” 
@ = e5*e1LMAGZ/KR 
bU TU 20 


x = ©5*#1MAGZ/U 


LF (Q#BRANCH( 1) 4+R*EBRANCH(c)) 3094U040 
CSQKI = CMPLXl-de-k) 

RE TUKN 

CSORI = CMPLX( eR) 

RETURN 

LNY 


SUBROUTINE LEW (AgdeNEUS*NSOLNS es LArlB VET) 

LINCAR EeHUATIUNS SULUTIUNS wRITICN AT AEC COMPUILNG CENICR (vTU) 
SULVt A SYSTEM OF LINEAR CQOUATIUNS OF THE FORY vAAkad AY A AVILIFIED 
GAUSS ELIMINATLON SCHEME 

NEQS = NUMRER UF EQGUATIINS ANU UNANOWNS 

NSOLNS = NUMBER OF VECTUR SOLUTIUNS DtolrReJ 


LA = NUMBER OF ROWS OF A AS VEFINCU BY DIMENSION SIATCMENI cNERY 
13 = NUMBER UF RUWS OF ob AS VEFINCU SY OIMENSIUON STATEMCNE cwlery 
AVE! = UCTERMANANT OF Ay AFTeR eXdl FRUY Lea 


ULMENSION A(LAsIA)+s B( 18018) 
VWSIZ = Neus 
VS8SIZ = NSOLNS 
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10 


20 
30 


40 


50 


60 


70 


80 


90 


NURMALTI2Zce EACH ROW BY IIS LARGERI CLEMENT 
FURM PARTLAL VETERMINANT 

UtT = 1.U 

VO 4¥U I = 1eNS12 

SL1G = Allel) 

Le (NSIZeLE.1) SO 10 130 

YO 10 J = @eNSI2Z 

1F (ABS (BIG) GE.ASS(A(I+4d))) GO 10 10 
316 = Alley) 

CUNT LNUE 

BG = 1.0/381G 

UO 20 J = 1.NS1Z 

A(IeJ) = A(1l+J) €8S 

vO 3U J = 1,N35SI2Z 

Bl(IeJ) = B(IeJ) *8S 

VET = DEI *EIG 

CUNTINJE 

START SYSTEM ReCDUCTIUN 

NUMSYS = NSIZel 

VO 120 I = leNUYSYS 

SCAN FIRST CULUMN OF CURRENT SYSTEM FOR LARGEST CLEMEWT 
CALL THE ROw CUNTAINING THIS ebLcMeNT, ROW VBGRW 
NN = [4] 

51G = Allel) 

N8GRwW = 1 

VU 50 J = NNeNSIZ . 

LE CABS(G LEC) eGe CASS (A(IeL))) GO TO 50 
B1G = Aldel) 

NBGRW = J 

CONTINUE 

B66 = 1.0/s1G 

SWAP ROW [ WITH ROW NBGRW JNLESS L=N35Kq 
tr (NBGRWeEQeL) GO TO 80 

SWAP A=MATRIX KOQWS 

VU 60 J = IyNSLZ2 

'eEMP = AINBGRWeJ) 

A(NBGRWed) = AlIeJ) 

A(IeJ) = EMP 

UcT = -OUcT 

SAAP BeMATRIX KOWS 

V0 70 J = 1,Nd3512 

'EMP = BI(NBGRWeJ) 

BINSORWed) = SII eJ) 

B(IeJd) = TEMP 

CLIMINATe UNKNJUWNS FROM FARST COLUMN OF CURRENT SYSTEM 
UJ 110 K = NNeNSIZ 

CUMPUTE PIVOTAL MJLTLPLIER 

PMULT = -A( Kel) «BG 

1F (PMULI.EQe-Ue-) GO TO 110 

APPLY PMULT TO ALL COLUMNS OF THE CURRENT A-MATRIX ROw 
VU 90 J = NNeNSIZ 

A(KyJ) = PMULT#*A(I +d) +A(KedJ) 

APPLY PMULT TU ALL CULUYNS OF MATRIX 8 
VO 1U0 L = 1eNdSIZ 

Bi(KelL) = PMULI*6(IebL)+6(KeL) 
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AAANAOG 


110 CONTINUE 
120 CONFILNUE 

VU 3ACK SUBSTLIUILON WITH 3g-MATRIX COLUMN = NCOLS 
130 YU 170 NLULS = leNBSIZ 

VO FUR RUW = NROW 

YU 160 1 = 1leN51Z 

NROW = NSIZ+1le-l 

i1tMer = 0-0 

NUMBER OF FREVLOUSLY COMPUTED UNKNOWNS = NAS 

NVXS = NSIZ-NROW 

AKE WE DVING [exe 3ZO0TTOM RKUW 

Ir (NXS.E0.0) GO FO 150 

VO 140 K = 1eNXS 

KK = NS1Z4+1-4 
140 $'tMP = TeMP4+8( KK, NCOLB) FAINRUW SAK) 
159 BI(NROweNCOLS) = (3INROWeNCOLB) -TEMP) SAINRUASNRUW) 

AAVe wE FINISHED ALL ROwS FOR BeMAIRIX COLIYN = NCULS 
160 CONTINUE 

HAVE wE JUST FINISHED WIIld B=MAIKIX COLUMN NCOLB=NSI17 
170 CUNTINUE 

NOw FINISH COMPUTING THE JETERMINANT 

U0 180 IT = 1eNSIZ 
TAo0 VEEP = DE1enchel) 

KETUKN 

CND 


SLOCK DATA 

LUMMON/O/S NisNe2sNGeN7e9N3eNG_1 

LCUMMUN SL/S SF ae SILZEeANG A XMAX eS YMAXVKURG VUR POSITS 

CUMMUN /J/ PO(4U3) eXMING IMINO XFACHTFACINCAtLSZeNELNe TwOe JL NO ACH 
VATA NL aN2eNSeN4eNGeN7T/1L eee 5e4eSeT/ 

eee CHANGE NCH TO 4 FOR 154%3600%*** 

**e*CHANGE LSZ TO 12U OV 18M 36044 

JATA KMINGYMLNG XFAC Ce YFACeNCHeLSZeNLINJ m1Le2emeGeDVe te 25eelUei due 31S 
VATA SF eSLZE ANS se XMAA Se YMAA re XUN eC VURePGSLZ/DeDeel4eUestlesllesuesies 
1 12-00/ 

tND 


**eee*ewR IeUP FUR ECNCUUVEKUCCOUE STAIEXMCNIS#& eee 

cNCUVE/UCCOCDE STATEMENTS 

tVCUVE AS SIMILAR TO A FURMATTEO wRITe STAICMENT 

UeECQUE AS SIMILAR 10 A FORMAIICO RtAU STAIEMENI 

VO TAPES ARE ACTUALLY INVUL VEO SINCe THe IRANSFERS Tant PLACL 
WITALN CURE 

IME FIRST ARGUMENT Is THe NUMdtR VF CHARACTERS IV de KeEAU UR 
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OA OC OGeEe 


aRITFten 

(dE SeCUND ARGUMENT IS THe FORMAL NUMB3ER 

tHE IHIKU ARGUMENT IS THt ADURESS 160 WHICH OR FRIM WHICH THE 
LAOLLISTO WICES Se WRETITEN JR REAY 

A MORE OtTPILeV VESCRIPTLUN CAN Be FOUND LVN THe FULLOWING YaANUAL 
CUNIRUL VATA 6400/65U0/56U0 COMPUTER SYSTC%S FURIRKAN KeFERC NCE 
MANUAL, SECTION 10.6(ENCUSE/UECODE STAIEMENTS), PAGES 15-17 (1969) 


Kee eKWR ET CUP FUR CALCOMP XUUIINES*¥#*« 

STMBUL 

SUBROUTINE TO WRITE MESSAGES ANU SPECIAL CHARACTERS Ow THe 
BRAPHLC UOJIPUT, 


USAGE 
CALL SYMBOL (Xs VeHe TEXT e ANGLE se NLHAR) 
WHERE 
Xa ARt IHe COORDINAIcS ON JIHte GRAPPA PAPER Ur THE FIRKSI 


CHARACTER OF THE YESSAGE, IF NOHAR IS PIUSLITVE, THE 
PEN SIARTS AT THE COWERL ECE T CUORNECRY OFUIAE ERR ST 
CHARACTER, IF NCHAAR IS NEGATIVE, THE PtEV STARTS 
AT THe CeNTER OF JHE FIRSI CHAKACIER,. 
fa 1S THe HEIGHT OF & CHARACTER IN (REAL NUMBER) INCHES, 
CHARACTERS FROM THE BIeOCIAL EWUIVALENCeE [ABLE 
ARE 4/7 AS WIUE AS THEY ARt AlGd, THE srPace 
BETWEEN CHARACTERS IS 2/7 IHAT OF THEIR ALIGHT. 
'HUSe THe TOTAL SPACE REGUIREV FOR N CHARACTERS 
CQOLALS (6/7)*(HI¥*IN) « 
AS THe CHARACTERS Ake STORED IN A 7 BY 4 MATRIX, 
AN INIEGER MULTIPLE OF 207 INCHES SIVES IHE Abts 
RESOLUTLONV. 
IcXxl LS-THE= VESSAGE JSED. IF NCHAR IS- positive URS ZEKO 
(A HYULLERITH CONSTANT |= NH MesSAhot). ON 
LS ONce SYMBOL TAKEN FROM IHE INTEoeER EQVULVALENCE 
'ABLE LF NCHAR IS NEGATIVEe (Ste 32L0w) 
ANGLE 1S THe ANSLE AT wAICH THE LINE OF CHARACTERS LS URAWN 
IN GeoRELS 
VCHARD>O GIVES IAaet INTEGER NUMBER OF CHARATIERS LV Text TU ge 
WRITTIcN 
=0 WRITtS CHARACTERS IN TEX! UNTIL & O CHARACTeR Ls 
ENCOUNTERED 
€O ONLY ONc CHARACTeR IS PRINIEDe THe CHARALI&ER BperiNo 
TAKEN rROM THE INTEGER EQUIVALENCE SYMBOL TARLEL, 
IHE FOLLUYUWING LNTcCGER EQUIVALENCES ARt JUStJ 
INTEGER 1 ORAWS AN OUCTOGON 
INTEGeR 3 OrAw>s A + 
INTEGER 4 ORKAWS AN X 
INTEbeERK 5 DRAWS A DLAMONU SHAPE 
INTEGeR 11 ORAWS A « 
INTEGER 12 DRAWS AN X WITH BARS ABOVE ANI seLUW LT 
INTEGeR 15 ORAWS A = 
INTEGER 1656 DRAWS A VERTICAL LINE 
INTEGER 20 DRAWS AN ARROW POLINIING [TO THe RIGHI 
INTEbER 25 DRAWS AN O 


INTEGER 63 OrAwd A PLUS WLlIH A 30X AKOUNY LI 
POT 
SUBROUTINE TO 4YOVE PEN OW PLUTI cr. 
PURPUSE 
A) CRrATe AN ORIGIN 
5B) PUSITLON [He PEN AND UU THE URAWING 
C) FERMINATE THe JSE OF IAE PLUT PACKAGE, 
PLOT IS CALLEU 3Y ALL ROUTINES THAT VO GRAPHIC wrlT ino. 


JUSANokt 
CAEL PLORG Ke VoelPEN) 
Xv (HE Pen IS MUVED FROM ITS PRESENT COORDINATES 10 


(He COURVDINATES (X,Y) ReLATIVe VO FAE CUXRENT 
ORIGIN. X AND Y ARe REAL NUMBERS MEASURCU LY INCHES, 
LPEN=1 WHILE MOVING TO (AeyY) THe PEN KeMALNS UN UR OFF THE 
PAPEKe AS BEFORE 3SEING CALLED, 
=? THE Pen MARKS THe PAPER wHlLe MOVING FRUY LIS 
FRESENT* POST TON” TOP Oxy) 
=4 THE PEN voOtS NOT “ARK JHE PAPER WAHLLE MUVING FROM 
ITS PRESEN SPOSIFLON T.0 “€X <7) 
="Ler2erd JHE PEN TS MOVED RESPECTIVELY AS? FOR +1042, UR +35 
MOKEQOVER, THE FINAL PUSIILUN JF THE PEN (Xe7Y) 
BeECOMcS THt NEW URIGIN (UereOe)e 
=99y TERMINATES SHE’ PEOT FILES - 1He> UseR BUSI CALL THs 
AT THe ENJQ OF THt PROGRAM IF IHE oRAPHIC UUT?PUI 
1S NOT 10 BE DESITROYeD. 
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This report was prepared as an account of 
Government sponsored work. Neither the 
United States, nor the Administration, 
nor any person acting on behalf of the 
Administration: 


A. Makes any warranty or representation, 
express or implied, with respect to the 
accuracy, completeness, or usefulness of 
the information contained in this report, 
or that the use of any information, 
apparatus, method, or process disclosed 
in this report may not infringe privately 
owned rights; or 


B. Assumes any liabilities with respect to 
the use of, or for damages resulting from 
the use of any information, apparatus, 
method, or process disclosed in this 
report. 


As used in the above, "person acting on behalf 
of the Administration" includes any employee 
or contractor of the Administration, or 
employee of such contractor, to the extent 
that such employee or contractor of the 
Administration, or employee of such contractor 
prepares, disseminates, or provides access to, 
any information pursuant to his employment or 
contract with the Administration, or his 
employment with such contractor. 
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This book may be kept 


FOURTEEN DAYS 


A fine will be charged for each day the book is kept overtime. 
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